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REMARKS 

Claims 1-51 are pending, with claims 1, 12, 29, 31, 33, 36, 42, 46, 48 and 50 being in 
independent form. 

Claims 29, 45, 46, and 50 have been amended to remedy informalities. No new matter 
has been added. 

An Appendix with four (4) exhibits is included with this Amendment. 
Reconsideration and allowance of the above-referenced application are respectfully 
requested in light of the following remarks. 

Rejections under 35 U.S.C. SI 12 First Paragraph (Best Mode') 
Claims 1-51 stand rejected under 35 U.S.C. §112, first paragraph, because the best mode 
contemplated by the inventor has allegedly not been disclosed. The rejection is respectfully 
traversed. 

The Examiner cites passages and phrases m the Specification, such as "m various 
embodiments" and "the order of description should not be construed as to imply that these 
operations are necessarily order dependent," as alleged evidence of concealment of the best 
mode and as examples of language in the Specification that raised questions in the Examiner's 
mind as to whether the best mode has been set forth. Applicant respectfully asserts that such 
statements in Applicant's specification are not evidence of concealment. 

The best mode inquiry involves two parts. First, a determination is made as to whether, 
at the time of the application, the inventor knew of a mode of practicing the claimed invention 
that the inventor considered to be better than any other. Second, if it is determined that the 
inventor knew of a mode that he considered better than any other, a determination is made as to 
whether the disclosure is adequate to enable one skilled in the art to practice the best mode. See 
MPEP 2165.03; Chemcast Crop. v. Arco Industries Corp., 913 F.2d 923 (Fed. Cir. 1990); Eli 
Lilly & Co. V. Barr Laboratories, Inc., 251 F.3d 955 (Fed. Cir. 2001). Satisfaction of the best 
mode requirement is a question of fact subject to the clearly erroneous standard, but the absence 
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of evidence to support the basic tenets of a best mode violation satisfies that standard. See 
Bruning v. Hirose, 161 F.3d 681, 687 (Fed. Cir. 1998). Because the best mode inquiry often 
requires information that is unavailable to examiners during ex parte examination, examiners are 
instructed to assume that the best mode is disclosed unless evidence is presented that is 
inconsistent with that assumption. See MPEP 2165.03. 

Firstly, the examiner's evidence does not show that at the time of the application, the 
inventor knew of a mode of practicing the claimed invention that the inventor considered to be 
better than any other. The Applicant respectfully asserts that the cited language fi-om the 
specification merely suggests the possibility of other embodiments beyond those explicitly 
disclosed in the specification. Such language is not inconsistent with a proper disclosure of the 
best mode in the specification. Secondly, assuming arguendo that the asserted evidence satisfies 
the first prong of the best mode inquiry (which it does not), the evidence does not show that the 
disclosure is inadequate to enable one skilled in the art to practice the best mode. The Examiner 
has not provided any evidence that shows that the quality of the disclosure in the instant 
application is so poor as to effectively result in conceahnent. See MPEP 2165.04. Moreover, 
one cannot conclude that there is a concealment of the best mode merely from a possibility that 
there may be other embodiments beyond those explicitly disclosed in the specification. 
Accordingly, the Examiner has not met the requisite legal burden for showing the lack of best 
mode. 



Rejections under 35 U.S.C. SI 12 First Paragraph (Enablement^ 
Claims 1, 11,21, 23, 25, 30-33, 36, and 47-50 stand rejected under 35 U.S.C. §1 12, first 
paragraph for allegedly lacking enablement. The rejection is respectfully traversed. 

Any analysis of whether a particular claim is supported by the disclosure in an 
application requires a determination of whether that disclosure, when filed, contained siLfficient 
information regarding the subject matter of the claims as to enable one skilled in the pertinent art 
to make and use the claimed invention, ^ee MPEP 2164.01. The standard for determining 
whether title specification meets the enablement requirement was cast in the Supreme Court 
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decision of Mineral Separation v. Hyde, 242 U.S. 261, 270 (1916) which postured the question: 
is the experimentation needed to practice the invention undue or unreasonable? That standard is 
still the one to be appUed. In re Wands, 858 F.2d 731, 737, 8 USPQ2d 1400, 1404 (Fed. Cir. 
1988). 

In order to make a rejection, the examiner has the initial burden to establish a reasonable 
basis to question the enablement provided for the claimed invention. According to In re Wright, 
999 F.2d 1557, 1562, 27 USPQ2d 1510, 1513 (Fed. Cir. 1993), the examiner must provide a 
reasonable explanation as to why the scope of protection provided by a claim is not adequately 
enabled by the disclosure. A specification which contains a teaching of the manner and process 
of making and using the invention in terms which correspond in scope to those used in 
describmg and defining the subject matter sought to be patented must be taken as in compliance 
with the enablement requirement unless there is a reason to doubt the objective truth of the 
statements contained therein which must be reUed on for enablement. See id. ; In re Marzocchi, 
439 F.2d 220, 223-224 (CCPA 1971); In re Wright, 999 F.2d at 1 561-62. 

Claims 1 and 36. Claims 1 and 36 recite the feature of a "generation technique." To the 
extent that the Examiner's argument is inteUigible, it appears that the Examiner has alleged that 
this feature is uncles or undefined in the specification. As a prehminary matter, the Examiner 
has not provided a shred of evidence supporting the position that this feature is "unclear" or 
"undefined." In fact, the very evidence that the Examiner proffers shows just the opposite. The 
specification describes a "generation technique" and provides examples of such, namely a 
triangular technique and a quadrilateral technique. See Specification, p. 9 lines 13-29. The 
specification also describes how to perform the given examples of generation techniques. See 
Specification, p. 10 line 7 - p. 14 line 18. 

Secondly, the Examiner provides no evidence that there is a reason to doubt the objective 
truth of the statements in the specification. One of ordinary skill in the art who reads the 
specification will understand what AppHcant's claimed "generation technique" is and how to 
make and use such. Thus, the Examiner has not met the legal burden for showing a lack of 
enablement for this feature. 
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Claims 1. 25. 32, 33. 36, 49. and 50. Claims 1, 25, 32, 33, 36, 49, and 50 recite the 
features of "selecting," "applying," or "determining" within a computing environment. It 
appears that the Examiner has alleged that these features are unclear or undefined in the 
specification. The Examiner has not provided evidence supporting the position that this feature 
is "unclear" or "tmdefmed." Furthermore, the specification describes examples of the selecting, 
applying, and determining operations. See, e.g.. Specification, p. 10, lines 1 - 5; p. 16 lines 9 - 
12. The Examiner provides no evidence that there is a reason to doubt the objective truth of the 
statements in the specification. Thus, the Examiner has not met the legal burden for showing a 
lack of enablement for this feature. 

Furthermore, these features were well known to those of ordinary skill in the art as of the 
filing date. A specification need not teach, and preferably omits, what is well known in the art. 
See MPEP 2164.01. The Examiner has not alleged that such operations are not within the skills 
of a person having ordinary skill in the art. One of ordinary skill in the art who reads the 
specification will understand how to make and use the selecting, applying, and determining 
features. Because such operations are well known in the art, fijrther description in the 
specification is not necessary for purposes of enablement. 

Claims 11 and 23. Claims 1 1 and 23 recite the feature of a "non-regularized" Boolean 
operation. It appears that the Examiner has alleged that this feature is unclear or undefined in the 
specification. The Examiner has not provided evidence supporting the position that this feature is 
"unclear" or "imdefined." 

Furthermore, "non-regularized" Boolean operations were well-known in the art as of the 
filing date. See, e.g., Rossignac, "Issues on Features-Based Editing and Interrogation of Solid 
Models," Comput. & Graphics, Vol. 14, No. 2, pp. 149-172 (1990) (Exhibit A); Rossignac and 
Requicha, "Constructive Non-Regularized Geometry," Computer-Aided Design, Vol. 23, No. 1, 
pp. 21-32 (1991) (Exhibit B). These references describe non-regularized Boolean operations 
well before the filing date of the instant application. A specification need not teach, and 
preferably omits, what is well known in the art. See MPEP 2164.01. Because "non-regularized" 
Boolean operations are well known in the art, fiirther description in the specification is not 
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necessary for purposes of enablement. Thus, the Examiner has not met the legal burden for 
showing a lack of enablement for this feature. 

Claims 21. 30. 31. 47. and 48. Claims 21, 30, 31, 47, and 48 recite the feature of a "wire 
body." The Examiner asserts that the term "wire body" is not disclosed in the specification. It 
appears that the Examiner has alleged that this feature is unclear or imdefined in the 
specification. The Examiner has not provided evidence supporting the position that this feature 
is "unclear" or "undefined." Also, an enablement rejection that is made merely because a term is 
recited in the claims but not in the specification is improper. The claimed limitation in and of 
itself may enable one of ordinary skill in the art to make and use the claim containing the 
limitation. See MPEP 2164. The Examiner has not alleged that the term "wire body," as recited 
in the claims, is not enabling in and of itself Thus, the Examiner has not met the legal burden 
for showing a lack of enablement for this feature. 

Furthermore, Apphcant submits that "wire body" was well-known in the art as of the 
filing date. See Lambda Research Corp., Illuminations newsletter, vol. 2, issue 1, Jan. 2000, 
pp. 1-4 (Exhibit C); Schild et al, "Open Data Exchange with HP PE/SolidDesigner," Hewlett- 
Packard Journal, October 1995, pp. 35-50 (Exhibit D). The references describe a wire body as 
synonymous with a wire frame. Because "wire body" is well known in the art, fiirther 
description in the specification is not necessary for purposes of enablement. 

Claims 25. 32. 33, 49. and 50. Claims 25, 32, 33, 49, and 50 recite in part, "valid ones of 
said faces." The Examiner asserts that the term "valid" is imclear because the Specification 
discusses "invalid" faces but not 'Valid" faces. The Examiner has not provided evidence 
supporting the position that this feature is "unclear" or "undefined." Also, an enablement 
rejection that is made merely because a term is recited in the claims but not in the specification is 
improper. The claimed limitation in and of itself may enable one of ordinary skill in the art to 
make and use the claim containing the limitation. See MPEP 2164. The Examiner has not 
alleged that the term "valid," as recited in the claims, is not enabling in and of itself Thus, the 
Examiner has not met the legal burden for showing a lack of enablement for this feature. 
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Furthermore, Applicant submits that the term "valid" is sufficiently clear in the 
specificaticHi in the specification. One of ordinary skill in the art will reasonably understand that 
"valid" is the opposite of "invalid." The specification describes various factors for determining 
whether a face is invalid. See Specification, p. 16, lines 13-21. One of orduiary skill in the art 
will reasonably understand that a "valid" face is a face that is determined as not invalid, based on 
the factors described in the specification. Thus, the term "valid" is sufficiently clear. 

Rejection under 35 U.S.C. §1 12 second paragraph 
Claims 1 and 36 stand rejected tmder 35 U.S.C. §112, second paragraph for allegedly 
being indefinite. The rejection is respectfiilly traversed. 

The Examiner alleges that the feature "generation technique" in claims 1 and 36 is 
unclear and that it is difficult to determine a definition for this term because the Specification 
read "[i]n various embodiments, the generation techniques being considered may include, but are 
not limited to . . . ." AppUcant submits that the definition of "generation technique" is definite in 
hght of the specification. The Specification describes a "generation technique to generate a data 
representation to model a fillet weld bead," provides examples of generation techniques, and 
describes how to perform such. See Specification, p. 9 line 13 - p. 14 line 18. In Ught of the 
Specification, and as addressed above, one of ordinary skill in the art will reasonably discern that 
"generation technique" refers to a technique used to generate a data representation to model 
geometiies, such as a data representation to model a fillet weld bead. Accordingly, tiiis rejection 
should be withdrawn. 

Rejections under 35 U.S.C. §103fa) 

Rejection over Subrahmanyam in view of CAI 

Claims 1-4, 7-1 1, 12, 15-39, 42, and 45-51 stand rejected under 35 U.S.C. §103(a) as 
allegedly being unpatentable over Subrahmanyam et al. (erroneously written as "Subrashekar" in 
the Office Action), "Feature Attributes and Their Role in Product Modeling," Sohd Modeling 
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'95, Salt Lake City, Utah, (1995) (hereafter referred to as "Subrahmanyam"), and further in view 
of Computer & Automation Institute, "PROARC, No. 7831, CAD-Based Programming System 
for Arc Welding Robots in One-Off Production Runs," ESPRIT, (Jan. 16, 2001) (hereafter 
referred to as "CAI"). The rejection is respectfiilly traversed. 

Claim 1. Claim 1 recites in part, "selecting ... a generation technique based at least in 
part on the result of said examining; and applying ... the selected generation technique to 
generate a data representation of the fillet weld bead." 

The Examiner asserted that Subrahmanyam teaches a method to generate a solid model 
that includes selecting a generation technique based at least in part on the result of said 
examining, and applying the selected generation technique to generate a data representation of 
the solid model. For this proposition, the Examiner relies on sections 4.2 - 4.4 and Table 1 of 
Subrahmanyam. In the cited portions of Subrahmanyam, operations are performed on a solid 
model and then attributes are processed after the operation. How an attribute is processed may 
be based on the operation performed and the attribute in question. For example, Subrahmanyam 
shows that a face normal vector attribute may or may not change, depending on the operation. 
However, there is no teaching or suggestion that a data representation of a shape is generated by 
processing attributes. Subrahmanyam merely shows how the attributes of a shape may be 
affected by operations on the shape. On the other hand, claim 1 recites selecting a generation 
technique and applying the generation technique to generate a data representation of a fillet weld 
bead. Therefore, Subrahmanyam does not teach or suggest this feature of claim 1. Accordingly, 
claim 1 is in condition for allowance for at least this reason. 

Claims 2-4. 7-11. Claims 2-4 and 7-1 1 are dependent from claim 1 and are allowable for 
at least the reasons set forth above. 

Claim 12. Claim 12 recites in part, "generating ... a tool based at least in part on the 
constructed profile." The Examiner asserts that Subrahmanyam teaches "generating with the 
computer environment, a tool based at least in part on a profile," citing p. 117, col. 2, line 5 and 
Figure 4 (erroneously written as Figure 5 in the Office Action) of Subrahmanyam. The cited 
portion of Subrahmanyam shows a starting block with attributes, a slot with attributes, and a hole 
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with attributes. There is no teaching or suggestion in the cited portion of Subrahmanyam of 
generating a tool based on a constructed profile. Assuming for the sake of argument that 
attributes are profiles, the cited portion of Subrahmanyam shows that the starting block, the slot, 
and the hole have attributes; there is no teaching that they are generated based on the attributes. 

The cited portion of Subrahmanyam also shows subtracting the slot and the hole fi-om the 
starting block. The Examiner relies on Subrahmanyam 's showing of subtracting of the slot and 
the hole fi-om the starting block as teaching generating a tool and conditionally trimming the tool. 
It appears to Applicant that the Examiner is treating the starting block as the "tool," because the 
Examiner is citing the subtracting as teaching trimming the tool and the object being "trimmed" 
in Subrahmanyam is the starting block. However, there is no explicit showing in Subrahmanyam 
of how the starting block is generated. Furthermore, as describe above, Subrahmanyam merely 
shows that the starting block has attributes. It does not teach or suggest generating the starting 
block based on attributes. Thus, Subrahmanyam also does not disclose "generating ... a tool 
based at least in part on a constructed profile." Accordingly, claim 12 is in condition for 
allowance. 

Claims 15-28. Claims 15-28 are dependent from claim 12 and are allowable for at least 
the reasons set forth above. 

Claim 29. Claim 29 recites in part, "locating . . . one or more bodies referred to by a 
plurality effaces of a plurality of components." The Examiner asserts that Subrahmanyam 
teaches locating one or more bodies referred to by the faces. However, the cited portions of 
Subrahmanyam (Subrahmanyam, Section 3.7 lines 1-6, 10-11; Figure 4) make no mention of 
locating bodies referred to by faces. Subrahmanyam Section 3.7 shows attributes representing 
relationships between entities and applying a parallel attribute to a pair of faces, but not locating 
bodies referred to by faces. An attribute is defined in Subrahmanyam as "a characteristic quality 
or property which associates meaning to an entity, significant to a particular stage in the life 
cycle of a product," and gives as examples of attributes the color of a face, type of thread, type of 
feature, or relationships between two faces. See Subrahmanyam, Section 1 (Introduction). This 
definition and the given examples suggest that attributes are not the same as bodies. Nor is the 
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locating of bodies apparent from Figure 4 of Subrahmanyam. Figure 4 of Subrahmanyam shows 
a design view that includes a starting block with attributes, a slot with attributes, and a hole with 
attributes. There is no teaching in Subrahmanyam Figure 4 that the starting block, the slot, and 
the hole are located or that they are referred to by faces. Thus, the Examiner has not shown that 
Subrahmanyam teaches or suggests "locating . . . one or more bodies referred to by a plurality of 
faces of a plurality of components." Accordingly, claim 29 is in a condition for allowance. 

Claim 30. Claim 30 is dependent from claim 29 and is allowable for at least the reasons 
set forth above. 

Claim 31. Claim 3 1 recites in part, "conditionally forming ... a wire body, if, data 
representations of more than one edge are replicated." The Examiner asserted that 
Subrahmanyam teaches conditionally forming a wire body if data representations of more than 
one edge are replicated (Subrahmanyam, Section 4.4, lines 12-17; Figure 6(b)). However, the 
cited portion of Subrahmanyam does not teach conditionally forming a wire body. The cited 
portions of Subrahmanyam show copying attributes of faces. Copying attributes of faces does 
not form a wire body in and of itself Nor is it apparent that the shape depicted in Figure 6(b) of 
Subrahmanyam is a wire body. Thus, the cited portions of Subrahmany^ does not teach or 
suggest forming a wire body. Accordingly, claim 31 is in a condition for allowance. 

Claim 32 . Claim 32 is dependent from claim 31 and is allowable for at least the reasons 
set forth above. 

Claim 33. Claim 33 recites in part, "copying and extending ... the selected valid ones of 
said faces into bodies." The Examiner asserts that Subrahmanyam teaches copying and 
extending, within the computing enviroiunent, the selected valid ones of said faces into bodies. 
However, the cited portions of Subrahmanyam (Section 4.4, lines 12-15; page 121, "Merging 
Operations," ^11 1-4) do not teach or suggest extending a face into a body. Subrahmanyam shows 
handling of attributes, such as copying or deleting attributes, between faces after a split or a 
merge. Such handling of attributes does not extend a face into a body. A face is not extended 
into a body merely because its attributes may have been changed by the copying or deleting of 
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attributes. Thus, Subrahmanyzun does not teach "copying and extending ... the selected vahd 
ones of said faces into bodies." Accordingly, claim 33 is in a condition for allowance. 

Claim 34-35. Claims 34-35 are dependent from claim 33 and are allowable for at least 
the reasons set forth above. 

Claim 36. Claim 36 is directed to an apparatus having a storage medium. The storage 
medium has stored thereon programming instructions designed to enable the apparatus to, in 
part, "select a generation technique based at least in part on the result of said examination; and 
applying the selected generation technique to generate a data representation of a fillet weld bead 
of the fillet welding operation." As set forth above with respect to claim 1, Subrahmanyam has 
no teaching or suggestion of selecting a generation technique and applying the selected technique 
to generate the shape. Accordingly, claim 36 is in a condition for allowance. 

Claims 37-39. Claims 37-39 are dependent from claim 36 and are allowable for at least 
the reasons set forth above. 

Claim 42. Claim 42 is directed to an apparatus having a storage medium. The storage 
medium has stored thereon programming instructions designed to enable the apparatus to, in 
part, "generate a tool based at least in part on the constincted profile." As set forth above with 
respect to claim 12, Subrahmanyam does not teach or suggest generating a tool based at least in 
part on the constructed profile. Accordingly, claim 42 is in a condition for allowance. 

Claim 45. Claim 45 is dependent fixjm claim 42 and is allowable for at least the reasons 
set forth above. 

Claim 46. Claim 46 is directed to an apparatus having a storage medium. The storage 
medium has stored thereon programming instructions designed to enable the apparatus to, in 
part, "locate one or more bodies referred to by a plurality of faces of a plurality of components." 
As set forth above with respect to claim 29, Subrahmanyam does not teach or suggest locating 
bodies referred to by faces. Accordmgly, claim 46 is in a condition for allowance. 

Claim 47. Claim 47 is dependent &om claim 46 and is allowable for at least the reasons 
set forth above. 
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Claim 48. Claim 48 is directed to an apparatus having a storage medium. The storage 
medium has stored thereon programming instructions designed to enable the apparatus to, in 
part, "conditionally form a wire body, if, data representations of more than one edge are 
rephcated." As set forth above with respect to claim 31, Subrahmanyam does not teach or 
suggest conditionally forming a wire body. Accordingly, claim 48 is in a condition for 
allowance. 

Claim 49. Claim 49 is dependent from claim 48 and is allowable for at least the reasons 
set forth above. 

Claim 50. Claim 48 is directed to an apparatus having a storage medium. The storage 
medium has stored thereon programming instructions designed to enable the apparatus to, in 
part, "copy and extend the selected valid ones of said faces into bodies." As set forth above with 
respect to claim 33, Subrahmanyam does not teach or suggest extending faces into bodies. 
Accordingly, claim 50 is in a condition for allowance. 

Claim 51. Claim 51 is dependent from claim 50 and is allowable for at least the reasons 
set forth above. 

Rejection over Subrahmanyam in view of CAI and Wang 

Claims 5, 6, 13, 14, 40, 41, 43 and 44 stand rejected under 35 U.S.C. §103(a) as allegedly 
being unpatentable over Subrahmanyam and CAI as applied to claims, 1, 12, 36 and 42 above, 
and Wang, et al., "The Design and Fabrication of Welded Tubular Joints Using Solid Modeling 
Techniques," 2nd ACM SoUd Modeling, 1993, hereafter referred to as "Wang." The rejections 
are respectfully traversed. 

Claims 5 and 6. Claims 5 and 6 are dependent from claim 1, and are in a condition for 
allowance for at least the reasons set forth above. Additionally, the cited portions of Wang do 
not teach or suggest a triangular or quadrilateral profile. Wang shows modeling of tubular joint 
welds. There is no teaching or suggestion in the cited portions of Wang (p. 430, last paragraph 
of Introduction; Section 2.4; Figure 8b) of a triangular or quadrilateral profile. There is no 
teaching or suggestion that a profile is of a particular shape. Furthermore, there is no indication 
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that Figure 8b of Wang depicts a triangular or quadrilateral profile. Thus, there is no teaching or 
suggestion of a triangular or quadrilateral profile in Wang. Accordingly, claims 5 and 6 are in a 
condition for allowance. 

Claims 13 and 14. Claims 13 and 14 are dependent firom claim 12, and are in a condition 
for allowance for at least the reasons set forth above. Claims 13 and 14 also recite in part, "a 
triangular profile" and "a quadrilateral profile," respectively. As set forth above with respect to 
claims 5 and 6, the cited portions of Wang do not teach or suggest triangular and quadrilateral 
profiles. Accordingly, claims 13 and 14 are in a condition for allowance. 

Claims 40 and 41. Claims 40 and 41 are dependent fix)m claim 36, and are in a condition 
for allowance for at least the reasons set forth above. Claims 40 and 41 also recite in part, "a 
triangular profile" and "a quadrilateral profile," respectively. As set forth above with respect to 
claims 5 and 6, the cited portions of Wang do not teach or suggest triangular and quadrilateral 
profiles. Accorduigly, claims 40 and 41 are in a condition for allowance. 

Claims 43 and 44. Claims 43 and 44 are dependent fi-om claim 42, and are in a condition 
for allowance for at least the reasons set forth above. Claims 43 and 44 also recite in part, "a 
triangular profile" and "a quadrilateral profile," respectively. As set forth above with respect to 
claims 5 and 6, the cited portions of Wang do not teach or suggest triangular and quadrilateral 
profiles. Accordingly, claims 43 and 44 are in a condition for allowance. 

Conclusion 

For the foregoing reasons. Applicant submits that all the claims are in condition for 
allowance. 

By responding in the foregoing remarks only to particular positions taken by the 
Examiner, Applicant does not acquiesce with other positions that have not been explicitly 
addressed. In addition, AppUcant's arguments for the patentability of a claim should not be 
understood as implying or conceding that no other reasons for the patentability of that claim 
exist. 
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Exhibit A 

Rossignac, "Issues on Features-Based Editing and Interrogation of Solid Models," Comput. & 
Graphics, Vol. 14, No. 2, pp. 149-172 (1990). 



Features and Geometric Reasoning 



ISSUES ON FEATURE-BASED EDITING AND 
INIERROOATION OF SOLED MODELS 



IBM, Research Division, T. J. Watson Research Center, P.O. Box 704, Yorktown Heights, NY 10598 

Abstract— Operations that create additive or subtractive volume features, such as bosses or slots, simplify 
the computer aided design of mechanical parts. Surface features, whether extracted automatically or selected 
interactively, group fimctionally related boundary elements, and thus provide an expedient intetfece between 
CAD systems and atmlysis or manufacturing applications. Despite much inogress in CAD, design remains 
an iterative process and involves error-prone modifications of previous solutions. Features should in principle 
offer a higih level vocabulary for characterizing errors and for specifyinghowth^ should be corrected. This 
paper points out the semantic ambiguities of simplistic feature-based commands for editing models. It 
recommends procedural models for editing volume features, and corrective volumes for editing surface 
features. It shows how space decomposition techniques and CSG expressions based on active zones reduce 
the cost of executing an editing command. Error detection may be automated by supporting intentional 
features, which correspond to the desired characteristics of the model, and by endowing them with domain 
dependent validity criteria ejqiressed in terms of associated geometric elements. The paper demonstrates 
that validity may be tested by simply interrogating a mixed-dimensional geometric structure which is used 
to reinesem not only the model, but also the interaction; between the geometric elements associated with 
intentional features. 



1. INTRODUCTION 

Solid modelling improves the efficiency of the design 
process for manufactured parts by supporting the geo- 
metric representations of these parts. It provides 
graphic feedback to the designer and offers interfaces 
to some analysis applications. Despite recent Progress 
in computer hardware and in interactive graphics, 
geometric design of three-dimensional shapes remains 
a complex and time consuming task. Since geometric 
representations of complex mechanical parts tend to 
be vettx)se, various abstractions, globally calle4 geo- 
melrufeatwes, are often used to characterize certain 
types of shapes or to refer to Portions of a part model 
that may be important to the designer or to an appli- 
cation program. 

Features are used in Computer Aided Design and 
Manufecturing for a variety of purposes: 

• Geometric features provide a concise description of 
the parts characteristics[ l]and thus simplify group 
technology and process planning. They also fecilitale 
the communication between designers and solid 
modelling systems. 

• Features provide a mechanism fijr attaching product 
or manufacturing information and various attributes 
to specific parts of a geometric model (see fijr ex- 
ample [2] for attaching tolerance information to 
features). 

• The properties intuitively associated with common 
feature types define many convenient shape altering 
operations [3] that attempt to create features of 
specified types and dimensions. 

• Access to the geometric elements that compose a 
feature simplifiesthe interrogation of shape by pro- 
viding a naming scheme for sets of boundary ele- 
ments [2], a convenient vocabulary for expressing 
relevant dimensions and positions [4] and for for- 



mulating validity checks that assess the compliance 
of the model with the designer'sintent. 
• Expressing and performing engineering changes or 
simply corrections of design errors may be eased by 
using geometric features [5] 

This paper focuses on the last two issues, namely 
the use of geometric features to automate— or at least 
simplify — the editing of the solid model and the 
checking of the model's compliance with fimctional 
requirements (validify conditions). 

Most CAD models of 3D manufactured parts are 
created by combining and incrementally modifying 
simple models. These combinations and modifications 
are often conveniently expressed in terms of Boolean 
operations. The primitive shapes from which the mod- 
els are constructed are often restricted to arbitrarily 
positioned and sized solid primitives (blocks, cylin- 
ders, spheres ...), generic volume features (holes, slots, 
bosses ..,), and linear or circular extrusions of 2D re- 
gions. The resulting part models can thus be repre- 
sented by a CSG (Constructive Solid Geometry) 
tree [6], which leads to certain algorithmic advantages 
(see [7-9] for examples) and to an obvious archival 
conciseness. 

CSG expressions may be complex, and the end-user 
often prefers to interact with a boundary model, which 
is algorithmically derived from CSG [10] and contains 
the list of feces and their adjacency graph [11], There- 
fore, it is important to develop techniques for inter- 
actively specifying validity conditions and modifica- 
tions in terms ofboundary elements (faces, edges, and 
their incidence graphs) rather than in terms of CSG. 
Domain-dependent features provide a particularly 
convenient vocabulary for accessing relevant groups 
of boundary elements. On the other hand, direct 
boundaiy editing is error-prone, and editingoperations, 
even if specified in terms ofboundary information, 
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should be translated into mathematically well-defined 
(nonambiguous) operations, such as Boolean set op- 
erations or global rounding and filleting opera- 
tions[ 12]. Besides, it is important to maintain a CSG 
representation of the model and to express validity 
conditions in terms of CSG so that the model can be 
parameterized, easily edited, and reused. 

This paper studies the translation process, which 
takes validity conditions or model modifications ex- . 
pressed in terms of features (and thus of boundary 
elements) and performs the appropriate model mod- 
ifications using CSG operations. Due to a lack of for- 
malism ofthe semantics of feature-based specification, 
automatic translation remains a challenging research 
goal (some pitfalls of a "naive" translation process are 
pointed out in Section 3). Several new or recently de- 
veloped techniques are discussed, \^ich do not always 
provide the correct translation, but at least increase 
the designer's vocabulary or can automatically generate 
a tentative solution, which may have to be finther ad- 
justed by the designer. Furthermore, the paper ad- 
dresses flie issue of efficiently performing the model 
modifications or the validity tests by using informa- 
tionally rich geometric structures and properties of 
Boolean expressions. 

The paper is organized as follows: 

• The basic concepts and terminology are introduced 
in Section 2. 

• Section 3 points out some of the limitations of a 
straightforward use of features for editing and inter- 
rogating solid models. This section focuses on con- 
cepts and techniques rather than on their historical 
evolution or implementation, (For a more formal 
survey ofthe literature on features, the reader should 
refer to [13, 14].) 

• The importance of procedural models for capturing 
the designer'sintent into a flexible parameterized 
sequenceis emphasized in Section 4. To edit a feature 
explicitly created by an operation, it may be simpler 
to "ask" the feature what operation created it, change 
the parameters of that operation, and reexecute the 
entire sequence. 

• Reexecuting the entire sequence amounts to evalu- 
ating the boundary of a CSG representation, and 
may be very costly for large CSG models. A new 
approach that reduces the reevaluation cost is pre- 
sented in Section 5. It derives a CSG expression for 
the regions that must be added to, or subtracted fi-om, 
the solid model. Furthermore, Section 5 also presents 
a recently developed mixed-dimensional geometric 
representation called SGC (Selective Geometric 
Complex). Algorithms for SGCs generate a subdi- 
vision of space imposed by the features. This sub- 
division can be used to reduce considerably the 
amount of geometric calculations and of logic 
expression evaluations necessary to perform the fea- 
ture-editing operations. 

• Some features do not correspond to a single opera- 
tion, and it may be too complicated to identify all 
the operations in the sequence that must be edited 



in order to rectify an "invalid" feature. Section 6 
demonstrates on some simple examples how correc- 
tive volumes, obtained by extruding feature faces, 
can be used to perform simple feature alterations 
without reexecuting the design sequence. In more 
complicated situations, these corrective volumes 
must be trimmed before they can be added to— a 
subtracted from — the model of the part. Without 
the trimming step, side-effects may appear, especially 
when several features interact or when compound 
features incrementally created by successive oper- 
ations, are edited. Trimming is best performed using 
Boolean operations, but producing a timming CSG 
expression may prove difficult and remains the de- 
singer's responsibility. 
• Section 7 addresses the problem of feature validity. 
Specifically, it shows how features may be efficiently 
tested by interrogating the corresponding SGC rep- 
resentation. 

2. BASIC CONCEPTS AND TERMINOLOGY 

This section clarifies the distinction between inten- 
tional features and their geometric embodiment, and 
between volume features and surface features. It also 
introduces the CSG notation used in this paper. 

2.1. Intentional features and their geometric embod- 
iment 

A distinction should be made between geometric 
features and intentionalfeatures. A geometric feature 
is a collection of geometric elements (for example, feces 
or volumes) that form a subset of Uie part'sinterior, 
boundary, and/ or complement. An intentional fea- 
ture [15] is an abstraction for accessing groups of geo- 
metric elements and for associating with them a fype 
and consequently certain properties defined for all tiie 
features of this particular type. For example, an inten- 
tional feature oftype slot may be associated with apart 
model. This association indicates that the designer in- 
tends to have a slot in the model, i.e., a void bounded 
on three sides by faces ofthe model. The intentional 
feature contains references to the corresponding feces. 
However, due to model manipulations, the referred 
faces may have been modified or even deleted from 
the modelsboundary. Whatever remains ofthem and 
of the associated void constitutes a geometric feature 
that may no longer exhibit the properties associated 
with a slot, 

Inconsistencies between intentional features and the 
actual geometry of the part are avoided by treating 
intentional features only as hints and by relating them 
to geometric elements through collections of uneval- 
uated references. It is acceptable that some, or all, of 
these references do not correspond to any geometiic 
element ofthe model' sboundary at some particular 
stage during the design process. Even if all geometric 
elements referenced in an intentional feature are pres- 
ent in a geometric model, their shapes and positions 
with respect to the rest ofthe model need not comply 
with tiie characteristics usually associated with the 
particular feature type. For example, an intentional 
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feature of type cylindrical hole could be associated with 
geometric elements (faces) that have been removed 
ftom the model'sboundary by some Boolean opera- 
tion, and thus do not correspond to a 'Valid" hole. In 
such situations, the intentional feature is said to be 
invalid, but should not be discarded, because the de- 
signer may have produced (intentionally or not) tem- 
poral situations, or instances of the model, where 
many previously defined intentional features are in- 
valid. The overall validity may later be restored by 
repositioning a subsolid or adjusting a parameter. The 
d^igner should not be required to redefine all inten- 
tional features that went through an invalid transition 

Furthermore, feature validity is very subjective and 
in fed depends on the role the featureplays with respect 
to a particular application. To take a simplified ex- 
ample, a cylindrical hole is a valid "detail feature" to 
be discarded for analysis purposes only if its radius is 
sufficiently small; on the other hand, it is a valid 
"manufacturing feature" for process planning only if 
it is empty and accessible.Feature validity criteria may 
be expressed in terms of validity rules, which are logical 
predicates defined in terms of the referenced geometric 
elements and of their existence, shape, and relation to 
other geometric elements of the model. Evaluating the 
model's geometric references is therefore necessary to 
establish the validity of an intentional feature with re- 
qjectto any one of the instances (or stages) through 
wliidi a solid model evolvesduring the design process. 
Consequently, the interrogation of invalid features 
plays an essential role in correcting design errors [4]. 
Typically, the presence of an intentional feature of a 
certain type, valid or not, implies some intention that 
Ifae designer has regarding the functionality of some 
portion of the part. Thus, intentional features may 
provide important hints for model alterations and 
manufacturing process planning. 

2.2. Constructive Solid Geometry (CSG) 

CSG (Constructive Solid Geometry)[ 6] refers to an 
unevaluated representation scheme for solids obtained 
by combining, in Boolean expressions, simpleprimitive 
shapes of arbitrary dimensions and positions. Solids 
specified in that way may conveniently be represented 
by a binary tree whose leaves correspond to primitive 
shs?)es, whose internal nodes correspond to Boolean 
operations and represent subvolumes, and whose root 
represents the final solid. Often, the primitive shapes 



are expressed as the intersection of closed half-spaces. 
Commonly used half-spaces (planar, cylindrical, 
spherical) are mathematically defined as the set of 
points for which the value of a simple linear or qua- 
dratic fimction is negative or null. For practical im- 
plementation reasons, solid models are often restricted 
to be r-sets (a subclass of closed, bounded three-di- 
mensional sets with no dangling boundary elements 
and wilfa a finite number of feces and edges) [ 16]. They 
are often represented in terms of their boundary, i.e., 
a list of their feces (in turn defined in term of their 
bounding edges) often structured in an adjacency 
graph [II]. To guarantee that results of Boolean op- 
erations are r-sets, a regularized version of these op- 
erations is used. It performs the standard operations 
and then removes the dangling and interior faces and 
edges and makes sure that a valid boundary is part of 
the pointset Theoretically, this cleaning operation 
amounts to taking the topological interior ofthe point- 
sets produced by the Boolean operation (this eliminates 
the exterior dangling faces, edges, and vertices), and 
then the topological closure of the result (which 
amounts to putting a tight boundary around the point- 
set). These transformations are illustrated in Figure 1. 
In practice, the faces and edges ofthe model are clas- 
sified using neighborhoods [ 1 0] .Only the elements that 
play the appropriate role in the boundary ofthe solid 
are kept. Throughout this paper, all Boolean operations 
are regularized, unless explicitly specified otherwise. 

The regularized Boolean operations will be denoted 
U for the union, n for the intersection, - for the dif 
ference, and © for the symmetric difference. Futher- 
more, the regularized complement of any solid X will 
be denoted .J. For simplicity of notation, it is assumed 
that the Boolean operators in Boolean expression are 
ranked by decreasing priority as follows: complement, 
intersection, difference, symmetric difference, and fi- 
nally union. 

Throughout this paper it is assumed that the part 
models, or solids, are created by a sequence of oper- 
ations that add or subtract material or move and com- 
bine subsolids through Boolean operations. Therefore, 
a CSG representation of such a model always exists, 
even though the explicit construction and use of the 
CSG tree may be avoided in certain cases. 

2.3. Volume and surface features 

An important distinction pointed out in[ 13] sepa- 
rates surfacefeatures, which are collections of feces of 
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Fig. 1. Regularization: Three blocks. A, B, and C, shown in (a) are combined through a nonregularized 
Boolean expression {A n B)U {C- B)to produce the pointset shown in (b), which has a dangling face/, 
and a missing lace/j. A regularized version of the pointset can be obtained taking its interior, which is 
an open set depicted in (c)that does not contain any of its feces, and then taking the closure ofthe result 
and thus adding to the model all its feces, shown in (d). 
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Fig. 2. Suriaoe feature: Addm tow bosses (volume features) to the part (left) aeales a slot (geometric 
feature) that can be leoogpizEdly the user and interactively assodated with an intentional surgce feature 
fer later references. 



a part model[ 17-19], such as the walls and the floor 
of a slot, from volumefeatures, wiiich are full-dimen- 
sional pointsets ofthe part or of its complement, such 
as bosses and holes [ 2 0, 2 1 ] .Very rarely are both types 
simultaneously supported in the same modelling sys- 
tem (as they are in the prototype system MA- 
MOUR [4 ] ) . Pratt [13] provides a detailed discussion 
of the historical motivations, current merits, and 
drawbacks ofboth types of features and concludes that 
all surface features should be converted to volume fea- 
tures, which, although slightly more complex to sup- 
port, offer greater flexibility for interactive editing and 
more information for driving analysis and application 
programs. The author believes that both volume and 
surface features are useful for editing and that a volume 
representation of a surface feature need not always be 
derived. 

2.3.1. Volume features. Design is often done in an 
incremental manner, by first laying out the overall 
shape, and then adding or modifying details by creating 
or editing features. The creation of a geometric feature 
is necessarily accompanied by a modification of the 
volume occupied by the part and io practice always 
corresponds to either an addition or a subtraction of 
material. This transformation may be expressedas the 
union or difference between the part and the volume 
feature. Because the volume feature may be viewed as 
a sophisticated parameterized primitive shape, this ap- 
proach is particularly effective in dual modellers which 
derive a boundary representation from a CSG tree. 

The computational expense of explicitly deriving the 
effect of a Boolean operation [ 1 0]has discouraged cer- 
tain developers of solid modellers firom evaluating the 
boundary ofthe part obtained by subtracting or adding 
a volume feature. Instead, implicit features (also called 
imevaluated) have been recommended [19]. 

A boundary representation ofthe feature is directly 
derived from the designer's specification without 
checking if this representation is geometrically correct. 



For example, an implicit feature of type slot may have 
been defmed by mistake as hanging in air, away from 
the part, or buried inside the part and not accessible 
from any side. This incompatibility problem does not 
occur when intentional features are used instead of 
implicit features because intentional features, although 
unevaluated, carry no assumption as to their corre- 
sponding geometry embodiment. 

2.3.2. Surface features. The volume features re- 
sulting from shape modifying operations do not always 
provide a sufficient set of abstraction tools for inter- 
acting with the model. For example, a slot feature of 
interest for manufacturing applications may have been 
created as a side effect of adding two parallel boss fea- 
tures (Fig. 2). The slot may provide a convenient ab- 
straction for expressing engineering changes (which, 
for example, modify its width) and thus should be made 
accessible to the designer through an intentional fea- 
ture. The use of such a posteriori identified features 
requires the association of intentional features with a 
subset of an existing geometry. Often such association 
is done by interactively selecting a collection of faces 
of the part model and treating it as a surface feature. 

Information provided by surface features may be 
sufficient for some applications, such as planning for 
surface finishing operations or as specifying and ana- 
lyzing dimensions and tolerances [2]. Other applica- 
tions, such as assembly or manu&cturing planning, 
heavily rely on the manipulation of voliune fea- 
tures [ 22]. Except for simple cases, the derivation of a 
volume feature that corresponds to a surface feature 
remains an open issue [13], because there is no unique 
mapping from surface features to volume features. 
Typically a selected set of closing faces is added to a 
surface feature in order to produce a valid two-dimen- 
sional shell that unambiguously defines a volume (Fig. 
3). Desirable, or even correct, closing feces may not 
always be obtained by extending existing adjacent feces 
(Fig. 4 ) . Methods or heuristics for automatically con- 




Fig. 3. The volume of a surfece feature: By adding a ckxsing face (center)to the laces of a surface feature 
(left), one obtains a valid bcwndaty of a oorrespondmg volume feature (right). 
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Fig. 4. Complex closing &ces: No simple set of closing &ces 
exists for the pocket surface feature. 



structing such faces are currently limited to simple sit- 
uations, and the designer' antervention may sometimes 
be required to generate acceptable solutions. 

2.4, Compovndfeatures 

Several volume or surface features may overlap. For 
example, two orthogonal slots (volume features) may 
hawe a common intersection volume (Fig. 5). 

Similarly, a boss and the adjacent slot may share a 
common vertical wall (Fig. 6). Furthermore, interior 
features, such as a boss on the floor of a slot (Fig. 7), 
may be used as modifiers of other features. It is not 
alw^ necessary to capture such feature interactions 
explicitly in the data structure. For example, a geo- 
metric element (face or volume) may be shared by 
several intentional features that are independently used 
by different applications. On the other hand, a hier- 
archical organization of intentional features may be 
usefiil to represent explicitly compound features (Fig. 
8) and also patterns of features (Fig. 9) when such 
situations reflect the designer' sintentions or are im- 
portant for applications such as process planning. The 
nature ofthe geometric and topological interaction be- 
tweai the individual features of a compound feature 
should be derived, when needed, fix)m the actual ge- 
ometiy of the faces referenced by the individual fea- 
tures. 

3. Pitfalls 

Because they provide an intuitive, domain depen- 
dent, high-level vocabulary, both volume and surface 
features are good candidates for facilitating the speci- 



fication of shape modifying operations and the expres- 
sion of validity Conditions, provided that one can make 
the specification convenient and unambiguous. 

For the designer' sconvenience, these specifications 
have to be imambiguous, so that the effect of shape 
editing commands can be clearly understood and easily 
predictable, and the validity rules must precisely char- 
acterize invalid situations independently of the veri- 
fication procedures employed. They also must be con- 
venient, so that the specifications correspond to pow- 
erfijl high-level operations that produce the desired 
effect and so that validity rules are simple to formulate 
and powerfiil enough to trap common design errors. 
Furthermore, procedures for executing the shape 
modifying comm ands and for evaluating validity rules 
must be available. 

This paper shows how extensions of several tech- 
niques may be integrated to improve the specification 
and the execution of unambiguous shape modifications 
using compound or isolated volume or surface features. 
It also shows how a rich geometric representation 
scheme can be used to simplify the expression and 
evaluation of validity rules. Most of these techniques 
have been made possible by recent developments in 
geometric modelling, which must now be integrated. 
These developments will be briefly summarized, and 
their potential applications to feature-based editing of 
solid models will be demonstrated. 

3.1. Limitations of simple shape modifying techniques 
To stress the need for the approaches such as those 
proposed in this paper, this section discusses the lim- 
itations of several simple schemes that come to mind 
as possible ways of using features to modify and test 
solid models. 

3.1.1. Implicit features. Formerly mentioned im- 
plicit features may be trivially edited by modifyingtheir 
parameters. For example, an implicit slot can be moved 
and enlarged by changing its position and its width. 
However, as pointed out earlier, the existence of an 
implicit feature with specified dimensions and position 
does not guarantee that the corresponding geometric 
feature with the expected characteristics is to be found 
on the part. Thus, to produce a reliable description of 
a part, implicit features should be treated as intentional 
features, and the corresponding geometric features 
should be constructed (if possible) and their validity 
(i.e., compliance with functional requirements) as- 
sessed. 
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6. Shared &oe: Adding a Inss (volume feature) to fte modd (left) and then subtracting a slot feature 
(center) cacales a modd in M\liich the boundaries oftwo adjacent featuresoverlap along a pation rf a &e 
(right). 



3.1.2. Procedural models. To further automate the 
generation of models and allow the designer to con- 
centrate on high-level design decisions, it is suitable to 
support the automatic derivation of CAD models from 
a set of functional constraints specified by the designer. 
The flinctional requirements specifying the geometric 
characteristics of intentional features could be consid- 
ered as constraints and combined with the geometric 
constraints describing the part. A constraint solving 
system would converge to a valid solution, if such a 
solution exists, or declare that the specification (i.e., 
set of constraints) is invalid. Such a scheme would 
have the considerable advantage of supporting incom- 
plete specifications of features. For example, an inten- 
tional featare of type slot could be defined and its di- 
mensions specified, but its position would not be pro- 
vided by the designer, except for one constraint: The 
slot should be abutting on a given face of the object 

Such an automated approach has been given a se- 
rious consideration [23,24], but it carries a high com- 
putational cost and requires that designers provide a 
complete set of consistent constraints before the CAD 
system can create a model and return some use&l 
feedback. Deriving and maintaining such systems of 
constraints, especially when additive (bosses) and sub- 
tractive (pockets, slots, holes) features interfere is a 
considerable endeavor, as pointed out in [25]. Fur- 
thermore, a simple indication that the set of constraints 
is incompatible does not provide useful hints as to what 
part of the specification should be modified to produce 
the correct result. A more practical approach is to build 
models incrementally by transforming and combining 
simpler models. A procedural rather than declarative 
approach, in which the designers specify a sequence of 



operations that transform a model in attempt to satisfy 
constraints, has been described by Rossignac in [26]. 
It relies on the designer'sability to decompose the 
problem into an ordered set of subproblems that can 
be solved one at a time. The procedural specification 
(i.e., the sequence of operations that solve the individ- 
ual problems) is saved and can be edited by the designer 
and reexecuted on demand. This technique could be 
used for feature-based editing by considering the im- 
plicit features as intentional features to be created in 
apredefined order. Reexecutingthe specificationwould 
attempt to create the geometric counterpart of the in- 
tentional features at specified positions and would re- 
port whether the creation was successful (i.e., whether 
valid features have been produced). 

Consideringthat a procedural model can be obtained 
simply by storing Hie designer'scommands into a log 
file and making the file available for modification and 
leexecution fails to address three important problems: 

1. Features successfully created during an execution 
of the procedural model can be invalidated by the 
subsequent creation ofother features. Therefore, to 
assess the validity of a design, intentional features 
created at an early stage of the specification must 
be preserved and methods for accessing the corre- 
sponding geometric elements (whenever they exist) 
and for testing the compliance of these elements 
with feature validity rules should be available. 

2. Feature parameters may be defined in terms of other 
features. For example, a boss may be intended to 
lie at the center of the floor of a rectangular pocket 
Although this relation may have been established 
at the creation of the boss, subsequent editing of 






Fig. 7. Nested features: Adding a boss featare in the middle of a slot feature (left) creates a model (center) 
with a nested feature (right): The boss is inside the volume of the slot 




Fig. 8. Compound features: A hiaardiical organization of features may be used to represent a compound 
feature made of a slot that contains two boses, one of wWch has a hde (left). "Ikse relations can be 
captured by creating additional intentional features of lype compound ftat refer to ote features instead of 
refening to geometric elements (right). 



the procedural model may alter the position of the 
pocket. To preserve the relation between the posi- 
tions of the two features, this relation must be cap- 
tured in the procedural model and used to position 
the boss at each execution. 
3. Surface features are typically selected by the user 
in an interactive mode, preferably using a graphic 
cursor to pick the appropriate geometric elements 
fium a picture of the model on a computer screen. 
Surfece features in a model may be used to attach 
tolerances or other manufacturing attributes to 
particular portions of the model and to serve as 
clues for process planning [ 5 ] , or simply to provide 
constraints on the position and size of volume fea- 
tures to be created later by the procedural model. 
It is therefore essential that surface features, once 
selected by the designer on one instance of a model, 
be selected automatically when a new instance of 
the model-is created. To support automatic rese- 
lection, operations that select surface features must 
be captured in the procedural model in such a 
manner that their execution produces the desired 
result, even when earlier parts of the procedural 
model are modified. 



Techniques for supporting procedural models to- 
gether with intentional surface features and for cap- 
turing relations between such features have been de- 
veloped and implemented by Borrel, Nackman, and 
(he author, and are described in [4]. They will be briefly 
reviewed in Section 4 and their applications to feature- 
based editing will be discussed. 

3.1.3. local boundary modifications. If no proce- 
dural model is available for editing and reexecuting, 
or if the execution of the procedural model is costly, 
the part model may sometimes have to be directly ed- 
ited, or "patched." Since a valid solid model is un- 
ambiguously described by its boundary, boundary 
"tweaking" seems an attractive technique for editing. 
For example, the four vertices of the floor of the slot 
in Fig. 10 could be raised to change the depth of the 
slot. A new boundary representation would be readily 
available if the floor and the adjacent faces were im- 
plicitly defmed by their vertices. However, such 
boundary tweaking techniques may require major al- 
terations of the boundary structure of the object. If 
polyhedral modellers, if the geometry ofthe edges and 
faces ofthe model are implicitly represented in terms 
of vertex coordinates, a &oe (for example, the floor of 




Fig. 9. Patterns of features: A pattern of hole features (left) may be repesoited by a compound feature 
referencing the intentiraial featoes of each hole (center). The compound feature has a descrqjtKSi of the 
pattern paramders (right) and can be used to access and intemogate ttie entire pattern. 
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F%. 10. Boundary tweaking: The floor of a slot (left)ni£[y be 
raised by moving up its volioes (right). 



a pocket) may be moved simply by displacing its ver- 
tices. However, the vertices of each face must remain 
coplanar. Coplanarity is guaranteed by forcingthe ver- 
tices of the face being moved to remain coplanar and 
to move along extensions of the edges they bomid. 

The situation is more complicated w4ien vertices 
correspond to intersections of curved surfaces. Suppose 
that the designer wishes to move only one face of the 
feature and wants to express this modification in terms 
of vertex displacement. Then vertices are constrained 
to move along extensions of abutting (curved) edges 
vrfiile simultaneously remaining on the edited surface 
as it evolves. The number of degrees of freedom, which 
specifies how many of the vertices may be moved in- 
dependently, is dictated by the nature of that surface 
and its acceptable deformations (scaling, rotations ...). 

In addition, the validity domain, which specifiesthe 
maximum amount of vertex displacement along each 
edge so that the boundary structure (adjacency graph 
between faces) remains constant, is defined by the ge- 
ometries and relative positions of the faces. Restricting 
the displacement of each vertex along the extension of 
an edge so that it does not exceed the intersection of 
this extension with other faces is not sufBcientto guar- 
antee that the resulting boundary will not be self-in- 
tersecting (as is the case in Fig. 11). 

An alternate approach is to modify the boundary 
structure or to alter several faces simultaneously. A 
simple example may be found in Fig. 1 1, where a ver- 
tical left wall is added to connect the lowered floor to 
the abutting cylindrical fece, but in general it is difBcult 
to devise procedures for specifyingthese modifications. 
Furthermore, detecting self-intersectingboundaries is 
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computationally difficult, and no unambiguous and 
useful semantics has been defined for specifying bow 
self-intersecting boundary representations should be 
corrected. 

In conclusion, the semantics of boundary tweaking 
operations is not well defined and the results may be 
invalid. This paper thus proposes techniques that use 
mathematically well-defined set theoretic operations 
to alter features. 

A possible approach would be to convert incorrect 
surface features into volume features (by the insertion 
of closing faces) and then delete these volume features 
and create new correct ones. As pointed out earlier, 
the set of closing faces is not unique for any surface 
feature, and even a single suitable set may be hard to 
produce. To overcome the "closing face problem,'^ 
technique based on corrective volumes, previously em- 
ployed by Requicha and the author for local blending 
operations[12], will be proposed in Section 6. Here, 
it uses a generalization of sweeps or extrusions to create 
solids that can be added tO'-or subtracted fixxn— the 
part in order to change the dimensions or positions of 
geometric features, 

3.1.4. Order dependency of volumetric alterations 
A volume feature, whether directly created by a Boo- 
lean operation or derived by closing a surface feature, 
could in principle be modified by deleting it and, if 
necessary, by creating a new volume feature with the 
desired characteristics (Fig. 12). 

Deletion of a volume feature may be obtained by 
using the inverse of the Boolean operation that could 
have been used to create it. For example, a subtractive 
feature of type slot could have been created by sub- 
tracting the corresponding slot volume from some pre- 
vious representation of the part. Therefore, adding tiie 
volume back, using a Boolean union, would delete (he 
feature. Unfortunately, this approach suffers from three 
major problems, which could be qualified as "unde- 
sirable side effects." 

1 . As demonstrated in Fig. 13, the lack of associativity 
properties of set theoretic Boolean operations do 
not in general permit to undo the effect of sub- 
tracting (respectively adding) a feature by adding 
(respectively subtracting) it back. Specifically: 
{A -^ B)^Bi'A, unless B C A, and similarly 
{A\J B) - B * A,\m\QSS B C A. 



Fig. 11. Validity problems: The depth ofttie slot (left) is altered ty lowering itsfloor face. In the resulting 
model, the boundary is disoMmected and self-inlersecting.To connect Ihe boundary m the right side of the 
slot, the right waU can naturally be extended to foDow die valices. On the left side, however, the wall is 
cylindrical and the kjweted vertioes do not lie on its extensicHi. They either have to be moved horizontally, 
or a new vertical face must be aealed. In any case, the boundary is self-intersectingon the right side of the 
slot, and therefore poilifflis of it must be eliminated, which invokes gecmetric calcuMcms. 
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2. When additive and subtractive features interfere, 
the order in which they have been created is im- 
portant and in general cannot be altered without 
producing undesirable side efifeds, such as the one 
illustrated in Fig. 14. 

3. Modifying features by deleting them and creating 
new ones could modify or destroy other features 
(see Fig. 15). 

For solids created by combining simpler solids and 
volume features through Boolean operations, correct 
results may be obtained, without reevaluatingthe whole 
Boolean expression, by confining the deletion of the 
dd feature and the creation of the new correct one to 
ihe active zone of the old feature. Active zones have 
been introduced by Voelcker and the author in [27], 
and would correspond here to the portion of space 
vUxse the shape of the feature is important. A formal 
definition and relevant properties of active zones will 
be summarized in Section 5, and new properties for 
applications to feature-based editing will be derived. 

Instead of computing the active portions of all fea- 
tures (i.e., their intersection with their active zones), 
^Hoe decomposition techniques [28-30] may be used 
to precompute and store the geometry of the active 
portionsofall features simultaneously, and thus to im- 
pnave the performance of algorithms that execute the 
editing operations. Then, each editing operation may 
be confined to the appropriate regions without re- 
peating expensive geometric calculations each time. 
Such an approach has been proposed by Pratt [ 1 3 Jusing 
an extension of the radial edge structure developed by 
Wdlof 28]. Section S briefly presents a different, more 
general, and slightly simpler structure called Selective 
Geanetric Complex (abbreviated SGC) described by 
O'Connoandthe author in [29]. Algorithms for sub- 



dividing SGCs may be used to decompose space into 
cells of dimension three or less, such that given any 
cell C and any volume (or even surface) feature F, C 
either lies entirely in the interior of F, entirely in its 
complement, or entirely in the boundary of F (when 
C is a fece and F is a volume feature). Each cell of an 
SGC "remembers" what features it belongs to and is 
associated with an attribute which defines whether the 
cell is part of the represented solid or not. Fig. 16 il- 
lustrates such a decomposition. The applications of 
the SGC structure to the deletion and modification of 
features within their active zones will be demonstrated. 



3.2. Limitations of simple interrogation techniques 

Detecting whedier llie geometry associated with an 
intentional feature satisfies the validify requirements 
explicitly associated with that particular feature-fype 
is extremely convenient for v^idafing a design and 
automatically verifying that changes produced by add- 
ing new features or modifying old ones did not create 
undesirable side effects. However, the validity of each 
individual feature may not be sufficient to assess the 
validify of a complex part, and sometimes, a relation 
between several features is also important. Most of 
these validity requirements may be explicitly attached 
as rules to single or compound features[4], and au- 
tomatic procedures for checking feature validify may 
be invoked, as described in Section 7. 

Some validity rules may be characterized in terms 
of topological relations between volume features and 
can be expressed in terms of Boolean operations be- 
tween a volume feature and the part it is supposed to 
modify. For instance, when a slot B is to be subtracted 
from a part model A (Fig. 17 center), one can detect 
situations vtliere the slot is too deep and where its floor 




Fig. 13. Side eflfeds: In the slanted top ice of fte model (left), a slot b oeated (center)by subtracting a 
bkxk. Adding block bade (right)does not restore the CHiginal model 
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...y: Nfodifying a soM (left)lw first creating a slot (top center) and then a boss (top 
right) produces a (fiffitent result than fist making Ihe boss (bottom center) and then the slot (bottom right). 



is missing in the boundary B (Fig. 17 right) by 
checking whether B - A represents the empty set or 
not (this works only if the floor is not coincident with 
the bottom face of A). 

When regularized Boolean operations are used, such 
simple tests are insufficient to distinguish between sig- 
nificantly different situations. For example, a distinc- 
tion betweenthe correct situation and the unacceptable 
configuration in Fig. 18 may not be obtained as pre 
viously suggested by considering B - A, because B 
- Ais empty in both cases. Clearly, in most cases, a 
more complex test involving Boolean combinations of 
auxiliary solids may be provided, but each situation 
and each test may require different types of auxiliary 
solids. For example, a "roof," thin blockB'over B 
may be used: B'HA is empty' in the correct situation 
of Fig. 19, but is not empty in the invalid case of the 
same figure. These tests are expensive to perform and 
may require constructing and intersecting complicated 
auxiliary volumes. Furthermore, this approach is lim- 
ited to volume features. An alternate approach is pro- 
posed in Section 7 which demonstrates that most com- 

t Testing whettier a Boolean combination of two or more 
solids is eii^ requires evaluating its boundaiy or peifomiing 
a Null Object Detection test on a CSG representation[31. 



mon validity criteria may be tested by querying the 
existence of lower-dimensional parts in an SGC rep- 
resentation of the space decomposition defined in terms 
of features. 

4. PROCEDURAL MODELS 
As pointed out in the introduction, procedural 
models are particularly convenient for trial-and-error 
geometric design, because they capture the designer's 
specifications in a form that can be easily understood, 
edited, and repeatedly executed. This section briefly 
presents a prototype system for procedural modelling, 
called MAMOUR, which has been implemented in 
the object-oriented language AML/X[ 32] and de- 
scribed in more details in [4, 5J.MAM0UR supports 
intentional features with validity rules and can keep 
trade of relations between different features. 

4.1. Sequences 

Such a system offers an adequate platform for in- 
tegrating the techniques presented in this paper. As the 
user of MAMOUR creates and transforms a part 
model, the system automatically constructs a proce- 
dural model composed of an ordered set of operations 
that measure geometric (orother)properties, create or 
move primitive or subsolids, define intentional features, 
or perform Boolean operations between objects. The 




J, . , Jn amodd with a throiffih-hole (a), a vertical slot (b) is cnsated 

_ a vwoi^ positicML Deleting the slot ^ a unim creates a solid (c) in which the through-hole feature is 
partly desSoyed. Mddng a new slot in me right poation (d)produces a valid slot, but leaves the through- 
hole invalid 
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Fig. 16. Space decomposition: In the model (a) a slot is made (b). Then two bosses are created inside the 
slot (c). Space is decomposed into 3D cells (d), so that each cell is either inside or outside the original 
model and any one of the ftatuies. 



execution of the procedural model constructs a geo- 
metric model and the associated set of intentional fea- 
tures. The procedural model can be edited by the user 
withoitt interrupting the design session, because it is 
stored in memory as an aggregate (list) of AML/X 
objects.fBach object corresponds to an operation, or 
design step. Fig. 20 top illustrates one such sequence. 
Each operation has its own internal variables and ex- 
ecution methods. A whole sequence of operations may 
be captured in another AML/X object of type SE- 
QUENCE and can be edited, executed, and included 
as a parameterized macro operation into other se- 
quences. 

4.2. Vnevaluated parameter expressions 

To provide a greater flexibili^ and multiple appli- 
cations of the same sequence, each operation and 

* An object is an entity that has a type (for example, 'Drill 
operation"),a set of internal variables (for example, expres- 
sions that define the radius and position of the hole), and a 

■ocedures, called methods C ' ' " 

sa CSG representation of i 
inderofthe ^propriate radius). 



therefore each sequence is parameterized. When several 
sequences are pieced together in different ways or when 
an early part of a sequence is edited, an individual 
operation may be executed on a model different from 
the model that has been used to specify that operation. 
For example, a make-rib operation may have been used 
to make a rib in the center of the boss in the original 
model. The parameters locating the rib were derived 
from the position of the boss. If operations that con- 
struct or modify the boss have been edited so that the 
shape of the boss is changed, the execution of the make- 
rib operation will produce a rib that is no longer in the 
center of a boss, utdess make-rib can adapt its execution 
to the new geometry. This adaptability is possible in 
MAMOUR, because parameters of operations may be 
stored in an unevaluated form, thus capturing the de- 
signer's intents (for example, to position the rib in the 
center ofthe boss). To simplify the formulation ofthis 
constraint, an intentional feature is associated with the 
boss. Intentional features in MAMOUR are also AML/ 
X objects with internal variables and methods. The 
internal variables typically contain sets of unevaluated 
references to boundary elements ofthe model (for ex- 




Fig. 17. Validitytest: The position ofavolumefeatureBin relation to the partA is shown in(a). Subtracting 
B fiom A Ms to produce the desired slot feature (b) because the depth of B is too large. Sometimes, such 
situations may be detected by testing if the difference .8 - ^ is an empty solid or not (c). 



Jaroslaw R. Rossignac 




/ 




Fig. 18. Inadequate validity test: Subtracting fiom a part A (left) a correctly positioned slot volume feature 
B produces a valid geometric feature of ^pe slot (center). Subtracting an Ul-positioned slot B oeates an 
invalid geometric feature (right). The diffitence may not be detected by considering ? """"""""^ n-^ni-x*" 
" ^'inofAandB. 



ample, references to the feces ofthe boss). The methods 
can be used to evaluate these references and obtain a 
geometric description ofthe appropriate elements (for 
instance, ofthe floor face ofthe boss). These geometric 
elements are also AML/X objects with methods for 
computing their geometric properties (for example, 
the center of a face). Thus, if boss 1 is the name of 
the intentional feature that corresponds to the geo- 
metric feature of type boss, the position of the rib 
may have been specified in make-rib using 
bossl.floorO.centerO, which defines the correct posi- 
tion, as long as the geometry referenced in bossl is a 
valid boss, or at least has a floor face. 

Until now, M AMOUR has been interfaced wath only 
a 2D geometric modeller, and therefore intentional 
features contain only edge-references. Such an edge- 
reference is an unevaluated AML/X expression defined 
in terms of: 

• the structure of a CSG representation ofthe partic- 
ular model (if available), 

• the adjacency information typically available in a 
boundary representation, and 

• thenamesof operations that created or modified the 

Tools for automating the construction of expression 
that consistently identify faces in various versions of 
a model are currently under ivestigation by Paul Bor- 
rel and by the author. 

4.3. Validity tests 

The validity of intentional features, and thus ofthe 
produced model, may be tested automatically using 
validity rules attached to features. To a particular fea- 
ture may be associated several rules (predicates), which 
measure geometric and topological properties of 
boundary elements. for example, any rule evaluates 
to FALSE, the feature is considered invalid. Note that 



the same feature may be valid for another application 
with different validity criteria. A rule can be any 
expression* in AML/X that returns a Boolean value. 
Typically such expressions involve references to inten- 
tional features, and therefore indirectly to correspond- 
ing geometric entities, and also calls to methods applied 
to these entities for computing their properties or de- 
riving specific measures. 
4.4. Editingfeatwes 

AAA. Volume features. Intentional features can be 
created in MAMOUR by executing a shape-modifying 
operation, which attempts to produce a geometric fea- 
ture through a Boolean operation and at the same time 
creates an object that represents the intentional feature 
and contains, in its internal variable, references to (he 
corresponding geometric entities. Note that these en- 
tities need not always exist in the boundary represen- 
tation ofthe model as it evolves. Since a volume rep 
resentation of the geometric feature is used to the 
Boolean operation, these features correspond to the 
volumefeatures discussed above. Such volume features 
may easily be used for editing the. model. To each 
boundary element and to each intentional feature in 
MAMOUR is associated a history attribute, which in- 
dicates what operation created the entity. After a fece 
is graphically selected, its history may be accessed by' 
the designer and the parameters* ofthe corresponding 



* In particular, it can combine through an OR operation 
seva:a] Boolean subexpressions, and thus provides a simple 
mechanism for expressingvalidity rules as conjunctive fernis 
or even more complex Boolean combinations. 

» It is the designer's responsibility to decide wWch param- 
. eters should be edited and how. The systems could povide 
: some help by maintaining a d^aidoicy griph that idates 
fioes to parameter e!p]BSSion.0mn1heseremians are sin^e, 
and vw have not investigafedsucfa fidlitiesfor fiirtherassisting 
the designer. 




Fig. 19. Auxiliary volumes: A thin blodc B' over the slot B may be used to test whether the slot B is made 
in the block A (left) is accessible from tiie top. hi the valid wmfiguration (center) tiie set S' n ^ is empty, 
hi the incorrect configuration (ri^t) B'r\Ais nC " 
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Fig. 20. Sequence of operations: The sequence (top) contains two slot-making operations, w4iich subtract 
two volume features, followed by a feature-selectionoperation which identifiesa suifece feature of ^ boss 
used in the last operation to position a rib volume feature. The execution of the sequence produces a model 
(center). Editing one of the slot-making operations and reexecuting the sequence applying it to a different 
starting object produces a different result (bottom), which still reflects the designer'sintent to center the rib 
on the boss. 



operation edited so that the correct feature is created 
by reexecuting the whole sequence. This reevaluation 
vns>/ be computationally expensive. Therefore, alter- 
nate techniques for editing the model are discussed 
below. 

4.4.2. Surface features. Intentional features may 
also be created without geometric modifications by a 
select operation which takes as parameters a feature 
type and a list of references to existing boundary ele- 
ments. §Typically there is no single she^ modifying 
operation responsiblefor creating such a sur&ce feature . 
and thus no explicit volume representation. Techniques 
for editing such surface features by adding another op- 
eration to the sequence are addressed in the next sec- 
tion. They complement facilities for editing operations 
already in the sequence, which are well suited for mod- 
ifying volume features when the cost of reevaluating 
the entire sequence is not prohibitive. 

5. EFTICIENT EDITING OF VOLUME FEATURES 
Active zones were introduced by Voelcker and the 
author in [27] to speed up certain geometric compu- 
tations over Constructive Solid Geometry (CSG)rep- 
resentations. The active zone associated with a CSG 
iMde is the region in which the shape of the set rep- 
resented by the node is important. An active zone is 
defined algebraically as the intersection of certain nodes 
of the whole CSG tree, and thus its CSG expression is 
ahvays available. The concept of active zones is applied 
in this paper to deal with editing feature. In this section. 



* Automatic extraction of existing 3D features is a difScult 
and expensive process[33] , especially when loose conditions 
are used to identify features. For example, a slot with a filleted 
bottom edge may still be a slot, although the floor fece is not 
direcdy connected to the wall. In MAMOUR, references to 
edges of features may be specified in a semiautomatic way 
using graphic interaction. These references are then integrated 
into unevaluated expressions, stored with the model, that will 
identify the corresponding feature in a &mily of part models 
genersied by reexecuting a modified version of the sequence. 



fixm the properties of active zones, a CSG expression 
is derived, which, in conjunction with spatial localiza- 
tion techniques, may be used to improve the perfor- 
mance of algorithms that update the boundary rep- 
resentation of a model when a volume feature is altered. 
These improvements are particularly interesting when 
a spatial decomposition scheme is used for the bound- 
ary representation. 

5.1. Active zones 

Let She a CSG representation of a solid encoded 
in a binary tree. (For simplicity, wte shall use the same 
symbol for the solid and its CSG tree.) Let A be any 
primitive (or any internal node) of S. The active zone 
ofA in S, denoted Z^, is equal to the Boolean differ- 
ence /5 - Ua, where and Ua are respectively the 
I-zone and the U-zone of A in S. Let A be a primitive 
or internal node of a subtree N m S, and let F be the 
parent node of N of another node B. The following 
properties provide a recursive formulation for incre- 
mentally constructing /i Ua, andZi. 

• li is the universe (Euclidean three-dimensional 
space) and Ui is the empty set. 

• When F = U B or /• = fl U M then y = I,", 
C/J = l/J U B, and Zj •= Zj - B, 

• When F = N O B OT F B H N, then V = X' 
n £/5 = um, and Z5 = Z^ n B. 

• When F = B - iV, then /J = B - U^, £/J = /J. 
and ZJ = ZiJnB, 

• When f = AT - B, then /5 = /2f n J, = Ul. 
andZ5=Z2'-B. 

It follows that CSG expressions for and Ua, and 
thus for Zj, may be computed by traversing the pathl 
in the tree d-ora A to S, and constructing the Boolean 
expressions for 1-zones and U-zones using intersections 



11 Thepatisdefinedasthe set of nodes traversed by moving 
fi-om S to A in tfie CSG tree, following the parent-to-child 
links. 
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Fig. 2 1. Chanses are confined to the active zone: Replacing the primitive A by X in Ihe CSG definition of 
S (left) can only affect Sin the active zone of A (center). The result is show (right). 



with other branching nodes" or with their complements. 
On the other hand, these expressions need not be ex- 
plicitly constructed, nor stored, for most applications, 
since efficient algorithms presented in [27] can perform 
calculations with respect to ^3 and t/5 using the ap- 
propriate branching nodes in the original tree S. 

Several applications of Active Zones are discussed 
in details in [27]. For instance, it is shown that if a 
node A does not intersect its active zone, it can be 
replaced by the empty set. Similarly, if a node A con- 
tains its active zone, it can be replaced by the universe 
without affecting the set represented by flie whole tree. 
(These results provide a new algorithm for testing 
whether any particular node is redundant .) Algorithms 
for the detection of interferences among solids defined 
in CSG and for the generation of shaded pictures di- 
rectly from CSG by ray-casting[34] or depth-buffer- 
ing[9] are based on facilities for classifyii^ subsets of 
the boundary of primitives against the CSG tree, The 
performance of these algorithms may be improved us- 
ing only the I-zones of the primitives, and not the whole 
CSG tree for the classification. 

5.2. Editing volumefeatures 

5.2.1. Localization to the active zone. Because a 
volume feature corresponds to an internal node ofthe 
CSG tree, it has an active zone, an I-zone, and a U- 
zone. It was established in [27] that all changes to a 
primitive or node A inside its active zone Z5 in Swill 
affect the shape of S and conversely that changes to A 
outside of Z5 will leave S unchanged. This property 
and related properties are used in this section to pro- 
duce trimming expressions for localizing changes that 
implement a volume feature modification. For ex- 
ample, let A be a volume feature subtracted in the se- 
quence (B - A) - C defining a solid S (see Fig. 21 ); 
its active zone Z5 is B - C. Only changes to A inside 
B - C need to be taken into account. 

5.2.2. Localization to the altered portion of the fea- 
ture. When a volume feature A is edited and replaced 
by X in a CSG representation of the solid S, a new 
boundary may be obtained without reevaluating the 



" The sjintol B in the above properties refes to vM is 
fixmally called a brandling node of A in S, if., a node that 
does not lie in the palfa ficm A to S, but v\iiose parent node 
does. J A and are defined as intersecdonsof such branching 
nodes or of their complements. 



entire boundary of the solid because changes are re- 
stricted to A © X,the symmetric difference* between 
A and X [35], Conseqijently, in the example Fig. 21, 
only alterations in " A') i i (5 - C) are needed in 
order to compute the boundary of the solid obtained 
by replacing A by X in the CSG tree of S. 

5.2.3. Combined localization. Let o and w represent 
respectively the empty set and the universe. Let Sx 
denote the set represented by the tree of S in which A 
was replaced by X. The two previous results may be 
combined to yield: Sx'-^.S®.Zin{X®A). A proof 
of this equation may be derived from Equation (4) 
in [36] , given that(A'©y=Z)*»(^*y©Z) holds 
for any three sets,X, Y, and Z. 

5.2.4. Cbssifying additive and subtractive parts. The 
symmetric difference A © Xmay be decomposed into 
two disjoint sets,JA - X and X - A, wdiich can be 
treated separately. For example, when a negative (sub- 
tractive) volume feature A is replaced by X, a subset 
of A - X must be added to S and a subset of X - A 
mustbe subtracted (see Fig. 22 for an example). This 
section derives new CSG expressions that characterize 
these subsets and are simpler than the CSG expression 
ofs§. 

Although Z 3 is the smallest region where changes 
to A affect S, the classification of A - X and X - A 
against the CSG expression of Z5 in general contains 
unnecessary steps which can be eliminated using the 
following property; 

Sx = s- {Ax - ui) 0 {At n 

v/hetc Ax is the portion of material added to A and 
where is the portion of material subtracted fixxn 
A during the same operation. 

Fu^t consider the case of an additive feature A . We 
lme:Ai =X-AandAi = A -X. For a proof, first 
we establish the foUovting: 



♦ The symmetric difference between A and X is defined as: 
A : X-{A-X)U{X-^A). 

* The tema disjoint is used here Ihree-dimensiOTal vol- 
umes vitose inteisection is empty or of lower dimension. 

I Note that each one of these CSG ejqaessions may be fiif- 
ther trimmed down b/ removing all primitives that are da^iMA 
ficmA -X(Kq)ectivelyX-A)by usingtediniquesdisaissed 
in [37]. 
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Fig 22. Changes are confined to the symmetric difference: The solid S (left) is represented by a CSG tree 
defined in terms of a subtractive volume feature A . Replacing A by X in the CSG tree of S may change the 
shape of S (center) only in the symmetric difference A © X (right). Note that A - X is added and that 
X- A is subtracted. 



/3 n c/J c /5 n u J, n t/3. 
/3nt/3c(/3u5.)ni/3. 

/J n C/5C5»n by Proposition 5 of [27], 
/In 17J C So, by Proposition 6 of [27]. 



Now the proof may be constructed as follows: 



S;f= 5<,U(ATl/3), by Proposition 1 of [27], 

5;f = u (.f n /3) u (-4 n /5 n c/3). 

since /Jnt/Jc^^ (see above), 

5^= (S,ni)U(5„nx)us„ 

u(^n/5)U(.4n/3n c/3), 

. since ^..nyiu 5on^e So, 

Sx = (So n ^) u {So n x) 

u (i',, n {;5) u (^ n /5) u (^ n /J n t/i), 

since 5„Cf/5 (see [27]). 

S;f= (5„n^)U(S,n<^r) 

u (50 i/S) u (^ n ATI /3 n i n Arn /5) 

U(^n75ni;3). since /iUi=»i', 

S;t ■? (So n i) g (5. n AT) u (s, n 1/3) 
u(^o/3ni)U(-4njrn/3) 

u(.in/3nz/3)U(inArn/3). 

S;f-(S.u^n/5)n(.?uA'u t;3) 

U (in ATI /3). factoring outs, and iUJT. 

S;f = s n (i u u t;3) u (i n AT n /3), 

replacing X with in Proposition 1 of [ 27 ] , 

Sx-s- [.{A - A) - 1/3) u {(^-yi) 0/3), 

using complementation. 

Now note that for a subtractive feature, one may 
oraisiderthe complements of A and X to be additive 
features, and the same proof holds. 



Since the interiors of i4 J and of are disjoint, 
the interiorsof4i - 1/3 and .4 J n /3 are also disjoint. 
Therefore, since Uand ~ are regularized, we also have: 

S;r-SU(^5n/3)-(^J-t/3). 

Consequently, subsets diA x need only be classified 
against /3; subsets inside /5 should be added to S, 
while subsets outside may be discarded.' Similarly, 
subsets of /I J need only be classified against C/3; sub- 
sets outside L/3 should be subtracted fi^om S, while 
subsets inside may be discarded." 

The CSG expression of Z5, defined as /3 - 1/3, is 
a complex as the combination of the CSG expressions 
for /3 and for C/3 . and classifying* a point witii respect 
to 7. A may often require classifying it against I," and 
against C/f and then combining the result. The above 
result suggests that A \ needs only be classified against 
I A and not against both /3 andt/5. Similarly, .4 J 
needs only be classified against 1/3 and not against 
both I A and Ua. Significant improvements for algo- 
rithms that perfonn feature modification result. To 
further improve the performance of these algorithms, 
tree pruning and other space subdivision tech- 
niques[39, 40, 35, 41] may be used in conjunction 
with classification against /3 and £/3.Note, however, 
tiiat this speed-up may suggest to add to S portions of 
i4 J that are already included in S , o r to subtract from 
Sportions of/ljr that are not in S. These unnecessary 
additions do not invalidate the result but may corre- 
spond to redundant geometric calculations unless the 
representation scheme presented below is used. 

5.3. Mixed-dimensional geometric model 

To represent a solid part together with its additive 
and subtractive volume features, one needs to extend 
the boundary representation scheme and to provide 
support for representing decompositions of solids and 
of their complements in terms of boundaries of poten- 
tially overlapping features (Fig. 16). Furthermore, the 
construction of closing faces for surface features may 



1 See Properly 13 in [27]. 
"See Property 12 in [27]. 

♦ "Classifyinga point against a set consists of determining 
whether the point lies inside or outside of the set[38]. 
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Rg. 23. Space deconpjatioi for editing a sur6ce feature: The slot surfece feature in the sdid B U C U D 
(left) may he modified by creating the oaresponding volume feature A and its modified vaaoi X, by 

' — ' — -n SGC reptesenlation of the lesullingqaDe subdiviaon (center)and by setting the status of 

odb of A - X to IN, and of cells of X - A to OOT (right). 



involve splitting faces of the part model and inserting 
internal faces in the part or creating external faces in 
its complement. Objects with internal structures and 
external dangling edges and faces are not supported by 
conventional solid modellers, and thus richer schemes 
for geometric representations must be used. 

Several data structures that support representations 
of unions of quasi-disjoint subvolumes have been pro- 
posed (see [28-30, 42] for examples). Such schemes 
may be used to decompose a volume feature B into its 
active and inactive portions with respect to various 
Boolean combinations of other features. The Selective 
Geometric Complex (abbreviated SGC) data structure 
permits representation of such decompositions together 
with the modelled part, whether the feature has been 
added or subtracted [29], 

SGCs provide a common firework for represent- 
ing objects of mixed dimensionality, possibly with in- 
ternal structures and incomplete boundaries. SGCs are 
composed oflinite coUectionsofmutually disjoint cells, 
which are open coimected subsets of n-dimensional 
manifolds and generalize the concepts of edges, faces, 
and vertices used in most solid modellers. The con- 
nectivity between such cells is captured in a veiy simple 
incidence graph, vAiose links indicate "is-a-boundary- 
of relations between cells. By setting the status (at- 
tribute) of certain cells to IN and others to OUT, one 
can associate various pointsets with a single collection 
of cells. When a cell'sstatus is IN, the cell is called 
"active," and the pointset spaimed by the cell is con- 
sidered as part of the object; otherwise it is considered 
as part of the complement of the object. Consequently, 
the pointset of an SGC object needs not be homoge- 
neous in dimension, nor even be closed or boimded. 
To support useful operations on SGCs. Boolean and 
other set-theoretical operations (closure, interior, 
boundary) have been decomposed into combinations 
of three fundamental steps for which dimension-in- 
dependent algorithms have been developed [29] : 



• a subdivision step, which makes two objects "com- 
patible"by subdividing the cells of each object at 
their intersections with cells of the other object; 

• a selection step, wdiich defines "active" cells, i.e., cells 
w*ose status is IN; and 

• a simplification step, which, by deleting or merging 
certain cells, reduces the complexity of an object's 



representation wifliout changing the represented 
pointset and without destroying decompositions that 
are marked as important for applications. 
Furthermore, combinations of these steps may pro- 
duce a variety of special-purpose operations whose ef- 
fect is controlled by simple predicates, or filters, for 
cell selection. 

The subdivision step may be used to create a space 
decomposition that reflects the geometry of the part 
and of all its volume and surface features. Then, feature 
modifications may be performed by modifyingthe sta- 
tus of cells that belong to the appropriate volume fear 
tures and tiie associated lorUzones. Fig. 23 illustrates 
how the result of a sequence of operations may be 
modified by selecting a surface feature and modifying 
it by replacing the volume feature A with X. 

5.3.1. Application to model updating. A space de- 
composition technique may be used to split A © X 
into connected cells that are entirely inside or entirely 
outside of both li and t/3 . Adding these cells to Sor 
subtractingthem from Sonly requires settingtheir sta- 
tus to IN our OUT. If such an approach is used, the 
main cost lies in the insertion ofXin the space decom- 
position and in the classification of each cell of A © X. 
The classification may simply be obtained by evalu- 
ating a Boolean expression. Using I a or t/f instead of 
the expression for Zj considerably reduces the cost of 
updating the model to reflect a feature modification 
(see Fig. 24). 

5.3.2. Volumefeature deletion. A particular case of 
editing is deletion. A volume feature could be deleted 
by altering a specification stored in a procedural model 
and by reexecuting the procedural model to create a 
new geometric model without the undesirable feature 
The same result could be obtained directly by changing 
the status of all cells that lie in the intersection of that 
volume feature with its active zone. 

For example, to delete a feature A in (Bu A) - C. 
it suffices to change the status of cells in A - (BU C). 
the active portion of A (Fig. 25). These changes can 
be done by traversing 'all the cells in A and classifying 

♦ Vik assume that eadi feature has a reference to all the 
cells that it includes or that tiiese cells may be efficiently iden- 
tified in the SGC. If no provision is made for such a direct 
access, all cells of spedfic dimoisiois may have to be traversed 
to see if they belong to the apfxopriate feature. 
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Fig. 24. Active zone in feature modification: Sdefined as (B U A ) - Cis shewn (left) superimposed 
^Hce deoomposidcm obtained tyjusin^ the 



shifted to 



is a misplaced additive volume feature m Sand must be replaced with a similar feature X sligh ., 

the fcfi and towards. To obtain the correct object, first cells that belongto A - X ate clasafied against Ui 
(here t/5 = B). Hie status of cells of A - X outside of !/5 is set to OUT. Note that one of these odk was 
already OUT (center left). Then, ceDs of X - A are dassified against I,"(here /5 = 0- The status of the 
cells of X - A inside /3 is set to IN (center right). Note that one of tee cells was alreaify IN. The result 
(cells whose attributes are IN) is shown (right). 



each one against B U C. Since each cell contains (in 
its history) a list of the features to which it belongs, 
this classification amounts simply to evaluating a 
Boolean expression. 

As discussed earlier, instead of changing the status 
of only the cells of A that lie in its active zone, one 
could simply set to OUT the status of cells of A that 
lie in its U-zone when A is an additive volume feature, 
or set to IN the status of cells of A that lie in its I-zone 
vkeci A is a subtractive volume feature. 

For example, in the solid (B U A) - C (Fig. 25) , 
the U -zone, U a = B, is simpler than the active zone, 
= sue, and thus a saving in classification time 
is achieved. This simplistic example does not reflect 
ttie importance of such savings, and one should con- 
sider that classification has to be performed for a large 
number of cells and that it may involve large Boolean 
expressions for both and I^. 

6. CORRECTIVE VOLUMES 

If a CSG tree representing the model is available, 
surface features could in principle be corrected by ed- 
iting the tree. Suppose that no single node of the tree 
represents the appropriate volume feature. Each face 
of the feature may still be associated with one or more 
half-spaces of the CSG tree, so that modifying the po- 
sition or shape of these half-spaces will affect the face, 
and thus the feature. 

Editing half-spaces directly in the original CSG tree 
has the following drawbacks: 



1. It may be difficult to recognize which half-space 
should be edited. Typically the boundary of more 
than one half-space coincides with each fece, and 
techniques proposed in [27] for classifying half- 
spaces against their active zones (to predict wiiich 
half-spaceswill affect which portion of a particular 
face) may not be sufficient. 

2. Editing half-spaces in a CSG tree to alter a particular 
portion of space will often produce undesired side- 
effects in other locations. 

3. Editing half-spaces and not solid primitives can 
produce unbounded subsolids and makes it very 
difficult to implement popular performance-im- 
proving techniques for CSG algorithms that use 
bounding boxes around CSG primitives. 

4. Designers' intentions expressed in a posteriori al- 
terations of half-spaces are difficult to specify and 
to capturein the sequence of a procedural model. 
Consequently, the edited specification is not well 
suited for reparameterization or reuse in a different 
context It is not even suited for fiirlher editing that 
involves a reevaluation of the sequence. 



J the CSG tree to regroup the relevant 
half-spaces of the feature into a single node that can 
be edited as a feature is not always possible since the 
result of evaluating a general Boolean expression is 
order-dependent. Consequently, this section proposes 
touse Boolean operations to edit sur&ce features (the 
previous section dealt with volume features). These 




Fig. 25. Ddedcm using the Active Zraie: Given the set Sdefined by (B UA) - C (left)flie additive volume 
feture A may be deleted by changrag the active attribute of all tiie cells that lie in flie active pcwtioi (f A 
(center), which isA - (BU C). The result is shown (right). 
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Fig. 26. Face extrusions: Extrudinga planar base fioe (left) produces a volume (right) that is bounded by 
ttie base face and its offeet by r, and by skte feces obtamed by sweepinglhe edges of the base fece. 



operations involve corrective volumes which are added 
to— or subtracted from— a solid model to alter a par- 
ticular feature. These corrective volumes could be de- 
rived from volume features themselves derived from 
surface features by generating closing faces, but, as 
pointed out in Section 2, the derivation of closing faces 
remains a research issue. Furthermore, once a volume 
feature A is derived from the surface feature, it is not 
easy to produce a modified volume feature X so that 
A - X and X - A may be used as corrective volumes 
and added to— or subtracted from— the model to alter 
the surface feature without side-effects. Instead, the 
author proposes to construct corrective volumes by 
extruding appropriate &cesofthe surface feature. Ex- 
amples of these extrusions are provided in the Mowing 
part of this section. 

Once collective volumes are computed, whether 
through closing faces and volume features or through 
fece extrusions, they must often be trimmed to avoid 
undesirable side-effects before they can be combined 
with the solid. Automatically computing flie correct 
trimming expression is impossible unless the expression 
"undesirable side-effect'ls formally defmed. Indeed, 
the correctness of a feature editing operation depends 
on the function of the edited feature and on the func- 
tion of its geometric relation with other features. 

6.1. Extrusion of faces 

A simplest corrective volume may be obtained by 
extruding a planar face called the base face along the 
normal to its supporting plane (see Fig. 26 for an ex- 
ample). Such a corrective volume may often be ade- 
quate to change the width or depth of features that 
have orthogonally oriented planar feces. 

The extrusion V,{F) of a planar face F by a distance 
r is a volume formally defined by 

yA.F) =: {p + tn/rwithO s t S r and p 6 F}i" 



where n/r is the normal imit vector to the surface con- 
taining F. 

For faces on curved surfaces, a normal extrusion 
will be used (Fig. 27). It is simple extension of the 
above extrusion. The normal extrusion K,(F) of a 
curved base face Fby a distance r is a volume formally 
defined as: 

VAF) = { p + rnf(p) with 0 s r :s r and p e /F}, 

where nHp) is the normal unit vector to F at point p 
and where iF isthe relative interior' of F with respect 
to its supporting surface. The orientation of the normal 
b chosen in a consistent manner throughout F. Note 
that, for the semialgebraic surfeces popular in solid 
modelling, n^p) is well defined forthe smooth portions 
of F (the relative interior of the base face), but needs 
not be well defined for the bounding edges of F, its 
cusps, or singulanties.f Therefore, V,(F) is obtained 
by extruding a nonclosed face and thus does not neces- 
sarily contain all its boundary; it needs to be regular- 
ized. iif(p) can often be computed from the cross 
product of partial derivatives when F is defmed in 
parametric form, or from a gradient when F is defined 
by an implicit equation. Such extrusion volumes are 
very simple to obtain for natural surfaces (see [43]). 

Applications of extrusions and normal extrusions to 
feature modification are illustrated in Figs. 28 and 29, 
where the extrusions are used as corrective volumes 
and added to the part to modify a feature. 

Such eq)plications are clearly limited to simple cases 
where, for example, the extruded base fece and 
abutting faces meet at a right angle. Fig. 30 shows a 

TTie interiorof a fioe is the feoe minus its boundingedges, 
cusps, or singular points, 
t The normal to a cone is not defined as its qjex. 




Vr{F) 



Fig. 27. Normal extnjsion fijr curved base feces: a nonnal extrusion of a curved face that lies on a cylinder 
(left) is a volume (right) bounded by Ite wiginal fice, its offiet by r, and by side feces that are subsets of 
mied surfeoes sustained by the ed^ of the onginal face. Note that the side feces are in the closure of V,(F) 
but are not in Vr(F). 
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Fig. 28. Application of extrusion: To change the depth of the slot feature (left), a corrective volume (center) 
may be generated 1^ extruding the floor of the slot and subtracted from the part (right). 



counterexample for which a corrective volume gen- 
erated by extruding a base face is inadequate. 

To circumvent such limitations, an extended cor- 
rective volume may be generated and then trimmed 
using the abutting feces (see Fig. 31). 

The extended corrective volume may be obtained 
by extruding an extension of the base face. A good 
candidate to use as base face extension is the entire 
surfacethat contains the base face, but in certain cases, 
to avoid side-effects in distant areas, it may be pref- 
erable to consider only a simple subset of that surface. 
This technique has been proposed by Requicha and 
the author in [ 1 2] for generating local fillets and blends 
by ( 1 ) growing and shrinking Boolean combinations 
of appropriate half-spaces, by (2) subtracting the result 
fiiom the original combination to obtain an extended 
blend, by (3 ) trimming the blend to the desired shape 
near its ends, and by (4) adding the resulting corrective 
volume to the part or subtracting it. 

Extended corrective volumes derived from a single 
base fece are convenient for modifying a single di- 
mension of a simpler feature (e.g., the vridth of a slot 
orthe radius ofa hole). In general, however a corrective 
volume involves more than one base fece (for example, 
wtei the depth of a pocket with an uneven floor is to 
bechanged). For such cases, instead of combining sev- 
aal corrective volumes, a single extended corrective 
volume may be obtained by extruding several faces 
alonga common direction (Fig. 32). Such a generalized 
extrusion, Vt(p) of a set effaces F along a direction 
u, may be formally defined as the Minkowski sum [44] 
of F with the line segment joining the origin 0 with 
the point O +u. And thus: 



K(F)= {p + ?u,withpeFforO:S/£r}. 



6.2. Automatic derivation of a default trimming 
expression 

Picking the base face and providing an oflfeet distance 
does not in general provide an unambiguous specifi- 
cation for the corrective volume. It is thus necessary 
to provide facilities for automatically creating supersets 
of the desired corrective volumes and trimming 
expressions that produce the desired subsets. Trimming 
operations may involve nontrivial Boolean combina- 
tions of the half-spaces bounded by the abutting faces. 
This Boolean combination may, in principle, be de- 
rived from the entire CSG tree, if available, but this 
derivation may be difficult and will oflen require hu- 
man intervention. A possible approach to assist the 
designer and si;%gest a default trimming CSG e?qires- 
sipn is to consider only half-spaces bounded by the 
feces of the solid that share edges with the base fece 
and to eliminate other half-spaces firam the tree. The 
correct Boolean combination of these half-spaces may 
be hard to derive, and the results may still be incorrect 
(see Fig. 33). The designer may have to provide an 
auxiliary trimming expression, but the system should 
be able to suggest a good default trimming expression. 

Clearly, corrective volumes should be confined to a 
region where they do not destroy the effects of other 
features. For example, they may be trimmed by some 
other (but not necessarily all other) volume features. 

Unfortunately, a surface feature, F, typically does 
not correspond to any individual node in the CSG tree 
and therefore is not associated with an active zone; the 
results derived in Section 5 may not be directly applied 
here. 

To provide an often reasonable default trimming 
expression, the concept of a virtual active zone may be 
used. The virtual active zone ofa surface feature F is 
defined as the active zone of the lowest node Tin the 




Fig 29. Application of normal extrusion: To change the radius of a cylindrical hole (left), a corrective 
volume (center) may be generated through a normal extrusion of the cylindrical base fece and added to the 
part (right). 
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CSG tree of Sat which the surface feature F appears. 
This technique may be automated, but is only effective 
if at T one could delete and create again the corre- 
sponding volume feature without changing the final 
result. In other cases, fine tuning "by hand" the au- 
tomatically generated virtual active zone may be nec- 
essary. 

7. VALIDITY TESTS 

Previous sections dealt with techniques for assisting 
the designer in performing object modifications using 
a feature-oriented syntax. It has been assumed that no 
automated solution exists and that human intervention 
is necessary to correct the side-effect of these editing 
operations. To further assist the designer, the system 
should support facilities for interrogating important 
properties of features. We shall refer to these properties 
using the global term of validity. 

The validity of a feature or of a compound feature 
greatly depends on the nature of the feature and on 
the function played by the feature in a particular ap- 
plication. Addressed here are only the qualitative (or 
discrete) validity issues that can be expressed id terms 
of the presence or absence of specific cells in SGC 
structures. (Many other quantitative criteria may be 
easily addressed through geometric measures, derived 
from features, as described in [4].) Cells of interest 
have a specified dimension and belong to the appro- 
priate combination of features. Note that this approach 
is not based on matching subsetsoflhe adjacency graph 
of a boundary representation, but on the interogation 



of a rich data structure that captures various geometric 
entities and their relations. 

Filters for cell selection in SGCs may be used for 
validity testing. The next section introduces the query 
operators, which provide the vocabulary necessary fii 
expressing the filters. The subsequent section dem- 
onstrates their application on a &w simple examples. 

7.1. Interrogation operators for SGCs 

Filters will be expressed in terms of the followmg 
queries that can be made to an SGC, O, or to a par- 
ticular cell, C, of 0. A list of typical filters, expressed 
as methods of cell or SGC objects, follows. 

• 0.cells(k ) returns the collection of cells of dimension 
kinO. 

• O.skeleton(k) returns the collection of cells of di- 
mension less or equal to k in 0. 

• C.boundary returns the collection of cells that 
bound C. 

• Cstar returns the collection of cells bounded by C. 

• C.dimension returns the dimension of C. 

• Given a cell D of Cstar such that D.dimension 
= C.dunension + 1, C.nbhd(D) returns leftdir, 
rightdir, or bothdirs. The value lefldu- means that, 
given the definition of "left" and "right" with respect 
to an orientation of C in the manifold containing D 
as an open subset, D lies on the "left" of C. Similarly, 
when C.nbhd(D) returns rightdir, C bounds D on 
the "right." bothdirs means that C is an interior 
boundary "surrounded" by D, or more precisely that 




Fig. 3 1. Extended corrective volumes: To raise the floor of the slot (left) without modifying the position of 
its walls, an extended corrective volume is generated and trimmed (center). The result is added to the part 
(right). 
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Fig. 32. Crarective swept volumes: To raise a slot in the model (left) without changing the sh^ of its floor, 
a corrective volume is obtained by swearing upwards all the &)es of the slot (center left), and then by 
trimming the result (center right) and adding it to Ifae model (right). The dimensions of the floor remain 
unchanged and IfaereKire the swqjt volumes of the diierait &£S have different Uiidaiess. 



C is contained in the topological interior (with re- 
spect to the manifold containing D) ofthe topological 
closure of D. 

• C.histoiy returns the set of features to which C be- 
longs. 

• Finally, C.status returns IN or OUT depending 
whether the cell is considered as actively contributing 
to the pointset of 0 or not. 

For more formal definitions of these operators, the 
reader should consult [29]. 



7.2. Exami. 

Validity criteria are domain dependent, and the goal 
of this paper is not to derive them but to illustrate a 
technique forexpressingthem. Three simple examples 
involving a block of material and two volume features 
will be used. Simple tests that characterize each situ- 
ation are proposed. These tests do not involve any geo- 
metric calculations, but simply search the SGC struc- 
ture for cells that satisfy appropriate selection criteria, 
h fact, to improve the interrogation performance, ref- 
erences to cells of the SGC could be organized in a 
data base and accessed using their history, dimension, 
or other characteristics by standard data base queries. 

The local inaccessibility from the top of a slot volume 
feature A in a part B may be detected by checking 
whether the "roof face, Fl, of A is connected on the 
outside (with respect to A) to a full-dimensional cell 
ofB (Fig. 34). The test may be performed by selecting 



cells D of Fl .star that contains B in their history and 
such that Fl.nbhd(D) = leftdir. (For simplicity, w 
assume that Fl is oriented so that the lefldir direction 
points toward the outside of A) 

Of course, the foregoing filter is not adequate for 
global accessibility, which may require performing 
Boolean operations on auxiliary volumes, such as the 
volume swept by the portion ofthe feature that is visible 
from the direction fi-om which the feature is to be ac- 
cessed. 

To find whether two slotsA and Care adjacent along 
a common face or not (Fig. 35), it suffices to inquire 
whether a 2D cell, or face, F2 exists such that it bounds 
a 3D cell of A on one side and a 3D cell of C on the 
other side. 

The search may be confined to the common 2D 
cells ofthe boundaries of A and C, which may be ac- 
cessed directly if a directory of cells that belong to each 
feature is maintained. 

To test whether a slot C is contained in aslotA (Fig. 
36), or more generally whether A and C interfere, it 
suffices to query if there is a three-dimensional cell 
whose history contains both A and C as well as R 
Again, the selection may be efficiently performed by 
standard data base queries on the directories of cells 
that belong to A, B, and C organized by dimension. 

8. CONCLUSION 
Interactive editing of CAD models may be simplified 
by the use of volume and surface features. Surface fea- 




Fig. 33. Incorrect trimming CSG: To raise the floor of fee slit in the model (left), an extajded corrective 
volume is computed (center left). CSG expressions defined only in terms of hafifspaces bounded by the 
&cesof the ori^nal model are not adequate KTtrimminglfae ooiectivevcdume. An example of an inoonectly 
trimmed canective volume is i&wn (centerri^t), and the resulting model is shown (right). 
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34. Testing top aocessibilily: To distm^iishbetwem the correct pc 
the incoirect one (right), one can simph[ inquire v\4iellier, in the oa 
space deconipositicn,the roof &ce El cffslotA isconneclBdonitsoi 



tures are interactively selected the user &om existing 
faces of a model. Volume features may be created by 
addition or subtraction of material, or derived from 
surface features by defining the necessary closing faces. 
This paper addresses the issue of interrogating such 
volume and surface features as to their validity and 
their relation to each other and to the final part. It also 
proposes techniques for editing the model by modifying 
or deleting its features. 

Interrogation techniques are based on a scheme for 
representingmixed-dimensionalpointsets with internal 
structures. In that scheme, space, i.e., the modelled 
part and its complement, is subdivided into connected 
cells (volumes, faces, edges, and vertices) such that all 
points of any given cell belong to the same set of fea- 
tures. Feature validity and particular relations among 
features may be easily characterized by the existence 
or the absence of cells of specific dimensions associated 
with specific sets of features. 

Model editing to alter or delete a volume feature 
may be performed by using a procedural representation 
of the designer's specification, locating in it the com- 
mand that created a particular volume feature to be 
modified, editing this command, and reexecuting the 
entire procedural model. Reexecuting means comput- 
ing the boundary of the geometric model, which may 
be an expensive process. When a CSG expression for 
the part is available, the reexecution may be Innited 
to a particular domain defined by the intersection of 
the volume feature with its active zone. This technique 



may be adapted to surface featuresby providing closing 
&ces which define corresponding volume features and 

by considering the role these volume features play with 
respect to the CSG tree. On the other hand, surface 
features may be directly modified by constructing cor- 
rective volumes and by combining them to the part 
model through addition or subtraction. The use of ex- 
tended geometric representations of sets with intemal 
and external structures permits calculation of the effect 
of feature deletion without further geometric calcula- 
tions. 

Global access to the geometric elements of a partic- 
ular feature is provided through intentional features, 
to which may be associated validity rules and methods 
for evaluating these rules or for measuring important 
properties of tlie feature. Intentional features may be 
created automatically when feature-based shape-mod- 
ifying operations are used or by interactively selecting 
existing faces. Intentional featuresdo not directly point 
to any geometric entity, but carry an unevaluated 
e3q)ression, constructed at feature creation or selection, 
vAach, when evaluated, returns appropriate geometric 
elements, if they exist. This indirect approach prevents 
inconsistencies between an abstract list of assumed 
features that could characterize some important aspects 
of a part and the actual presence and geometry oflhese 
features in the part. This point is essential if further 
shape modifications, done either by editing and re- 
playing the designer's specification or by creating new 
features, can alter the geometry of a previously defined 




Eg 35. Testing adjacency: For pnxsss planning, the two adjacoit slots A and C (left) must be treated 
difierently ficm the two disconnected slots (right). 




feature to the point that it no longer exhibits the ex- 
pected geometric characteristics. 
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ABSTRACT 

Solid modeling is concerned with the construction and manipulation of unambiguous computer 
representations of solid objects. These representations permit to distinguish between the interior, 
the boundary, and the complement of a solid. They are conveniently specified in CSG (Construc- 
tive Solid Geometry) by a construction tree that has solid primitives as leaves and ri^d body 
motions or regularized Boolean operations as internal nodes. Algorithms for classifying sets with 
respect to CSG trees and for evaluating the boundaries of the corresponding solids are known, at 
least for simple geometric domains. Emerging CAD applications require that we extend the domain 
of solid modelers to support more general and more structured geometric objects. This paper fo- 
cuses on dimensionally non-homogeneous, non-closed pointsets with internal structures. These 
entities are well suited for dealing with mixed-dimensional ("non-manifold") objects in R" that have 
dangling or missing boundary elements, and that may be composed of several regions. A boundary 
representation for such objects has been described elsewhere. In this paper, we propose to specify 
and represent inhomogeneous objects in terms of Constructive Non- Regularized Geometry 
(CNRG) trees that extend the domain of CSG by supporting non-regularized primitive shapes as 
leaves and by providing more general set-theoretic and topological operators at interior nodes. We 
also provide filtering operations that construct CNRG objects from aggregates of selected regions 
of other CNRG objects. We present a syntax and semantics of the operators in CNRG and outline 
some basic algorithms for classifying pointsets with respect to the regions of objects represented by 
CNRG trees. 

Keywords: Geometric Modeling, Solid Modeling, Constructive Solid Geometry, Inhomogeneous 
Objects, Non-Rcgularizcd Pointsets, Non-Manifold Topology, Internal Structures, Algorithms, Set 
Membership Classification. 
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I.O INTRODUCTION 

Solid modeling is concerned with unambiguous computer representations for physical solid objects, 
and with algorithms that process such representations. Earlier CAD systems represented solids by 
collections of edges, known as wireframes, which sometimes correspond to several solids— i.e., they 
are ambiguous or incomplete— or to no solid at all— i.e, they are invalid. Much of the initial impetus 
of solid modeling theory and technology was motivated by the desire to ensure the geometric 
completeness and validity of representations and the correctness of algorithms. Solids were modeled 
mathematically by r-sets, which are bounded, closed, regular, semi-analytic subsets of [1]. Some 
authors considered a more restrictive modeling domain consisting only of r-sets whose boundaries 
arc 2-manifolds [2,3]. (Informally, a 2-manifold is a collection of 2-D faces— for example, 
triangles— such that each edge is shared only by two faces and each vertex is the apex of only one 
cone of faces [4,5].) Although all solid modelers internally manipulate lower-dimensional geometric 
entities such as curves and surfaces, typically a user cannot access directly such entities, or use them 
to construct other non-solid objects. (Some of the earlier solid modelers such as BUILD [6] also 
accommodated I-D and 2-D objects, called shelb and wires, whereas others such as PADL-2 
[7] provided procedural facilities for manipulating such entities only as part of a solid's boundary.) 
The most useful schemes for representing solids have been, and still are, Constructive Solid Ge- 
ometry (CSG) and Boundary Representations (BReps) [1]. A CSG representation is a tree whose 
internal nodes represent either Boolean operations or rigid motions, and whose terminal nodes 
correspond to primitive solids such as blocks, cylinders, and tori. A BRep is a graph with nodes 
corresponding to the faces, edges and vertices that define a solid's topolo^cal boundary. Much has 
been said (correctly and incorrectly,..) about the virtues of each of these representations. We, and 
many other geometric modeling practitioners, believe that both representations are useful and 
complementary. For example, a BRep is excellent for graphic user interaction, whereas CSG is very 
concise and easily supports parameterized object families. Current solid modelers t3rpically contain 
both CSG and BReps in some form, and this trend is likely to continue indefinitely. 
Many CAD/CAM applications now being developed require a geometric domain that extends be- 
yond solids. The two following examples illustrate some of the requirements. Other examples may 
be found in [8,9, 10, 11] and range from feature-based design to stereolithography and fmite-element 
methods. 

Contact relationships are crucial in kinematics, assembly planning, grasping, and related areas. How 
can one tell if two solid objects A and B are in contact? What is their region of contact? Contact 
exists if the boundaries of the solids intersect and their interiors are disjoint. Alternatively, contact 
may be formalized by requiring that A and B intersect but that their regularized intersection be 
empty. (Recall that the regularized intersection is the topological closure of the interior of the 
standard intersection [I].) The region of contact is simply the intersection of the topological 
boundaries of A and B. This example shows that contact calculations involve regularized and 
non-regularized Boolean operations, as well as topological operations such as boundary and inte- 
rior. Observe that regularized operations may be expressed in terms of their non-regularized 
counterparts plus closure and interior. Therefore, a modeler that supports the standard 
Boolean— intersection, union, diiference and complement— and topological— interior, closure and 
boundary— operations on solids and lower-dimensional geometric entities is sufficiently powerful for 
studying contacts between objects. 

Many practically-important objects such as composite-material aircraft parts and semiconductor 
circuits are composed of several regions with different material properties. Observe that the object 
in Figure 1 (left) cannot be represented in a solid modeler as the union of its regions, because the 
union operator obliterates the internal boundaries between the different regions. Furthermore, the 
subtraction operation performed in Figure 1 (center) preserves the structure of the object. This is 
not the case for Boolean operations commonly supported in solid modelers. I'inally, the decom- 
position of the object's complement. Figure 1 (right), is available for validity checks. This example 
shows that geometric modelers should accommodate structured objects composed of several re- 
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gions, and that operators that manipulate simultaneously all the rc^ons of a structured object are 
very useful. 




Figure 1. Operations on structures: The four-layered object (left), A, is represented as an 
aggregate of regions. A hole through the three top layers of A is obtained by 
subtracting an object B (center). The decomposition of the model into layers has 
been preserved during the subtraction operation. Furthermore, a decomposition 
of the subtracted region B according to the layers of A is available (right) for 
verification. For instance, the intersection of B with A4 is lower-dimensional, 
which indicates that the surface of the region A* has been exposed, but that the 
region itself has not been modified. 



Several definitions have been proposed to extend the notion of solids to encompass more general 
geometric models. Many are ba^ on mathematical definitions of simplicial complexes or related 
concepts which support pointsets that need not be dimensionally homogeneous. These extensions 
seek to represent, in the same model: a wireframe (i.e., a set of curves), a collection of surface ele- 
ments (isolated faces), an aggregate of volumes, and even higher-dimensional sets. 

Several boundary-based data-structures for representing "Non-Manifold Topolo^cs" have been 
proposed in [12,13] for 2D objects, in [14,15,16,17,18] for 3D objects, and in [19,20,21] for 
higher-dimensional objects. These schemes are in general limited to the boundaries and internal 
structures of closed pointsets that are collections of vertices, lines .segments, polygonal faces, 
polyhedra, and their higher-dimensional counterparts. 

A representation based on "s-sets", i.e., finite aggregates of disjoint, open regularized (and thus 
full-dimensional) cells was proposed in [22]. These cells are equal to the interior of their topological 
closure. The set of such s-sets is closed under Boolean intersection, but not under the standard 
Boolean union. Indeed, the union of two s-sets needs not be an s-set (see Figure 2 and also Figure 
8 of [22] for example). Consequently, the construction of s-scts models of assemblies may use 
intersection, difference, and "gluing" operations (i.e., operations that perform a union of cells whose 
interiors are disjoint), but fails when two objects with overlapping cells must be combined through 
a union operation. Arbab therefore defines a regularized version of the union operator, which re- 
turns the interior of the topological closure of the pointset produced by the union. Unfortunately, 
when applied to an aggregate of cells, this operator destroys all internal structures (i.e., removes 
faces that separate cells). 

A more general scheme for representing arbitrarily-dimensioned non-closed pointsets that have 
internal structures and external "dangling" boundary elements was reported in [23], and is based 
on the concept of Selective Geometric Complexes (SGCs). An SGC object is an aggregate of mu- 
tually disjoint cells that arc connected, dimcnsionaily-homogcncous, rclativcly-open, singularity-free 
subsets of algebraic varieties. The vertices, edges, faces, and volume elements typically represented 
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A B C D 

Figure 2. S-sets: The union of the two s-sets, A and B, produces an open pointset, C, with 
an interior crack. C is not an s-set, i.e. the union of disjoint cells that are equal 
to the interior of their closure. Arbab's regularized union produces the 
corresponding s-set, D, by "fdJing in" the crack. 



in solid modellers, and their counterpart in higher dimensions, are cells. With each cell of an SGC 
is associated a flag, which specifies whether the cell is active or not. The pointset represented by an 
SGC is the union of (the pointsets of) its active cells. Furthermore, SGC cells may be tagged with 
attributes. The union of all active cells with the same attribute values forms a region. An SGC may 
therefore be decomposed into several re^ons, each region representing a dimensionally non- 
homogeneous pointset that need not be closed. High-level algorithms for performing Boolean and 
topological operations on SGC structures are briefly described in [23]. These algorithms maintain 
the decomposition of SGC objects into regions, while merging them through Boolean and other 
operations. 

SGC representations are extensions of boundary graphs, and are produced by performing sequences 
of user-specified construction steps that transform or combine SGC representations of other ob- 
jects. Editing commands, resultmg geometric features, and arguments to algorithms that compute 
properties of the resulting structures are best specified in terms of individual regions, or their 
topologically-distinguishable subsets, which, for practical purposes, must be unambiguously defined 
in terms of the original users' input (for instance, references to construction steps or to primitive 
shapes) [24,25]. 

To compare the topological domains of geometric modellers one should examine the topological 
constraints imposed on the pointsets supported by these modellers and on their internal decom- 
positions. Figure 3 shows characteristic examples of increasing topological domains. Only two- 
dimensional examples are used for clarity. 

• A is an r-set with manifold boundary. (Each vertex is adjacent to exactly two edges.) Note that 
manifolds need not be composed of simply connected regions (one component of A has a 
hole), but the various connected components of manifolds must be totally disjoint. 

• B is an r-set with non-manifold boundary. (A vertex is bounding more than two edges.) Note 
that regions whose interior are disjoint may share common vertices, but not common edges. 
Both A and B are equal to the closure of their interior. 

• C is an s-set, (i.e. union of disjoint open regularized regions). C is composed of two disjoint 
regions that share a common boundary edge. Tliis common edge is called "interior" because 
it is in the part of the boundary of C that lies in the interior of the closure of C. Obviously, 
closed sets do not have interior edges. S-scts cannot have non-separating interior boundary 
elements, such as interior vertices and cracks (i.e., interior boundary edges that do not separate 
two disjoint regions, but arc connected to the same region on both sides). 

• Non-separating interior boundary elements are possible with open sets, such as D, which do 
not contain their boundaries. Note that the boundary of D is non-manifold, i.e, there is no 
restriction on the number of edges bounded by any particular vertex. 

• While the sets A through D are all homogeneously two-dimensional, Fi is not. It has a dangling 
edge that is not bounding any two-dimensional region. 
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• F is a one-cUmensional non-manifold boundary. Note that F has components that do not 
bound any finite two-«Umensional region. Topologies, such as F, have been associated with 
the term "non-manifold" boundaries. Note that F implicitly defines two bounded two- 
dimensional regions, Fj and Fj. The term of "Non-Manifold Topology" has been used to 
popularize such topological domains in R' [26], where sets like F would be considered as the 
boundary of a non-manifold set. However, the union of F with F| and Fj is a closed set shown 
m E, whose boundary does not include all of F. Therefore, F cannot be the "boundary" of a 
two-dimensional pointset that would include these two regions. 

• G is the union of an open two-dimensional set with a subset of its boundary, which is not a 
manifold (as imbedded in the two-dimensional space) and is not even a one-dimensional r-set 
(some edges are in F but their bounding vertices are not). 

• 11 illustrates a general non-regularized pointset. It is a collection of disjoint connected open 
cells of dimensions 0, 1, and 2. This type of pointsets subsumes all previous types (A to G) 
illustrated in this figure. H exhibits cracks, interior vertices, and dangling edges and vertices that 
do not bound any rc^on. Note however, that to produce a crack or a separating edge, the 
corresponding edge has been removed from the pointset, so that it can be included in the 
boundary. 

• I is a closed non-manifold structure that subdivides the underlying pointset into two two- 
dimensional regions, having interior boundaries (iilled cracks and points) that are distinguished 
from the rest of the two-dimensional regions because they are explicitly listed as boundary el- 
ements. Objects such as I are commonly associated with the category of "non-manifold 
topologies", because the topological notion of a boundary has been relaxed to accommodate 
interior boundaries. Note that this type of decomposition is a generalization of simplicial 
complexes, as discussed in [23]. 

• J is an example of a Selective Geometric Complex [23], i.e.. a collection of disjoint relatively 
open cells— here, two-dimensional re^ons, edges, and vertices. J is more general than I, be- 
cause it needs not be closed. It is also more general than II, because it is a decomposition of 
space, and not the representation of a pointset and its boundary. Cells, as defined in [23], are 
required to be connected and dimensionaily-homogeneous. Such cells correspond to a de- 
composition of pointsets that is too fine for practical interaction with users. In the present 
paper, we describe a coarser decomposition, where cells are grouped into non-homogeneous 
regions according to their origin in the design process. 

In this paper, we propose a representation scheme that covers a rich set of inhomogeneous ge- 
ometric objects and operations. We call it Constructive Non-Regularized Geometry (CNRG). 
CNRG trees represent objects that are aggregates (i.e., unions) of mutually disjoint regions. Each 
region is a pointset in R" and needs not be connected, regular, or even dimensionally homogeneous. 
The leaves of a CNRG tree correspond to parameterized primitive shapes such as volumes, faces, 
curve segments, or points. Internal nodes correspond to intermediate CNRG objects and are asso- 
ciated with topological and Boolean operations. 

SGCs provide a model for representing non-regularized internally-segmented pointsets as a col- 
lection of connected open cells defined recursively in terms of their boundaries. CNRG trees yield 
an alternate representation for these pointsets as a collection of regions defined in terms of original 
primitive pointsets. Re^ons of CNRG objects need not be open nor connected and typically cor- 
respond to unions of SGC cells of various dimensions. CNRG models should be viewed as a pri- 
mary model for user interaction because they support a high level vocabulary for expressing 
operations and regions, and because the CNRG trees are, similarly to CSG trees, easy to edit and 
archive. SGC models should be automatically derived from CNRG representations, as boundary 
representations arc derived from CSG trees. 
The proposed scheme extends CSG in several ways: 

1. It uses standard (non-regularized) Boolean operations and topolo^cal operations— boundary, 
closure and interior. The regularized Boolean operations can be implemented in this scheme 
as three-operator sequences: standard Boolean, followed by interior and closure. 
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Figure 3. Topological domains: Two-dimensional pointsets are shown as dashed regions. 

Their bounding edges and vertices are shown as thin rectangles and small bullets, 
respectively. A is an r-set with manifold boundary. B is a general r-set. C is an 
s-set. D is an open set with non-manifold boundary. E is a closed set with 
non-manifold boundary. F is a one-dimcnsional non-manifold boundary. G is 
the union of an open two-dimensional set with a non-regularized subset of its 
boundary. H is a non-regularized pointset. I is not a pointset, but a closed 
non-manifold structure. J is a Selective Geometric Complex. 



2. It admits as primitives non-solid objects such as points, curves, or surfaces, and higher- 
dimensional objects. 

3. It mtroduces a new opoator, called aggregation, which constructs structured objects composed 
of several regions. The aggregation operator is a formally-defined and more sophisticated ver- 
sion of the "assembly" operator provided by modelers such as PADL-2. Structured objects 
are not pointsets. They are collections of point sets, much like the cell complexes of algebraic 
topology [4]. Structured objects, also called in this paper CNRG objects, or simply objects, 
also have underlying pointsets, as cell complexes do. The underlying set of a CNRG object is 
the union of all the regions of the object, and tliercfore it is a pointset with no structure or 
"internal boundaries". 

4. It defmes Boolean and topological operations on structured objects. These operations do not 
extend the domain of the scheme, since they can be expressed in terms of standard Boolcans, 
topological operators, and segregation. But they arc very convenient because they manipulate 
simultaneouiily all the regions of their structured-object arguments. The definitions ensure that 
the underiying sets are treated consistently. For example, the union operator for structured 
objects produces another structured object whose underlying set is the union of the underiying 
sets of its arguments. 

The remainder of this paper is organized as follow.-i. Section 2 contains the formal definitions of the 
CNRG operators. Section .1 develops representations for the individual regions of structured ob- 
jects. Section 4 discusses some of the fundamental algorithms for CNRG, and shows that they are 
very similar to their CSG counterparts. (A full-fledged suite of CNRG algorithms is still under 
development, and is outside the .scope of this paper.) The final section summarizes the paper. 
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2.0 DEFINITIONS AND NOTATION 



A CNRG object is a set of pairwise disjoint pointsets of R" called regions. 

A region is a possibly non-regular (semi-analytic) subset of R". Although a region of a CNRG 
object is a pointset and should be manipulated as a whole without farther decomposition, for 
practical purposes it may be internally represented in some modellers as the union of primitive ge- 
ometric elements, such as vertices, edges, faces, and other "cells" as defined iji [23]. 
The pointset, pA, of a CNRG object, A, is the xmion of the pointsets of its regions. 
Thus a CNRG representation not only defines a pointset, but also a decomposition (internal 
structure) of the pointset. Regions, and thus objects, need not be dimensionally-homogeneous, 
closed, or connected. 

A CNRG tree is a rooted directed acyclic graph that represents a CNRG object. Leaves of the tree 
are CNRG primitives, which may be composed of more thjm one region. For simplicity, we assume 
that each primitive is given in its absolute position, although in practice, rigid body transformation 
nodes may be used. Internal nodes represent intermediate CNRG objects obtained by applying 
Boolean, topological, simplification, or filtering operators to the CNRG objects represented by 
their child-nodes. 

To simplify the following discussion, we use "|" to denote a "gluing" or aggregation Wnary operator 
that takes two disjoint legions, aggregates of regions (all of which are pairwise disjoint), or disjoint 
CNRG objects and produces a CNRG object that aggregates all the regions. Note that e:q>ressions 
involving only aggregate operators are order-independent. 

In this paper, an uppercase letter denotes a CNRG object, and each of the object's regions is de- 
noted by the same letter with a subscript. For example, A, and A2 are two regions of the same 
CNRG object, A. If A is composed of only these two regions, wc write, A = {AilAj}. By defmition, 
A| and Aj are disjoint for i 9^ j. A single-region object, A, is considered distinct from its region A,. 
We write A = {A,}. Note that A, = pA for single region objects. 

In the following, we assume that S, A, B, C and D are CNRG objects. Furthermore, we often use 
C or S to denote the objects produced by applying to A (and to B) a Boolean or topolo^cal CNRG 
operator. 

A and B are disjoint, if and only if the intersection of their pointsets, pA and pB, is empty. 
Given two regions. A, and B,, A, n B, denotes their intersection and A|\Bj their difTerence in the 
standard set-theoretical sense. Note that A, n B, and A,\Bj are sin^e re©ons that may be empty, 
discoimected, dimensionally inhomogeneous, and not closed. 

A region is said to be contained in an object if it is contained in the pointset of the object. It need 
not correspond to the union of any subset of the object's regions. 

We now define a set of new operators for CNRG objects. These operators are called: union, 
intersection, difference, complement, interior, closure, boundary, and rcgularization. They are an- 
alogues to the standard pointset operators. We use a different notation to distinguish between the 
two set of operators. For example, we use "n", "U", and "\" for pointset intersection, union, and 
difference, and we use " + '\ and for their CNRG counterparts. The topological comple- 
ment of a pointset S is denoted S. The corresponding CNRG operator will be denoted "c". 
The names of these operators were chosen according to their effect on the pointsets spanned by the 
objects and not on the internal decomposition of these pointsets. We have chosen to include the 
CNRG version of the standard topological operators (complement, interior, closure, and boundary) 
because these operators offer a well accepted and convenient syntax for specifying pointsets and for 
analyzing their topological properties. Our Boolean operators (union, intersection, and difference) 
correspond to the standard operators used in CSG. Except for the symmetric difference of A and 
B, which can be expressed as (A - B) -f (B - A), these three operators cover all non trivial bounded 
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Boolean combinations of two pointscts, A and B. Unbounded combinations may be expressed in 
terms of their bounded complements and of the complement operator. 

The aggregate operator "I" and the pointset operator "p" were introduced in this paper to simplify 

the definitions. They need not be made available for use in CNRG trees. Instead, the union and 

simplification operators, precisely defined below, should be used. 

A filter operator that eliminates selected re^ons will be introduced in the next section. 

We define below the semantics for each of the CNRG operators. 

2.1 Aggregation 

Given n disjoint pointsets, A„ the aggregation operation, denoted "|", creates the corresponding 
CNRG object: A = {A, I Aj I . . . |A„}. 

2.2 Simplification 

The simplification, sA, of A is a CNRG object with a single re^on: the pobtset pA. Thus, 
sA = {pA}. 

2.3 Complement 

The complement, cA, of A is an object composed of a single region that is the pointset complement 
of pA. Hence, cA = {pA}. 

2.4 Union 

The union, A + B, is an aggregate of regions of the following three types: A| 0 Bj, A,\B, and Bj\A, 
for all combinations of regions, A,, of A and regions, Bj, of B. 

Potentially each region of A is split into two sets: the part in B and the part outside of B. The 
second set is a single region. The first set may be decomposed according to the decomposition of 
B into regions (see Figure 4 for example). 

Union produces a subdivision of (pA) U (pB) that is compatible with the decomposition of A and 
B into regions (see [23] for a more detailed discussion of subdivision). We call it "union" because 
the pomtset p(A + B) equals the set theoretic union, pA U pB. 
We discuss in Section 3 a technique for distingiMshing each region of A + B. 

2.5 Intersection 

The intersection, A*B, of A and B is the aggregate of all regions of the type Ai n Bj. 

A region, A,, of A is truncated to A, fl pB, and is subdivided according to the subdivision of B into 

regions. Consequently, p(A*B) = (pA) n (pB), which justifies the name for this operator. 

2.6 Difference 

The difference, A - B, is the aggregate of all regions of the type A|\pB. 

Clearly, p(A - B) = (pA)\(pB). It follows from the definition of union, intersection, and difference 
that A + B = {(A - B) I (A*B) | (B - A)}. 

2.7 Interior 

The topolo^cal interior, iA, of A in R" is the aggregate of regions. A, fl interior(pA),. that are the 
intersection of the regions of A with the topological interior of A in R". 

Note that, although iA is always full-dimcnsional, regions of iA need not be full-dimensional. 'ITius 
the interior operator returns a subset of the original object and preserves its internal decomposition 
into regions. The pointset, piA, of the interior of A equals the topological interior of pA. 
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Figure 4. Union: A single-region object (left), A, is unioncd with an object (center), B, 
that has two regions. In the resulting object (right), S = A + B, the unique region 
Ai of A is split into 3 regions: Si? = A,\B, which is disconnected, S31 = A, n B,, 
and S33 = A| n Bj. Each one of the two re^ons of B is split into two parts: their 
respective overlap with A (S31 and S3J ) and the parts out of A (Sj, = B,\A and 
Sjj = B2\A). The lo^c of our the subscript notation will be discussed in Section 
3.10. 



2.8 Closure 

The closure, kA, of A, is the aggregate composed of all the regions of A plus a single new region 
defined as the difference between the topological closure of pA and pA itself. Thus, we can write: 
k{A„ A„} = (A A„, I (closure(pA)\pA)}, and pkA equals the topological closure of pA. 

2.9 Boundary 

The topological boundary, dA, of A is defined as: 3A = kA - iA. The boundary operator does not 
necessarily return a single-region object. Note that p(3A) equals the topological boundary of pA. 

2.10 Regularization 

The regularized version, rA, of A is defined as kiA. The pointsct, prA, spanned by rA corresponds 
to a regularized solid as defined in [27] and equals rpA. 

Note that regularization does not imply amplification (i.e., a regularized object may be composed 
of many non-regularized re^ons). However, taking the boundary, C, of a regularized object, A, 
will produce a lower-dimensionaJ object whose pointset, pC, is the topolo^cal boundary of the 
pointset of the regularized object. 

2.11 An example 

The effect of interior, closure, and boundary are illustrated ui Figure 5. 



3.0 IDENTIFYING AND SELECTING REGIONS 



A CNRG expression combines primitive regions through the operators defined above. A CNRG tree 
is a graph structure obtained by parsing the CNRG expression and generating an internal node for 
each operator and a leaf for each primitive object. 
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Figure 5. Operations: The top row shows a CNRG object, A, composed of two regions, 
Ai and Aj, shown on the same line. The next row shows B, the interior of A, and 
its two constituent regions. The third row shows C the boundary of A, (labelled 
"bA" in the figure) and its three constituent regions: a subset of A,, a subset of 
A,, and a new region, outside of A. The bottom row shows D, the closure of A, 
and its three constituent r^ons: A,, A,, and the same new repon as in C. Note 
that p(lcA) = p(iA) U p(9A). The subscript notation, called the region's signature, 
and our convention for constructing it is introduced in section 3.10. 



Union, intersection, difference, closure, interior, and boundary operations on objects, when fol- 
lowed by simplification, behave as their set theoretic counterparts applied to pointsets. However, 
without simplification, they return aggregates of pairwise disjoint regions. Therefore, when 
"evaluating" a CNRG tree or classifying a pointset with respect to a CNRG object, we need to take 
into account the object's decomposition into regions. 

A CNRG tree defines an object, i.e., a pointset and its decomposition. But how can individual 
regions be identified? IIow can one test points for inclusion in a specific region? These questions 
will be addressed below. 

First consider a simple example. Let A, U, C and I) be single-region objects. The expression 
(A + B)*(C + D) produces an object composed of possibly as many as 9 regions, some of which 
may be empty. We provide a .set theoretic formulation for these 9 regions below. (For simplicity 
of notation, let the "|" operator have the lowest priority in CNRG expressions). 
From the definition of union, A + B = (A - B)|(A*B)|(B - A). Iixpanding the union in the object's 
expression yields 

(A + B)*(C + D) = ((A - B)|(A*B)|(B - A)) ♦ ((C - D)|(C*D)|(D - Q). 
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Distributing "♦" over 'T'l produces the following 9 regions: 

(A + B)*(C + D) = (A - B)*(C - D) i (A - B)*(C*D) | (A - B)*(D - C) | 
(A ♦ B)*(C - D) I (A ♦ B)*(C ♦ D) I (A ♦ B)*(D - C) | 
(B - A)*(C - D) I (B - A)*(C ♦ D) | (B - A)*(D - C). 

For generjJity, we formulated each region in terms of operations on objects, but, because the 
primitive objects are single-re^on objects, each region has an equivalent set-theoretic expression. 
For example, (A - B)*(C - D) is the singje region: (pA\pB) n (pC\pD). 

How many regions does a CNRG have? The answer depends on the pointsets (it may have zero 
regions, because all the re^ons may be empty), but we can still enumerate all the possible regions 
and count their maximum number. 

Suppose that a particular repon of A (and a pajriiculay re^on of B for the case of binary operators) 
can be selected. Are these selections sufficient to unambiguously define any particular region of C? 
The answer depends on the operator. As shown below, for some operators, selecting a region of 
A (and a region of B) is insufficient to define a specific region of C. Unambiguous region selection 
may be achieved by associating with operators a semantic interpretation, which determines a specific 
group or "type" of regions produced by the operation. The possible interpretations are listed below 
for each operator. 

Any region of C can be unambiguously defined by: (1) a region of the left child, A, of C, (2) a re- 
gion of the right child, B, of C (if it has one), and (3) a semantic interpretation for the operator 
associated with C- If the primitives (the leaves of the CNRG tree) are single-region objects, then 
any particular region of the resulting CNRG object may be specified by selecting a semantic inter- 
pretation for the operators in the CNRG tree. Depending on the CNRG expression, some se- 
lections may not be necessary. For example, s(A + B) has only one region, regardless of the 
interpretation of 

In the following sub-sections we assume that A and B are composed of m and n regions respectively 
and that C corresponds to a singile-operator expression involving a left child, A (and possibly a rigiht 
child, B, for CNRG operators that take two arguments). 

3.1 Simplification 

The simplified version, sA, of A has only one repon, and thus simplification has only one inter- 
pretation. No selection of a region of A is needed for specifying the re^on of sA. 

3.2 Complement 

The complement A of A has also only one re^on. These is no need for selecting a region of A, nor 
an interpretation of complementation. 

3.3 Union 

The union A + B has (n + m + nm) regions of three different types. Kach type corresponds to a 
different interpretation of " + ", which wc distinguish by associating the corresponding subscript 
number with the " + " sign; 

1. Interpretation A -f-,B yields the m sub.sets of the m regions of A that lie out of pB. 

2. Interpretation A +2B yields the n subsets of the n regions of B that lie out of pA. 

3. Interpretation A +3B yields the (nm) regions of type A, n Bj, that lie in both A and B. 

Any region of A + B may be selected by specifying the appropriate interpretation of " + " and by 
selecting a rc^on of A and/or of B. (For example, when the first interpretation is used, the spec- 
ification of a region in B is not necessary and will not change the result.) 

Expressions involving union operators may define several regions. If these expressions combine 
single-region primitives, any particular rcgjon of the resulting CNRCJ object may be defined pre- 
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cisely by selecting the interpretation (i.e. type 1, 2, or 3, as described above) of some or all of the 
union operators. 

3.4 Intersection 

The intersection A*B has (nm) regions of the form: Aj fl Bj. Any one of these re^ons may be se- 
lected by specifying a region of A and a region of B. Intersection has only one interpretation. 

Note that A*B and A +3B refer to the same thing, (i.e., identical pointset and identical decompos- 
itions). However, one uses "♦" and " + " operators to specify pointsets and their decompositions. 
A +3B is a convenient notation used in this paper to introduce a way of identifying regions of 
A +• B. As discussed in Section 3.10, an aggregate of the subscripts of the operators of a CNRG 
expression may be used to identify a region or an aggregate of regions. 

3.5 Difference 

The difference, A — B, has m regions of the type Ai\pB. Any one of these regions may be selected 
by specifying the corresponding region of A. There is only one interpretation for the difference 
operator and there is no need for specifying a region in B. 

Note that A - B and A +,B define identical objects. 

3.6 Interior 

The interior, iA, of A has at most m regions that are subsets of the corresponding regions of A. 
Any one of these regions may be selected by specifying the corresponding region of A. There is thus 
only one interpretation for the interior operator. 

3.7 Closure 

The closure, kA, of A has {m+ 1) regions: the m original re^ons of A, plus one new region. Any 
one of the first m regions may be specified by selecting the corresponding region of A and by se- 
lecting the "ki" interpretation of closure. The last region is obtained by selecting the "kj" interpre- 
tation. 

Note that p(kjA) = clo5ure(pA)\pA. Wc have decided to consider all the pieces of k2A as part of a 
single region, independently of tfieir relations to re^ons of A, because some pieces of kjA may si- 
multaneously bound several regions of A (see for example D2? of Figure S on page 9). Note that 
the region kjA needs not be fuU-dimensional or closed. 

For example, given two single-region primitives, A and B, the expression k(A + B) defines four 
regions specified as follows: k,(A -f-jB), k|(A +2B), k|(A +3B), and kjCA -t- B), which is independent 
of the interpretation of 

3.8 Boundary 

5A has at most (m + 1) regions. The m bounding regions that arc subsets of rc^ons of A correspond 
to the "9," interpretation of "d". The only re^on not contained in A corresponds to the "di " in- 
terpretation of "5". 

Note that the operator interpretations "62 " and "K " arc equivalent. 

The effects of interior, boundary, and closure, and the resulting regions are shown in Figure 5 on 
page 9. 

3.9 Reguiarization 

rA also has (m+ 1) regions. The m regions contained in A correspond to the "r," interpretation, 
and the region outside of A corresponds to the "rj" interpretation of "r". 
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3.10 Region signature and sub-objects 

The subscript notation for the semantic interpretation may not be practical for users. It should ci- 
ther be hidden or replaced by a more intuitive notation in a modeler's user interface. But the no- 
tation is well suited for internal use in a modeller and has the advantage of providing a compact 
specification for a region of a CNRG object. 

For example, let S be defined in terms of single-region primitives as (A*(B + C)) + D. S has seven 
regions shown in Figure 6. Each one can be associated with a particular interpretation of the " + " 
operators. We define the r^on's signature as the concatenation of the interpretation numbers that 
correspond to the region, in the order in which their operators appear in the expression of the 
CNRG object. (We use numbers only for operators that have multiple interpretations.) We use the 
regions signature as a subscript of the object's name to defme the corresponding region. For in- 
stance, S„ = (A*(B +,C)) -t-|D. 

Note the use the wild card, "?", in Figure 6 when the semantic interpretation of a particular oper- 
ator is irrelevant. For example, the region 8,3 does not depend on the interpretation of the fust " + ", 
since it lies outside of p(A 4- B). 

A re^on signature is complete when it unambiguously defines a single region, i.e. when all relevant 
semantic operators have been specified. An incomplete signature may be obtained by replacing by 
"?", in a complete signature, one or more semantic operator numbers. An incomplete signature is 
a filter that selects regions and defines a new CNRG sub-object decomposed into regions that cor- 
respond to all the possible combinations of the values that the "?" wild cards can take. 
A sub-object is a valid CNRG object and may be used as such in CNRG expressions. 

The wild card mechanism provides a simple notation for specifying sub-objects. For example, S,i 
in Figure 6 is composed of the three regions: 

{ (A*(B -H.O) +,D I (A*(B-f,C))-|-,D I (A*(B -K,q) H-,D }. 

For any sub-object, T, of a CNRG object, S, we have pTcpS. Furthermore, any region of T is a 
region of S. In particular Sn . . . ? is equivalent to S. 

For a CNRG object defined in terms of multiple-re^ons primitives, a region signature concatenates 
(1) the numbers associated with the operators' semantic interpretation and (2) the numbers associ- 
ated with the primitives' regions. Both types of numbers are sorted according to the ordered of the 
corresponding symbols in the object's expression. Figure 4 on page 8 shows another example of 
the use of subscripts for selecting re^ons of primitive objects. 

A filter operator, f,ig„„„„ , is associated with a valid signature of its CNRG object argument. For 
example, f,ip„tureA returns— as a CNRG object— the sub-object A5ij„„ure. Results of filter operations 
on the same object may be combined using the union operator to produce filters that cannot be 
expressed with a single signature of the object. 

3.11 Constructive representation for regions 

The CNRG operators are defined in Section 2 by specifying their eflccts on the regions of the ar- 
guments. All of the operators produce regions that can be expressed in terms of standard set the- 
oretic and topological operations on the regions of the operands. 

For example, the re^ons of C = A - B arc of form C, = Ai\pB and pB = [JB,. Therefore, each C, 
can be expanded in terms of Aj and Bj through standard set theoretic diffebnce and union. If A 
and B arc themselves composite CNRG objects, wo can represent their regions constructively, and 
apply recursion until wc reach the primitives. It follows that any region of a composite object can 
be represented constructively in terms of standard set-theoretic and topological operations on the 
regions of the primitives in the object's CNRG tree. The constructive representation for a region 
is easy to obtain from the CNRG tree and the interpretation attached to some of the operators. 
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Figure 6. Signatures: The expression S = (A*(B + C)) + D, yields 7 regions: S,i, 8,3, Sj„ 
S». Sji, S33, and Stj. 



The existence of constructive representations for regions is an important property of CNRG. 
Constructive representations for regions aie similar to CSG and can be processed by similar algo- 
rithms. Whereas CSG algorithms deal only with regularized set operations on solids, CNRG algo- 
rithmis must process non-solid regions through non-regularized set operations as well as topological 
operations. Thus, CSG algorithms caimot be used directly for CNRG, but must be suitably mod- 
ified. Work on CNRG algorithms is under way. Preliminary results are reported in the next section. 



3.12 Regions in the complement 

For some applications, such as the analysis of the validity of design features [11], regions outside 
the pointset of the object (for instance slots removed from a part model) and their interactions may 
be of interest. 

Consider the object S = (A - B) - C., of Figure 7, obtained by subtracting from the primitive A 
two slots, B and C. Although regions of B -J- C are outside of S, they must be analyzed to derive 
the topolopcal relation between the two slots. For example, the existence of an intersection 
pA n pB n pC of B and C inside A (Figure 7 a and b) indicates that, in order to reduce manufac- 
turing costs, both features should probably not be machined independently. 

Most outside re^ons of interest are bounded, lie inside the union of the primitives' pointsets, and 
can be formulated as CNRG expressions that combine the ori^nal primitives. 
In particular, replacing specific "-" and "♦" operators of the CNRG tree with " + " operators 
produces a useful decomposition of the original object and of the relevant portion of its comple- 
ment. 

Let S be the original CNRG object. The operators in the tree of S are numbered according to the 
order in which they appear in the CNRG expression of S. \xX I be the number of a "-" operator 
and J the number of a "♦" operator in S. I^t S' be the CNRG object obtained by replacing in S 
the I and J operators by " + ". S is the sub-object of S' defined by using in the signature of S' a 
"1" in the I'" position, a "3" in the J"" position, and a "?" in all other positions. Indeed, is 
equivalent to "-" and "+3" is equivalent to In the example of Figure 7, where 

S = (A - B) - C, we can define S' to be (A + B) + C, and obtain S = S'„. 

Using filters that combine other interpretations for operators I and J yields sub-objects of S' that 
are outside of S. Regions of these sub-objects may be constructed through filters and unions, as 
described in Subsection 3.10. Their topological relations with other regions may often be expressed 
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in terms of the existence and dimensionality of regions of CNRG objects obtained by applying 
topological and Boolean CNRG operators to these sub-objects. 

For example, the existence and maximum dimension of selected regions in these sub-objects indi- 
cates contact or interference between features that have been subtracted by intersection or differ- 
ence. In the above example, the dimensionality of S'33 indicates whether the portions of the two 
features, B and C, that lie inside A overlap (see Figure 7a), share a common face (see Figure 7b), 
or are disjoint (see Figure 7c). 

4.0 ALGORITHMS FOR CNRG TREES 

The most fundamental tool for interrogating a CNRG tree is an algorithm that classifies a point 
with respect to the CNRG object represented by the tree. 

On a standard CSG tree, this operation is called PMC (for Point Membership Classification) [28] 
and returns a scalar that can take three values: IN, ON, or OUT, depending on whether the point 
is in the interior, the boundary, or the complement of the represented CSG solid. For points that 
lie on the boundary of several primitives, the result of PMC (with respect to sub-solids) is not only 
one of these three values, but also the neighborhood of the classified point with respect to the 
pointset represented by the corresponding node of the CSG tree [28,29]. The neighborhood is fiJl 
if the point is in the topological interior of the solid or eny}ly if the point is outside the solid. For 
points on the boundary of the solid, the neighborhood representation is more complex. The 
neighborhood of points that lie in the (2-D) interior of a solid's face may be represented by the 
surface that contains the face and by a single bit of information specifying whether the solid is on 
one side of the surface or on the other. For points that lie on an edge of the solid, but not on a 
vertex, an arrangement of faces around the edge, plus relative orientation information defmes the 
neighborhood. Neighborhoods must be combined according to the CSG e:q)ression. 
Since a CNRG object is typically composed of several regions, the notion of PMC must be ex- 
tended to distinguish these regions. As e;^lained earlier, any individual region of a CNRG object 
may be selected by attaching semantic interpretations to some of the operators of a CNRG tree. 
PMC (and other functions) may be computed independently for each region, or for all regions in 
parallel. As in CSG, PMC may be done by classifying the point with respect to the primitives and 
merging the results of these classifications upward in the CNRG tree. We do not discuss here 
pomt/primitive classification, and simply assume that it is available for all the primitives of the 
domam. (This is clearly the case for primitives defined as algebraic inequalities.) Assuming that 
PMC results (with respect to a particular region) are known for the operands (children) of a CNRG 
node, we define below how these results should be comUned for each interpretation of the CNRG 
operators. We limit the discussion to two cases: interior points, when points are in general position 
(not lying on any curve, face, ... hyperface), and hyperface points, when points lie on a single 
hyperface. For points on more than one hyperface, one needs n-dimensional neighborhood re- 
presentations and techniques for performing Boolean operations on such neigliborhoods. 
Neighborhood representation for CNRG are more complex than for CSG because CNRG objects 
may be non-regular sets. Observe that a point on the boundary of a regular set must be adjacent 
to a portion of the object's interior, but this need not be true for CNRG objects. The neighbor- 
hood representation must be able to distinguish between points that are ON the boundary and 
belong to the object and other points that are ON the boundary but do not belong to the pointset. 

4.1 Interior points 

Consider a point that docs not lie on any boundary (with respect to R") of any primitive. Such a 
point must be either in the topological interior of a region or in the interior of the object's com- 
plement. Therefore, the point's classification has only two pos.sible values that distinguish whether 
the point belongs or not to the region. We represent these values, respectively, by Boolean true and 
false values. PMC for such points may be computed as follows. Ixt Ka (respectively K„) be the 
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Figure 7. Outside regions: Two slots B and C have been subtracted from A. The 
corresponding objects is defined by the CNRG expression: S = (A — B) — C- 
Different situations may be characterized by the dimensionality of S'33, where 
S' = (A + B) + C. In Figure a, S'33 is fuU-dimensional, which indicates that part 
of feature C will be removed while machining B. In Figure b, S'33 is 
lower-dimensional, which indicates that the two features are connected and should 
be considered as a compound feature. In Figure c, S'jj is empty and both features 
may be machined independently. 



classification of a point P with respect to a particular region of A (respectively B). These regions 
have been selected by specifying, when necessary, the appropriate semantic interpretations for op- 
erators in the sub-trees of A and B. Also let Kp^ (respectively Kpo) denote the classification of P 
with respect to the pointset, i.e. union of all the regions, of A (respectively B). Let Kc be the clas- 
sification of the point with respect to a region of the parent node that corresponds to an operation 
on A, or on A and B. can be expressed as a Boolean formula in Ka, Kp^, Kg, Kps- We treat each 
interpretation of each operator as a separate case: 

Case A: Kc = not Kp^ 

CasesA: Kc = KpA 

CaseiA: Kc = Ka 

CasekiA: Kc = Ka 

Case kjA: Kc = false 

Case 3iA: Kc = false 

Case djA: Kc = false 

CaseA*B: Kc = Ka and Kr 

Case A - B: Kc = K^ and not Kp„ 

Case A +, B: Kc = K^ and not Kp„ 

Case A +j B: Kc = K„ and not KpA 

Case A -{-j B: Kc = Ka and Kb 

These formulae follow directly from the definitions of the operators. Note that classification of an 
interior point with respect to a region defined by kj, 3,, or 5j is always false. 



4.2 Hyperface points 

I^t us now consider the case of a hyperface point, P, i.e., a point that lies on a single 
(n-l)-dimensional hyperface of a primitive or possibly on the (n- 1 ) -dimensional intersection of se- 
veral overiapping primitive hyperfaccs. P does not lie on any (transversal) intersection of two 
primitive hyperfaccs that are not subsets of the same hypcrsurfacc. Wc also assume that P does not 
lie on any cusp or singularity of the hypcrsurface. Ixjt 11 denote the hypersurface that contains the 
primitive hyperfacc(s) containing P. 
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The open n-dimensional neighborhood ball, N, of a hyperface-point, P, is divided into three sets 
by the hypersurface H: Nj„H, the (n-l)-dimensional intersection of N with II, and the two open 
disjoint sets, and Nr of N\H, one on each side of H. Because P is a hyperface point, the radius 
of the ball can be chosen such that these three sets are disjoint from any other hypersurface that 
bounds the primitives of the CNRG tree. Consequently, ^ven any region of the CNRG object, 
each of these three sets is either entirely inside or entirely outside of it. Therefore, given P and II, 
the neighborhood of P with respect to any CNRG region may be represented by three "flags" (bits 
of information), each indicating whether the corresponding parts of B (N,„h> N^, and Nr) lie inside 
or outside that region. 

Assuming that H is consistently oriented, we refer to the two sides of II as left and right. The 
Boolean symbols K^c and Krc define whether there is material (of that particular region) around 
the classified point. K^c refers to the left side of the hypersurface and Krc to the right side. For 
example, K^c = true indicates that the region of C under consideration includes a hypervolume on 
the left side of H, in the neighborhood of P. 

The symbols Kc, Ka, Kb, Kp^, and Kpg are defmed as above. Furthermore, Kla. Kra, Klb, and 
Krb specify whether there is material of the selected regions of A and B around P on the corre- 
sponding side of H. The example of Figure 8 illustrates the meaning of these symbols. 

Klpa. Krpa, Klpb. and Krpb have analogous meanings, but with respect to the entire pointsets of 
A and B. 

Case sA: Kc = KpA, K^c = K,pA, Krc = Krpa- 

Case A: Kc = not KpA, Klc = not Ki^a. Krc = not Krpa- 

Case A - B: Kc = Ka and not KpB, K^c = K^a and not K^pB, Krc = Kra and not Krpb- 
Case iA: Kc = Ka and K^pA and Krpa. K,c = K,,a, Krc = Krpa. 
Case A*B: Kc = Ka and Kb, Klc = K^a and K^e, Krc = Kra and Krb. 
Case A +, B: Kc = Ka and not Kpe, Klc = K^a and not Klpd, Krc = Kra and not Krpb. 
Case A +j B: Kc = Kb and not KpA, K^c = Klb and not Klpa, Krc = Krb and not Krp^. 
Case A +3 B: Kc = Ka and Kg, K^c = K^a and K^b, Krc = Kra and Krb- 



Case 5, A: 
Case dik: 
Case k,A: 
CasekjA: 



Ic = Ka and not (Klpa and Krja), K,x = Klpa, Krc = Krpa. 
Kc = not Ka and (Klpa or Kr^a). Kj^ = K^,a. Krc = K^^. 
Kc = Ka, Klc = Kla> Krc = Kr^. 
Kc = not Ka and (K^a or Kra), Klc = K^pA, Krc = Krpa- 



5.0 CONCLUSION 



Dimensionally inhomogeneous pointsets of R" having missing boundary elements and with mtemal 
structures may be defmed using a set of new Boolean and topological operators that construct and 
maintain a structural decomposition of pointsets into aggregates of disjoint regions which are non- 
regularized pointsets. These operators do not have all of the properties of their standard counter- 
parts (for example, the union does not distribute over intersection). But they treat the pointsets that 
underlie the structured objects in a manner consistent with the corresponding standard operations. 
For example, p(A -I- B) = pA U pB and p((A + B)*C) = p((A*C) + (B*C)). 
The specification of inhomogeneous geometric objects may be stored in a CNRG (Constructive 
Non- Regularized Geometry) tree. Geometric algorithms ori^nally developed for processing CSG 
trees in three-dimensions can be adapted for CNRG trees. 

Since CNRG objects are aggregates of mutually disjoint regions, algorithms for processing CNRG 
trees must be able to distinguish among these rc^ons. Any particular region of a CNRG object 
may be specified by its signature, i.e., a selection of the semantic interpretations for some operators 
of the CNRG expression. The signature indicates how the corresponding region can be expressed 
in terms of standard set-theoretic and topological operations applied to primitive pointsets. 
Sub-objects, i.e. aggregates of regions, may be specified using a wild card in the object signature. 
These sub-objects are CNRG objects and thus can be used in Boolean and topological operations 
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Figure 8. Neighborhoods: The object A + B„ combines the two single-region objects, A 
and B. The point P lies on a bounding hyperface conunon to A and B. Let C 
denote the region A +,B. The neighborhoods of P with respect to A, B, and C 
are shown using a small disk around P. For A (left), K/v = true, 
Ki-A = Klpa = true, ajid Kra = Krpa = false. For B (center), Kb = true, 
Ki,B = Klpb = false, and Krb == Krpb = true. Consequently, for d, Kc = false, 
Klc = true, and Krc = false. 



and for region selections. Sub-objects may be obtained through a filter operation parameterized by 
the object's signature. Union of sub-objects suffice to express any combination of an object's re- 
^ons. 

A point can be classified with respect to a region of a CNRG object by a recursive descent, divide 
and conquer algorithm, which classifies the point with respect to the CNRG primitives (i.e., the 
leaves of the tree), and combines the results according to the CNRG expression of each region. 
Formulae for combining point classifications are presented in this paper for points that do not lie 
on the transversal intersection of two or more primitive-bounding hypersurfaces. 
Explicit, boundary-like representations of CNRG objects are supported by Selective Geometric 
Complexes (SGC) structures and other "non-maiufold" representation models. 
Algorithms for evaluating SGC representations incrementally from a CNRG tree have been pro- 
posed [23] as compositions of operations that subdivide the cells (i.e. connected manifold com- 
ponents of the regions) of two objects, select and mark the appropriate components, and/or simplify 
the structure by deleting unmarked components and merging marked ones when appropriate. An 
implementation of these algorithms is currently under way. 

Algorithms that directly evaluate boundary elements of the individual regions of 3-D CNRG ob- 
jects and that perform Set Membership Classification with respect to regions are also under study. 
They involve more elaborate versions of the classification procedures discussed in this paper. 
In summary, the CNRG methods discussed in this paper extend CSG repre.sentations and algo- 
rithms beyond solids to the domain of inhomogcncous objects which have "non-manifold 
topologies". 
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Lambda Research's Latest Version 
of TracePro 



Lambda Research Corporation released Version 2.0 in 
September '99. This release raytraces, on the average, 
10 times faster than its previous version. Further feature 
enhancements include a wizard to convert raw BSDF 
data into TracePro format (later article in this newsletter 
describes in more detail), over 3000 catalog lenses, a 
new bulb catalog, STL import capability, new expanded 
material catalogs and an optional healing husk to fix 
imperfect geometry. 

In additional incremental updates since then, features 
have added specifically for the lighting design and 
biomedical markets. These features include bulk 
scattering for biomedical applications, updated 
luminaire design and analysis, a ray generation tool for 
bitmap images and a more powerful user interface. 

The new features for the medical equipment market 
include ray tracing volume scatter in human tissue and 
general bulk scatter through any volume. TracePro 2.0.4 
provides enhanced support of standard illumination 
output formats for the Illumination Engineering Society 
and Eulumdat, and DXF output for irradiance and 
illuminance maps. The new user interface displays 
additional data on-screen for multiple or single 
wavelengths at the user's discretion. And the optional 
bitmap module wizard translates any standard BMP, 
JPEG, GIF, or PNG file into a TracePro source file 
which can be ray-traced through any geometry and then 
displayed on any surface as an irradiance or illuminance 
output. 




January 2000 



In the first quarter of 2000, TracePro Version 2.1 will be 
released featureing thin film stacks and ray tracing 
through gradient index materials. 




What one of our customer's has to say 

about Version 2.0 - Fred Bushroe, President of 
Inov Inc., with over 14 years experience building optical 
hardware and developing products and systems, 
commented - "TracePro 2.0 is incredibly fast compared 
to TracePro 1.4. I saw a speed increase of 26 times, 
raytracing a project containing reflective conic 
surfaces. Analysis that took an hour to trace now takes 
about two minutes. What's great about TracePro is the 
well thought out User Interface, and the fact that I 
routinely import and export complex solid models with 
ease. lean bring a radical 3D model or assembly in 
from SolidWorks, and focus my efforts on tracing rays 
and getting results, rather than spending days wrestling 
with model translation issues. TracePro's user interface 
is intuitive. This is important to me because I do many 
different things in my business. Jean come back to the 
pro-am after a month, and get right to work without 
having to releam the user interface. Assigning optical 
properties to lenses and mechanical parts is a breeze. " 



Our Mission Statement: Establish a technol- 
ogy bridge between optical and mechanical 
engineers by providing state-of-the-art opto- 
mechanical design and analysis software, services 
and technical support 
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SolidWorks Partnership 

Lambda Research has just signed on as a Research 
Partner with SolidWorks Corporation. Our objective is 
to improve the data transfer between Solidworks and 
TracePro. 

How do I create a thin sheet using a 
macro in TracePro? 

A thin sheet is a bounded area with zero thickness. In 
TracePro this may be build using the following ACIS 
Solid Topology: edges, wires, sheets and bodies. The 
boundary is build by a series of connecting linear edges 
which fully enclose the desired area. The loop of edges 
must not intersect. 

In the Macro language we build up aa list of edges and 
convert them into a wire body. A wire body is a wire 
frame like solid without any surface data. 

Once we construct the wire body which bounds our 
desired area we cover it with a surface. This becomes a 
single sided solid object. Single sided objects only have 
one surface normal (or can be thought to be infinite in 
thickness). To be useful to TracePro we can make the 
sheet double sided, with surface normals on either side 
of the now infinitely thin object. 

Once constructed the thin sheet may be extruded and 
revolved to give it a non-zero volume. 

The following example will construct a square boundary 
in the XY plane of size X wide and Y high. 

We can assign some variables, X&Yin this case: 
(define X 2) 
(define Y 3) 
(define edge-list 
(list 

(edge:linear (position 0 0 0) (position X 0 0)) 
(edge:linear (position X 0 0) (position X Y 0)) 
(edge:linear (position X Y 0) (position 0 Y 0)) 
(edge:linear (position 0 Y 0) (position 0 0 0)))) 

The boundary is enclosed so we can create the wire 
body: 

(define the-wire-body (wire-body edge-list)) 
Create a sheet body from the wire-body: 

(define the-sheet-body (sheetxover-wires the-wire- 
body)) 

Change the sheet body into a double-sided sheet body: 
(sheet:2d the-sheet-body) 



Done! 

This can now be modifies within TracePro. To change 
the boundary you simply add more edges to the list. 

The macro with minimal comments follows: 



(define X 2) 
(define Y 3) 

; Construct a list of the edges 

(define edge-list 
(list 

(edge:linear (position 0 0 0) (position X 0 0)) 
(edge:linear (position X 0 0) (position X Y 0)) 
(edge:linear (position X Y 0) (position 0 Y 0)) 
(edge:linear (position 0 Y 0) (position 0 0 0)))) 

; The boundary is enclosed so we can create the wire 
body 

(define the-wire-body (wire-body edge-list)) 

; Create a sheet body from the wire-body. 

(define the-sheet-body (sheet:cover-wires the-wire- 
body)) 

; Change the sheet body into a double-sided sheet body 
(sheet:2d the-sheet-body) 

You can find more helpful hints on our web site in 
the TracePro App Notes Section. 

TracePro's BSDF Wizard 

The BSDF Wizard is a small 32-bit Windows utility, 
which can read ASTM files (as generated by a Schmitt 
Measurement Systems Scatterometer) and plot the graph 
of BSDF against p. p is calculated from 0 using the 
Harvey-Shack Linear, Shift-Invariant Representation. A 
synthetic curve can then be matched to the data graph, 
thus giving an approximation of A, B, and g. 

An ASTM file can be opened either by dragging the file 
into the graph rectangle area, or by clicking the "Open" 
menu item in the "File" menu. If there is a signature file 
specified in the ASTM file, a dialog box will ask to 
confirm opening the signature file. 
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The signature file is a second measurement that was 
taken with the instrument in an equivalent configuration, 
with no surface to be measured. This allows the back- 
ground from the instrument to be subtracted out of the 
measurement. If the ASTM file does not specify a 
signature file, or specifies a non-existent or invalid 
signature file, a different signature file can manually be 
loaded by choosing the "Load Signature" menu item in 
the "File" menu. 

There are several view options available in the "Graph" 
menu to aid in matching the synthetic curve to the data. 

The Logarithmic P and Logarithmic BSDF selections 
control how the two axes are displayed in the graph. A 
log-log plot is best for mirror-like surfaces, whereas a 
linear-linear plot is best for rough surfaces. For mixed 
surfaces, such as glossy paint, a semi-log plot is usually 
best. 

There are two display states for the input data: 
connected and non-connected. In the connected state, the 
data is represented as a continuous line. In the non- 
connected state, each data point is represented by a 
diamond shaped marker. The unconnected state helps in 
determining the distribution of the data, the connected 
state is useful for observing how closely the synthetic 
curve matches the data curve. The display state can be 
toggled by selecting "Connected" from the "Graph" 
menu. 

The "Subtract Signature" selection in the "Graph" menu 
can be used to include or ignore the differences the 
signature makes. The show signature button can show 
or hide the signature graph. 

A legend in the lower right hand comer of the BSDF 
Wizard dialog box, which indicates the colors of the 
data, signature, and synthetic curves. 

At the upper and lower right hand side of the graph, 
there are two pairs of up- and down- arrow buttons that 
adjust the range of BSDF that is displayed in the graph 
window. Each pair of buttons increases or decreases the 
upper or lower limits of the graph. The synthetic curve 
can be adjusted to match the data points by dragging the 
two black control points until the sample curve appears 
to fit the data. If the control points are not within the 
viewable portion of the graph, the "Normalize" button 
will adjust the sample curve and graph limits so that all 
data points and control points are visible. For fine- 



tuning, or to compare curves, the roll-off point and "A" 
can be specified manually by entering numbers in the 
edit boxes and clicking the "Apply" button to apply the 
changes to the synthetic curve. As the curve is adjusted, 
the "Integrated Scatter" (lower right hand side of dialog 
box) is calculated to show the total scatter of the 
synthetic curve. 

The preferences can be adjusted by clicking on the 
"Preferences" menu item in the "File" menu. At this 
time, the only option is whether to confirm loading the 
signature file if detected. If this option is disabled, the 
signature file will automatically be loaded if it has a 
valid format. The "Properties" menu item in the "File" 
provides a view of additional information from the 
header in the data file and signature file. 

To create an ASTM file manually: The first line of the 
file may contain anything at the beginning of the line. 
ASTM files generated by the SMS Scatterometer 
contain the file name, and the date and time the file was 
created. If there is a signature file associated with this 
ASTM file, then the first line should end with the words 
"signature file" (in all lowercase) followed by a colon 
and space, then the name of the signature file. The lines 
following the first line may contain comments about the 
contents of the ASTM file. Next, it is necessary to create 
the header block, which contains information about the 
instrument settings at the time of the measurement. The 
only header that is required by the BSDF Wizard is the 
"Incident Angle" field. The incident angle line should 
begin with the exact words "Incident Angle (deg)" 
followed by a colon, space, and the incident angle in 
degrees. The last item in the file, after the header block, 
is the BSDF measurement data. Immediately before the 
data should appear a line reading "Angle Data Follows 
[angle (deg), BSDF]" followed by a colon. The next 
line after this line should begin the data. Each line in 
the data section is one data point. Two numbers 
separated by spaces specify one data point. The first 
number is 8 in degrees, and the second number is the 
BSDF. Scientific notation may be used. 

Here is a minimal ASTM file: 

signature file: SIGNATURE. C02 

This is the comment section, may omit 

Incident angle (deg) : 3 

Angle Data Follows [angle (deg) , BSDF] : 
-92.0000 5.632E-06 
-91.5000 1.428E-05 
-91.0000 1.664E-05 
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Illuminations 



Volume 2 Issue 1 



Another Lambda Product-APEX™ 

APEX is an automatic interferogram analysis program 
for optical testing. It analyzes images from an optical 
interferometer to determine the aberrations present in 
the surface, element, or system under test. Pass/fail 
tolerances can be set in the program for go/no-go 
testing. Detailed reports can also be generated, along 
with graphical presentation of the aberrations, as well as 
Modulation Transfer Function (MTF) and Point Spread 
Function (PSF). 

UPCOMING TRAINING 
Illumination Analysis 
Using TracePro 

March 27-28, 2000 
Two-Day Intensive Course 

Stray Light Analysis 
Using TracePro 
March 29-31,2000 
Three-Day Intensive Course 

Call for registration information (978)486-0766 
Watch for other course dates for Illumination Analysis 
and Stray Light Analysis 



Lambda Research represents the 
SOLEXISy family of products. 

Recently Lambda Research has conunenced reselling the 
database products of coatings and materials compiled by 
Stellar Optics Research International Corporation 
(SORIC) over seven years. This is the world's only data 
resource and selection tool for black, white, reflective 
and transmissive surfaces and materials for ground and 
space based applications. 

Lambda Presented at SPIE Annual 
Meeting in Denver-July 1999 

Our topics were: 

• Translation and Interchange for Mechanical CAD 

and Lens Design Data 

• Polarization Models for Monte Carlo Ray Tracing 

• End to End Electro-Optical Modeling Software 

• Edge diffraction in Monte Carlo Ray Tracing 

If you would like to read the papers in their entirety, go 
to our Website and find them under TechNotes. 



Distributors 

Taiwan: InfoTek Information Systems Co., Inc. 
No. 19, L-128 

Chong-Shan N. Road. Soc. 2 
Taipei, 104, Taiwan 

Contact: Susan Lo (Susan@infotek.com.tw) 
Phone: (81)866 252 35377 
Fax: (81)886 251 12360 

Japan: Future Instruments Trading Inc. 

1-3 Nihonbashi-Ohdenmacho Chuo-Ku 

Tokyo, 103, Japan 
Contact: Tom Yamamoto 

(tsutom@green.ocn.ne.jp) 
Phone: (81)3 366 67100 
Fax: (81) 3 366 70094 

Korea: Modern High Tech Co., Ltd. 

Lifecomby Building, 61-4, Suite 418 
Yeidodong, Yungdungpo-Ku, Seoul, So. Korea 
Contact: H.B. Seo (modhitec@nuri.net) 
Phone: (82)278 60447 
Fax: (82)278 60814 

Europe: LIGHT TEC 

2 avenue Olbius Riquier 

Hyeres, 83400, France 

Contact: Yan Comil (sales@lighttec.fr) 

Phone: (33) 4 94 01 37 47 

Fax: (33) 4 94 01 37 46 



Watch for Upcoming Announcements! 



See Lambda Research at these 
Upcoming Shows in 2000 

Jan. 25-27 - Photonics West, San Jose, CA 

Booth # 1209 
March 6-9 - SAE, Detroit, MI 
March 19-23 - LightFair, Frankfurt, Germany 
April 10-11 - OptoSouthwest, Albuquerque, NM 
May 9-11 - Lightfair, New York, NY 
May 16-18 - SID, Long Beach, CA 
June 27-30 - Optatec, Frankfurt, Germany 
August 1-3 - SPIE Annual Meeting, Denver, CO 
September 10-15 - CLEO, Nice, France 
November 6-8 - Photonics East, Boston, MA 
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Schild et al., "Open Data Exchange with HP PE/SolidDesigner," Hewlett-Packard Journal, 
October 1995, pp. 35-50. 



open Data Exchange with 
HP PE/SolidDesigner 



Surface and solid data can be imported from HP PE/MESO and 
exchanged with systems supporting the IGES, STEP, and ACIS formats. 
Imported data coexists with and can be manipulated like native data. 

by Peter J. Schild, Wolfgang Klemm, Gerhard J. Walz, and Hermann J. Ruess 



HP PE/SolidDesigner supports the coexistence of surface 
data with solid data and provides the ability to import and 
modify surface and solid design data from a variety of CAD 
systems. Backward compatibility with HP PE/ME30 preserves 
the investment of existing HP customers. Using improved 
IGES (Initial Graphics Exchange Standard) import capability, 
both surface and wireframe data can be imported. Surface 
data and solid data can also be imported and exported using 
the STEP (Standard for the Exchange of Product Model Data) 
format. Once imported, this data can coexist with HP PE/ 
SolidDesigner solid data. It can be loaded, saved, positioned, 
attached to, managed as part and assembly structures, de- 
leted, and used to create solids. Attributes such as color can 
be modified. If the set of surfaces is closed, HP PE/SolidDe- 
signer will create a solid from those surfaces automatically. 

HP PE/SolidDesigner 3.0 also allows solid parts and assem- 
blies to be exported to ACIS-based systems using Version 1.5 
of the ACIS® SAT file format. This feature provides a direct 
link to other ACIS-based applications. 

From PE/ME30 to PE/SoIidDesigner 

HP PE/ME30 is a 3D computer-aided design (CAD) system 
based on the Romulus kernel. To preserve the investment 
of existing customers it was required that the transition from 
HP PE/ME30 to HP PE/SolidDesigner be as smooth as pos- 
sible. Therefore, an HP PE/ME30 file import processor is a 
integral component of HP PE/SolidDesiger. 

In HP PE/ME30, 3D objects are built from analytic surfaces 
like cylinders, cones, spheres, planes, and toruses. The inter- 
sections of these surfaces can be represented as explicit 
analytic curves such as straight lines, circles, and ellipses, or 
implicitly by describing the surfaces involved and providing 
an approximation of the intersecting arc. Parabolic and 
hyperbolic intersections are represented implicitly. 

HP PE/ME30 Native FUe Organization 

HP PE/ME30 supports the Romulus textual transmit format. 
The transmit file is not intended to be read by humans but 
the general structure can be examined. The file contains 
only printable charaaers, and real values are represented as 

* A kernel is the heart of a rnodelirig system. Currently, three kernels are 
systems. These are Romulus from Shape Data, Parasolid. an extension of Ronuihjs, and the 
ACIS Kernel from Spatial Technology. 



text strings. The full format of a transmit file consists of six 
different sections. These will be described using the example 
of a single cylinder positioned at the origin of HP PE/ME30's 
coordinate system with base circle radius 10 and height 20. 

The first section, the header section, describes the environ- 
ment, the machine type, the user login of the file creator, 
and the time and date when the model was created. 

©•AOS 

€>- Machine ^pe HP-UX 

@' Transmitted by user j(yz on 27-May-94 at 13-06 

The second section contains index and counting information 
related to the schema described in the third section. The 
schema defines the data structures used to represent the 
objects. It consists of a collection of record definitions. The 
following is an example of a record definition: 

SH 2 FS 3 -1 BK 3 1 One record definition 

tl t Property definition 1 
I Property definition 2 
Number of properly definitions 

' Name of record definition 

The following is an example of a property definition: 

FS 3-1 One properly definition 

I I ^ Length field: If -1: variable length 

' Type field 

' ID name 

In the second section of the transmit file the number of 
record types, the numbers of record instances and property 
instances, the name of the schema, and its version and 
update number are supplied. The record instances and 
property instances contain the concrete data describing the 
model. The semantics and the sequence of data entities have 
to conform to the format specified by the corresponding 
record definition and property definition entities. 

The information in the cylinder example file says that 11 
instances of record definitions are supplied to describe the 
schema for the instance of the cylinder. For the actual object, 
23 record instances built out of 115 property instances are 
used. 
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11 

23115 

ROMDSCHMA74 
16 

The third section, the schema section, contains the definition 
of the data structures used to represent the model. This sec- 
tion consists of the subset of record definitions ftom the 
HP PE/ME30 internal data struaure schema that are needed 
to represent the model. The schema sections of files repre- 
senting different models will be different. The schema sec- 
tion for the example cylinder is: 

BY 19 UP 4 -1 SE 3 -1 TX 5 -36 Fl 4 -6 CI 4 -12 PI 4 -8 
Gl 4 -1 SI 4 -3 Tl 4 -3 

RA 2 1 RN 2 1 Zl 4 -2 FN 1 1 CN 1 1 PN 1 1 IN 1 1 
SN 1 1 ZN 1 1 NM 1 1 
SH 2 FS 3 -1 BK 3 1 

FA 8 UP 4 -1 AK 3 -1 RV 1 1 SF 3 1 SX 2 1 VR 3 -1 
HA 2 -3 SL 3 1 

VR 4 PT 3 1 BE 3 1 BV 3 1 FE 3 1 

ED 2 CU 3 1 RV 1 1 

CU 3 UP 4 -1 AK 3 -1 TR 3 1 

TR 6 UP 4 -1 AK 3 -1 BK 3 1 EQ 2 -7 IS 3 2 TY 1 1 

PT 4 UP 4 -1 AK 3 -1 C0 2-3GP3 1 

GP 5 UP 4 -1 AK3 -1 BK 3 1 CO 2 -3 PX 3 1 

SF 7 UP -1 SD 3 -3 AK 3 -1 BK 3 1 EQ 2 -7 SU 3 -5 TY 1 1 

UA 3 OW 3 1 CL 1 1 II 1 '1 

The fourth section contains, for each record type defined in 
the schema section, the number of data objects used for the 
transmission of the model. The sequence of numbers is 
identical to the sequence of record definitions used in the 
schema section. In the cylinder example, the object consists 
of one body built of one shell built of three faces. Four ver- 
tices, four coedges, two curve geometries, two edges, two 
points with two geometric point definitions, three surfaces, 
and one attribute are needed to represent the cylinder ob- 
ject. The file contents are: 

1 1 3 4 4 2 2 2 2 3 1 

The fifth section, the data section, contains the data structure 
instances. The contents of all records needed to represent 
the object are found in this section. To every record an inte- 
ger record label is assigned. This number will be used in 
other record instances to point to the instance. In general 
the instances in the file appear in the order in which they 
are referenced by other entities. The data of an entity in- 
stance is not split. If forward references are contained in the 
inistance definition the next instances can be found in 
exactly the same sequence as referenced. Because this rule 
applies recursively, newly referenced entities can be found 
first in the physical file sequence. If all references of an en- 
tity are resolved completely the next reference of the next 
higher level will be resolved. For the cylinder, the data 
section is: 



1125Color1 203 3F04F15F2214E015E1218P019 PI 0 0 0 0.000001 
0.00000000001 0 3 3 2 0 0 0 0 251 1 1 16777215 2 0 1 3 0 0 0 22 0 1 6 0 2 22 
0 00 0 5 00 0 00 -1 01 61810610180 00 20 200003 010 0181014 
014001616 0007 00000 -1 100 02 4 001 23 01 70 2 23 0 0006 
0 0 20 0 0 -10 1 7 1911711 19 0 0 0 2121 0 0 0 3 0 10 20 1911 151 15 0 
0 17 17 0 0 0 7 0 0 20 0 0 -110 0 0 2 5 0 0 0 24 0 2 8 9 0 2 24 0 0 0 0 7 0 0 
0001100281812812121419191391313150 

The sixth and last section contains, for each top-level object 
transmitted in the file, the corresponding root entity and its 
name. In the cylinder example only one objea is trans- 
mitted. HP PE/ME30 supports user-named objects, but in this 
example an HP PE/ME30 default name, BO, has been used 
for the cylinder. 

BO 

Analyzing the Transmit File 

Because the information content of an HP PE/ME30 file can- 
not be understood by simply looking at the file, several 
internal analysis tools are used to extract the information. 
Statistics showing the number of different curve and surface 
types give a first hint of the complexity of the file. A graphi- 
cal presentation of the data instances of a file can be 
generated. 

|_BY |_UA |-> 

l-SH ^-> 

l-FA ^_SF_ 

|_VR I P T I GP 

|_ED LCU LTR 

l-FA USF_ 

|_VR I PT I GP 

UED_|_CU_|_TR_ 

UFA |_SF_ 

LVR 1-> 

LED_|-> 
LVR 1-> 

UED„|-> 

This reference structure can be read easily. The (cylinder) 
object in the file is a body (BY) which consists of a shell 
(SH) and three faces (FA). Shell and faces share the same 
hierarchy level. Each face consists of a reference to a surface 
(SF) and a start vertex (VR). Each start vertex is based on a 
geometric definition of a point (PT) and serves as the anchor 
vertex of an edge loop. A loop is not represented explicitly 
in the HP PE/ME30 exchange file. The implicit connection is 
done by a reference bom a start vertex to the next and pre- 
vious vertices in the loop. The edge (ED) entity represents 
the topological direction of the edge with respea to the 
loop. The curve (CU) entity is an intermediate instance on 
the way to the curve's geometry (TR). 

If complete information from the data section is needed a 
translation tool is available that maps the data section to a 
format much more useful for human readers. The following 
extraa describes how one of the faces and the correspond- 
ing surface component of the cylinder example are repre- 
sented. The mapping from the data section to the readable 
format is also supplied. 
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For this component from the data section: 



the corresponding translated part is: 



= FA (Face owning (anchor) ve 
UP is ElVIPTY 
AKis EMPTY 
RV: INTEGER = 0 
SF: P01I\ITER = 24 
SX: REAL =0 

1 VR: POINTER = 8 

2 VR: POINTER = 9 
HA Is EMPTY 

SL: POINTER = 2 



24 = SF (Surface of face), ttie properties are : 



.. Surface of 
.. Hatching pitch 
.. Anchor of face 
...Anchor of face 
...Hatch direction 
...Shell efface 



List of permanent universal attributes 
Surface supporting this surface 
definition 

Backpointer from element of feature 
Backpointer from assembly or body to 

token 
Geometry definition 
Geometry definition 
Geomeuy definition 
Geometry definition 
Geomeuy definition 
Geometry defintkxi 
Geom^ definition 



REAL = 
REAL = 
REAL = 
REAL = 

5 EQ: REAL = 

6 EQ: REAL = 

7 EQ: REAL = 
SUIs EMPTY 

TY: INTEGER = 2 (CYLINDER) ...Geometry type 
Import Module 

The HP PE/ME30 to HP PE/SolidDesigner import interface 
is linked direcdy to the HP PE/SolidDesigner code. In 
HP PE/SolidDesigner's user interface it simply adds a button 
to the external filing menu. If a file name is specified, the 
processor is activated. Internally, several C++ classes are 
added to HP PE/SolidDesigner to represent the schema and 
instance entities of the HP PE/ME30 file. For every sup- 
ported HP PE/ME30 record definition entity a class derived 
from a generic record instance object is defined. The most 
important member function of each of these classes is the 
convert function. This function performs the mapping of the 
HP PE/ME30 file objea to the corresponding HP PE/Solid- 
Designer entity. 

The three main components of the HP PE/ME30 to HP 
PE/SolidDesigner processor are a lookup table, a schema 
manager, and a set of classes to represent the supported HP 
PE/ME30 file entities. 

The lookup table is part of the interface to an HP PE/ME30 
file. The main task of this table is to manage the mapping of 
HP PE/ME30 file entities to already created corresponding 
HP PE/SolidDesigner entities. A lookup table is generated 

for every open HP PE/ME30 file. 

A schema manager is initialized if a new HP PE/ME30 file is 
opened. It contains the schema section information found in 
the newly opened file. For every open file a corresponding 
schema manager is available to control the interpretation of 
the entities of the file. 



The record instance class builds the third basic data struaure 
of the processor. Record instances are generic containers to 
store all of the data objects that can be expressed by valid 
record definitions. The coristructor of the record instance 
class calculates the entity type from the reference number 
and then allocates memory and reads in the properties from 
the file corresponding to the property definitions of the 
schema. For every supported HP PE/ME30 entity a separate 
C++ class is derived from the record instance class, but the 
generic constructor is used for all subtypes. The main differ- 
endator between the classes is the convert function. 

Conversion Process 

The convert function of the record instance class itself is not 
called by the conversion process. Rather, every derived class 
implements its specific conversion function (in this sense the 
convert function is purely virtual in C++). The individual 
conversion function converts Itself to an HP PE/SolidDe- 
signer entity. 

Conversion and the creation of new derived instances of the 
record definition class constitute a recursive process. If dur- 
ing an active conversion an unresolved (not already con- 
verted) reference is found the corresponding HP PE/ME30 
file entities can be found as the next entities in the physical 
file (see the description of the data section). The conversion 
module then creates a new derived instance of the record 
instance class and forces the translation of this entity to a 
HP PE/SolidDesigner entity that can be used to complete the 
conversion of the current entity. The algorithm is as follows: 
A reference to an HP PE/ME30 file entity is found: 

Already 'converted"? (lookup table search) 
YES: Use the available conversion result 
NO: Create the n&N derived class of record instance 
Call the convert function 
Attach the conversnn result to the k)okup table 
Delete the Instance to free the memoiy used 
Use the newly generated conversion result to continue the 
conversion. 

Nonanalytic Intersection Curves 

The conversion for intersection curves is not done on the fly, 
but by a postprocessor after the rest of a body is converted 
completely. The convert routine for an intersection o-ack 
simply collects the two intersecting surfaces and all available 
additional information found in the file to represent the in- 
tersection. The completion of the intersection curves is done 
by the convert fijncUon for HP PE/ME30 bodies. After a first 
intermediate topology of the new HP PE/-SolidDesigner 
body is calculated and all analytic surfaces and analytic 
curves are attached to the created body, the calculation of 
the intersections begins. 

The topology of the intersection between two surfaces in 
HP PE/SolidDesigner is not always the same as in HP 
PE/ME30 because different constraints on topology and 
geometry exist in the two modelers. For instance, it may be 
necessary to represent the single segment found in HP 
PE/ME30 as a sequence of different curves. In such cases the 
original topology has to be modified and some edges may 
be split. To find the appropriate 
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PE/SolidDesigner is mainly a selection process. In many 
cases two surfaces intersea at not only one but several 
distinct sections. 

Consider the intersection of a cylinder witii a torus in the 
case of perpendicular axes. Four possible intersection curves 
may be part of the model (see Fig. 1). In the HP PE/ME30 
file additional help points are supplied to allow the correa 
selection. The direction of the intersection curve (the tangent 
to the curve) is not guaranteed to be the same in HP PE/ 
SolidDesigner as in HP PE/ME30. Therefore the correct fit to 
the model is calculated and the resulting direction is 
reflected in the topology of the imported model. 

Quality and Performance 

To test the quality of the HP PE/ME30 import processor a 
laige HP PE/ME30 test library has been compiled. It now 
contains more than 2300 examples of parts and assemblies. 
All of the test cases used during HP PE/ME30 development 
and support are included along with new user models con- 
sisting of recently acquired data from internal and extemal 
HP PE/ME30 users. An additional test matrix subtree was 
developed by creating base parts with critical features. In 
particular, all possible surface-to-surface intersections and 
various special cases have been generated. 

The regression test procedure is to import HP PE/ME30 
models from the test library part by part and perform the 
HP PE/SolidDesigner body checker operation on each. The 
loading time and the body checker result are coUeaed in a 
reports file. A reports file can be analyzed by a shell script to 
supply a statistical summary of the current quality of the 
HP PE/ME30 interface. Because of the laige amount of test 
data a complete test takes a long time. Therefore, an inter- 
mediate test is available. The complete test performs the 
basic load and check test on all currently available test mod- 
els of the library directory. The intermediate test examines 
the reports file of the latest complete test and repeats all 
reported problems. It also repeats a random selection of the 
successful tests. At this time over 99% of the complete test 
conversions are classified as successful. 

The performance of the import process for HP PE/ME30 files 
is mainly dependent on three variables: the size of the 
schema, the number of entities, and the number of intersec- 
tions that have to be calculated: 




Fig. 1. Intersection of a torus and a cylinder. 
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Load Time = Size x kl + Entities x k2 + Intersections x k3, 

where kl < k2 < < k3. The size of the schema section 
does not vary very much between different files and is nor- 
mally relatively small compared to the size of the data sec- 
tion. The number of entities and the file size are strongly 
related. The calculation and selection of the nonanalytic in- 
tersection curves fitting the model is a relatively expensive 
component of the processor because a completely new rep- 
resentation of the data structure has to be generated. 

Data Exchange Using IGES 

An important task in computer-aided design is the transfer of 
the completed model to downstream applications and other 
CAD applications. These applications vary from finite element 
analysis and numerically controlled (NC) manufacturing to 
visualization and simulation. HP PE/SolidDesigner currently 
uses IGES 5.1 (Initial Graphics Exchange Specification) for 
file-based data exchange. 

Because of the broad variety of receiving systems an IGES 
interface must be flexible so that the contents of the output 
file match the capabilities of the receiving system. It must be 
possible to transfer whole assemblies keeping the informa- 
tion on the parts tree, or only specific parts of a model, or 
even single curves or surfaces. This is achieved by a mixture 
of configuration and selection mechanisms. 

An analysis of the IGES translators of many different systems 
showed that it is possible to classify them in four main cate- 
gories: 

Wireframe Systems. These systems are only capable of im- 
porting curve geometry. This is typical for older CAD sys- 
tems or 2D systems with limited 3D capabilities. 
Surface Systems Using Untrimmed Surfaces. These systems 
are capable of importing untrimmed surfaces and indepen- 
dent curve geometry. This is typical for low-end NC systems 
that need a lot of interaction to create tool paths and define 

Surface Systems Using Parametrically Trimmed Surfaces. 
These systems are able to handle trimmed surfaces. Trim- 
ming is performed in the parametric domain of the surfaces. 
Periodic surfaces are often not handled or are incorrectly 
handled. Each surface is handled independently. This is typ- 
ical for surface modelers and sophisticated NC .systems. 
Topological Surface Systems and Solid Modelers. These sys- 
tems are able to handle trimmed surftices using 3D curves as 
trimming curves. They are able to handle periodic surfaces, 
nonplanar topology, and surface singularities. Connection 
between adjacent trimmed surfaces is maintained and the 
normal to the trimmed surface is important for inside/out- 
side decisions. This is typical for advanced surface and solid 
modelers. 

HP PE/SolidDesigner's IGES interface is designed to work in 
four output modes: wireframe, untrimmed, trimmed paramet- 
ric, and trimmed. Each output mode represents one of the 
categories of receiving IGES translators. This has the advan- 
tage of giving as much information about the solid model as 
possible to high-end systems (trimmed, trimmed parametric), 
without burdening low-end interfaces with too much infor- 
mation. For some modes (trimmed parametric) more config- 
uration parameters allow fine tuning to specific systems to 



e the transfer rate. Each mode has a specific entity 
mapping that describes which IGES entities are used to de- 
scribe the model (see Tables I, II, and III). Users can specify 
additional product related data and arbitrary comments for 
the start and global sections of the IGES file direcdy via the 
IGES output dialog box. Specific configurations can be saved 
and loaded so that the configuration has to be determined 
only once for each receiving system. Fig. 2 shows the IGES 
dialog menu. 

To allow maximal flexibility in what is translated, the user is 
allowed to select assemblies, parts, faces, and edges and 
arbitrary combinations. All selected items are highlighted and 
the user can use dynamic viewing during the selection 
process. If the user selects assemblies, the part tree is repre- 
sented with IGES entities 308 and 408 (subfigure definition 
and instance). Shared parts are represented by shared 
geometry in the IGES file. 



Table I 
Curve Mapping 



iif 


IGES Settings M 




Output Mode 


♦ Trimmed 


Trimmed Parametric 


Untrimmed 


Wireframe 


Accuracy 


0.001' 


Trimmed Parametric Switches 




if. S&ms .. i^dnaatla Arc 




s lines 


'VX\ 






Coordinate System 


Global 


Local 




Output File Format 


♦ HP-UK 


DOS 




Global Section 


Author 




Organisation 




Sending PID 




Receive PID 




Start Section 






Config File 


, Config Name 


default.cfg' 


Save Conflg 





HP PE/SolidDeslgner 

Straight 
Circle 

B-spline 

Intersection curve 
: curve 



IGES 3D Entity 

Une (110) 

Qrcular arc (100) 
with transformation 

Rational B-spline 

curve (12^ 
Rational B-spline (126) 
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Trimmed Mode 

The trimmed mode is the closest description of the internal 
B-Rep (boundary representation) data stiucture of HP PE/ 
SolidDesigner. It uses tiie IGES bounded surface entities 143 
and 141 as the top element of the model description. Each 
selected face of the part maps to one bounded surface (en- 
tity 143) containing several boundaries (entity 141). Trim- 
ming of the surfaces is performed by 3D model space 
curves. To fulfill the requirements of the IGES specification 
of entities 141 and 143 some minor topological and geomet- 
rical changes of the HP PE/SolidDesigner internal model 
have to be made. Vertex loops are removed, propedges on 
toruses are removed, and intersection curves are replaced by 
B-spline approximations. 

Because the IGES bounded surface entity 143 does not have 
any information about topological face normals, the surfaces 
are oriented so that all geometrical normals point to the out- 
•side of the part (Fig. 3). Thus, enough information is put 
into the IGES file that a receiving system can rebuild a solid 
model fi-om a complete surface model. 

Untrimmed Mode 

The untrimmed mode contains basically the same informa- 
tion as the dimmed mode. For each face the unoimmed 
surface plus all trimming curves are translated. But instead 
of explicltiy trimmming the surfaces with the appropriate 
entities, surface and trimming curves are only logically 
grouped together. This usually requires manual trimming in 
the receiving system, and is only suited for some special 
applications. 




Fig. 3. (left) Solid model, (right) Surface model with normals. 



Kg. 2. HP PE/SolidDesigner IGES output dialog m 
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Trimmed Parametric Mode 

The trimmed parametric mode uses the IGES trimmed para- 
metric surface entity (144) and the curve on parametric sur- 
face entity (142) as representations of a trimmed surface. 
These entities have been established in the IGES standard 
for a longer time than entities 143 and 141 or the trimmed 
mode. For this reason they are more commonly used. The 
main difference fom the trimmed mode is that the trimming 
is performed in the parametric domain of the surfaces. Each 
surface must have a parametric description that maps a point 
from the parameter domain D (a rectangular portion of 2D 
space) to 3D model space: 

S(u,v) = (X(u,v), Y(u,v), Z(u,v)) for each (u,v) in D. 

D = {all (U,v) with Umin S u < Umax, V^in ^ V S VmaJ. 

The following conditions apply to D; 

There is a continuous normal vector in D. 

There is a one-to-one mapping from D to 3D space. 

There are no singular points in D. 

Furthermore, tiimming curves in 2D space must form closed 
loops, and there must be exactly one outer boundary loop 
and optionally several inner boundary loops (holes). Fig. 4 
illustrat&s parameter space trimming. 

These restrictions make it dear that there will be two prob- 
lem areas when converting HP PE/SolidDesigner parts to a 
parametric trimmed surface model: periodic surfaces and 
surface singularities. 



On full periodic surfaces like cylinders, HP PE/SolidDesigner 
usually creates cylindrical topology. There will not necessar- 
ily be exactly one outer loop. Furthermore, 3D edges can 
run over the surface seam (the start of the period) without 
restriction. This leads to the situation that one edge may 
have more than one parametric curve (p-curve) associated 
with it. Also the p-curve loops may not be closed even if the 
respective 3D loop is closed. Fig. 5 illustrates this situation. 

HP PE/SolidDesigner avoids this problem by splitting peri- 
odic surfaces along the seam and its antiseam. The seam and 
antiseam are the isoparametric curves along the parameters 
Umin and Umin+Uperiod/2. Thus, one face may result in 




Fig. 4. Trimniing in parameter space (p-space). 
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Fig. 5. Cylinder topology in 3D and p-space. 

two or four parametrically trinnmed surfaces (u- and v- 
parametric surfaces (toruses)) in the IGES model. Fig. 6 illus- 
trates this situation. 

Another problem with parametric trimmed surfaces are sur- 
face singularities. Singular points are points where the sur- 
face derivatives and normal are not well-defined. For such 
points there is not always a one-to-one mapping from 2D 
parameter space to 3D model space. This means there is an 
infinite set of (u,v) points in parameter space that result in 
the same 3D model space point. Such singularities are easily 
created by rotating profiles around an axis where the profile 
touches the axis. Examples are cones, spheres, degenerated 
toruses, triangular spline patches, and so on (see Fig. 7). 

HP PE/SolidDesigner is designed to handle singularities as a 
valid component of a model. They are marked with a vertex 
if they are part of a regular loop or with a special vertex 
loop if they are isolated from the remaining loops. However, 
it is not possible to express singularities in trimmed paramet- 
ric surfeces legally in IGES. 

To resolve this issue we reduce the singularity problem to 
the problem of the valid representation of triangular sur- 
faces. The splitting algorithm just described is applied so that 
all singularities are pan of a regular loop. Thus, we are 
always faced with the situation illusti-ated in Fig. 8. 

Each singularity of a face is touched by two edges, one 
entering and one leaving the singular vertex. Knowing how 



Hg. 6. Periodic surfaces in 3D and p-space after splitting. 

triangular surfaces are handled in potential receiving sys- 
tems, we offer four ways to export this kind of geometry. 
These are the four possible combinations of closed or open 
parameter loops and avoiding or using singularities. 

Some systems do not need closed p-space loops, while 
others stricdy expect them. If the closed option is chosen, 
the endings of the p-curves are simply connected with a 
straight line. 

Geometrical algorithms usually become unstable near singu- 
larities. Some systems are not prepared to handle this situa- 
tion and will fail. To avoid this, it is possible to shorten the 
parameter curves when entering or leaving a singular vertex 
and connect them at a numerically safe distance. This dis- 
tance is measured in 3D space and is also configurable. It 
usually varies between 0.1 and 0.001. This will result in a 
surface where the region around the singularity is cut out. 
Fig. 9 illustrates the four possible singularity representations. 

WireCrameMode 

For the wireframe mode HP PE/SolidDesigner also avoids 
the cylindrical topology, because in some cases information 
about shape would be lost (e.g., a full surface of revolution). 
After applying the face splitting algorithm all edges of the 
selected faces and parts are translated. No surface informa- 
tion is contained in the resulting IGES file. 




Fig. 7. Examples of surface 
singularities in parameter 
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Fig. 8. Triangular surface situation. 

Extracting Solid Information from Sur£!ice Models 

IGES surface data from solid modelers often contains all sur- 
faces of a closed volume or a connected fece set, However, 
the connectivity between adjacent faces is lost. If the surfece 
model fuUfiUs some specific requirements it is possible for 
the receiving system to recompute this missing information. 
The following describes these requirements and shows how 
connectivity between faces can be reestablished. This 
method can be used to create a solid model from HP PE/ 
SolidDesigner IGES output. 

Automatic comparison of all boundary curves on coinci- 
dence or reverse coincidence would be a very time-consum- 
ing and numerically unstable task. However, it is common 
for the endpoints of the trimming curves of adjacent faces to 
be coincident within a very small accuracy. This makes it 
possible to identify trimming curves that share common start 
points and endpoints. If the two faces of these trimming 
curves have the same orientation one can try to connect the 
faces to a face set. For this task one must try to find a geom- 
etry for a common edge that fulfills the following accuracy 
constraints (see Fig. 10): 
The curve is close enough to surface 1. 
The curve is close enough to surface 2. 
The curve is close enough to curve 1. 
The curve is close enough to curve 2. 

The first candidates for such a curve are the original trim- 
ming curves, curve 1 and curve 2. If either satisfies all four 



requirements it is incorporated into both face descriptions 
and the connection is established. If neither curve can be 
used, one can try a combination of the two, or reduce the 

receiving system's accuracy. 

This method fails if the face orientation is inconsistent or if 
adjacent faces do not share common start points and end- 
points. 

Importing IGES Wireframe Data 

IGES wireframe data can be easily imported into HP PE/- 
SolidDesigner, since HP PE/SolidDesigner's kernel supports 
wire bodies. The modified wire data can be saved in HP 
PE/SolidDesigner's data format. Possible uses for this capa- 
bility include migration from old-line systems to HP PE/Sol- 
idDesigner, interaction with different sources and suppliers, 
and communication with manufacturers. 

In HP PE/SolidDesigner a wire is defined as a set of edges 
connected by common vertices. A body consisting only of 
wires is called a wire body. IGES 3D curve data is used to 
generate the edges of a wire body. This includes lines, 
circles, B-splines, polylines, and composite lines. IGES sur- 
face data such as trimming curves of trimmed surfaces are 
also used to generate edges. To simplify later solid model 
generation the axis and generatrix of a surface of revolution 
are also transformed into edges for the wire body. Since 
only edges have to be generated for a wire body, there are 
no accuracy problems as described above for IGES surface 
importation. On the other hand, information on B-spline 
surfaces is lost. 

Wire data imported from an IGES file is collected into an 
assembly. The assembly gets the name of the IGES file. Any 
substructure of the IGES file like grouping in levels is trans- 
formed into parts within the assembly. Thus, hierarchical 
information contained in the IGES files is maintained within 
HP PE/SolidDesigner, The generated parts can be handled 
like any other part in HP PE/SolidDesigner. To distinguish 
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Fig. 9. Four possible singularity 
representations. 
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Fig. 10. Finding a common edge for adjacent faces. 



wire parts, they can be colored. The options of HP PE/Solid- 
Designer's show menu work for the parts as well as the 
settings of the part container. A wire part can become the 
active part. The edges and vertices of a wire part are display- 
able, all browsers work with wire parts, and wire parts can 
be moved or become members of an assembly. 

To build a solid model from a wire body, the edges and 
vertices of the wire body can be used to position a work- 
plane. Then edges of the wire body can be selected and 
projected onto the workplane. The resulting profile can then 
be used to create a solid, for example by measuring an edge 
length needed for an extrude operation. 

Fig. 11 shows an example of an IGES wireframe mode! with 
four parts and the resulting solid model. Automatic genera- 
tion of solids from wires could be implemented but freeform 
surface information would probably be lost. The real benefit 
of wireframe import is for reference purposes. 

STEP-Based Product Data Exchange 

Manufacturing industries use a variety of national and indus- 
trial standards for product data exchange. These include 
IGES for drawing and surface exchange (international), VDA- 
FS for surface exchange (mainly the European automotive 
industry), and SET for drawing and surface exchange 
(France and the European Airbus industry). This variety of 
different incompatible standards causes a lot of rework and 
waste of valuable product development time which cannot 
be afforded if companies are to survive in the competitive 
marketplaces of tomorrow. Today's standards, originated in 



the eariy 1980s, are no longer satisfactory for product data 
description and exchange. Standards like IGES or VDA-FS, 
which are limited to surface or engineering drawing ex- 
change, do not adequately handle other explicit product 
data categories such as product structure or assemblies or 
geometric solid models. 

Industry trends today are charaaerized by internationaliza- 
tion of manufacturing plants which are spread over the con- 
tinents of the globe, and by lean production in which many 
parts are subcontracted or bought from local or international 
suppliers. National standards and incompatibilities between 
existing standards are obstacles to these trends and will have 
to be replaced by international standards. 

Laige companies in the aerospace and automotive industries 
in the U.S.A. and Europe have now taken the offensive to- 
wards the implementation and use of STEP (5ifandard for the 
Exchange of Product Model Data) as an international stan- 
dard for product data exchange and access, starting in 1994. 
Companies such as BMW, Boeing, Bosch, General Motors, 
General Electric, Daimler-Benz, Pratt&Whitney, Rolls Royce, 
Siemens, and Volkswagen have been using STEP prototype 
implementations in pilot projects with promising results. 

Ultimately, STEP is expected to meet the following require- 
ments for an international product data exchange standard: 
Provides computer interpretable and standardized neutral 
product model data. Neutral implies compatibility with any 
CAD or CIM system that best fits the design or manufactur- 
ing task. 

Implements the master model concept for product data. The 
entire set of product data for a product with many single 
parts is kept in one logical master model which makes it 
possible to regenerate the product as a whole at a new 
manufacturing site. This means that product assemblies, 
including administrative data and bills of material, are 
handled. 

Provides completeness, conciseness, and consistency. This 
requires special data checking and validation mechanisms. 
Provides exchangeable produa data without loss. The 
product data must be exchangeable from one CAD or CIM 
system to another without loss of data. 




Fig. 11. Imported wire-body 
and the solid model constructed 
by HP PE/SolidDesigner. 
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• Provides long-term neutral data storage and interpretability. 
Product data is an important asset of a manufacturing com- 
pany. The piodua data should be retrievable and interpret- 
able by any CAD or CIM system after a long period of time, 
say 10 years or more. This is a significant challenge. 

These requirements cannot be satisfied immediately. The 
STEP program also has shorter-term priorities for standardiz- 
ing specific subsets of the produa data. These include: 

• The complete 3D geometric shape in the form of a 3D 
boundary representation solid model (B-Rep solids) 

• Surface model and wirefi:ame model data 

• Product structure and configuration data. 

Another priority is product documentation. An important 
goal is consistency of the engineering drawing with the 3D 
product geometry. 

STEP Overview 

STEP, the Standard for the Exchange of Produa Model Data, 
is the ISO 10303 standard. It covers all produa data catego- 
ries that are relevant for the product life cycle in industrial 
use. STEP describes product data in a computer interpretable 
data description language called Express. Hie STEP standard is 
organized in logically distina sections and is grouped into 
separate parts numbered 10303-xxx (see Fig. 12). 



The resource parts of the standard describe the fundamental 
data and product categories and are grouped in the Ix, 2x, 
3x, and 4x series. The Express data description language is 
defined in part 11. AH other product description parts use 
the Express language to specify the product data characteris- 
tics in the form of entities and attributes. In addition to the 
product description parts there are implementation resources 
which are given in part 21, the STEP product data encoding 
scheme (the STEP file), and part 22, the Standard Data 
Access Interface (SDAT), which provides a procedural 
method for accessing the produa data. There are different 
language bindings for part 22, such as C or C++ program- 
ming languages. The 3x series parts specify conformance 
requirements for STEP implementations. 

Examples of STEP-standard resource parts are the funda- 
mentals of product description and support (part 41), the 
geometrical shape (part 42), the product structure (part 44), 
material (part 45), the product presentation (part 46), toler- 
ances (part 47), and form features (part 48). The application- 
spedflc resources are grouped in the Ixx series. Examples 
are drafting resources (part 101), electrical (part 103), finite 
element analysis (part 104), and kinematics (part 105). On 
top of the resource parts and application resources are the 
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Fig. 12. Architecture of ISO 10303, Standard for the Exchange of Produa Model Data (STEP). 
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application protocols (AP) which use the underlying re- 
sources in a specific application context, such as mechanical 
design for discrete part manufacturing, and interpret the re- 
source entities in the application-specific context. STEP im- 
plementations for CAD or other computer-aided systems are 
based on application protocols. Application protocols are 
under definition for application areas like basic drafting, 
associative drafting, mechanical design, electrical design, 
shipbuilding, piping, architecture, and others. Here, we high- 
light just two examples, AP203 and AP214. 

AP203: Configuration-Controlled 3D Design. AP203 was 
developed under the leadership of PDES Inc. It covers the 
major requirements for U.S.-based industries such as the 
aerospace industry for government and industrial manufac- 
turing contracts. The produa data covered in AP203 includes 
geometric shape (B-Rep solid models, surface models, wire- 
frame models), product structure, and configuration manage- 
ment. AP203 is die underlying STEP specification for many 
CAD and CIM system implementations. 

AP214: Core Data for Automotive Mechanical Design. 

AP214 has been developed by the automotive industry and 
covers product data categories relevant for the design and 
manufacturing of automotive parts and products. AP214, 
initiated in Gemiany and internationally supported, is still 
under finalization in parallel with its industrial implementa- 
tion in CAD and CIM systems. The implementations have 
been coordinated and harmonized in the European ProSTEP 
consortium and the implementation is focused initially on 
the geometrical product descriptions (solid models, surfece 
models) and product structure. However, all other kinds of 
product data categories relevant for mechanical design in the 
automotive industry (e.g., form features, materials, toler- 
ances) are within the scope of AP214 and are going through 
the standardization process. 

Initial Release 

The initial release of STEP parts focuses on the most 
uigently needed kernel definitions of the standard, which 
cover the geometrical shape description, including all topo- 
logical information, the product structure, and the configura- 
tion management data. Basic product documentation in the 
form of low-level engineering drawings is also covered. The 
parts included in the initial release are parts 1, 11, 21, 31, 41, 
42, 43, 44, 46, 101, 201, and 203. The first two application 
protocols to become standards are AP201: Explicit Drafting 
and AP203: Configuration-Controlled 3D Design. 

Upcoming releases of STEP will cover the next priorities in 
the area of drafting, such as AP202: Associative Drafting, 
materials, tolerances, form features, and parametrics, and 
other application protocols such as AP204: Mechanical De- 
sign Using B-Rep Solid Models and AP214; Core Data for 
Automotive Mechanical Design. 

HP Involvement in STEP 

HP has been working on the standardization of product 
model data since 1989 and has focused on the emerging in- 
ternational standard STEP for 3D product data. The produa 
data focus has been on 3D kernel design data, completeness 
of topology and geometry, B-Rep solid models, and product 
structure and assemblies, as well as on associative drafting 
documentation. HP is an active member in organizations that 



have an impact on the ISO STEP standard, and contributes to 
STEP through national standards organizations in the U.S.A. 
(e.g., NIST, ANSI) and Europe (e.g., DIN In Germany). Of 
particular interest are the oiganizations PDES Inc., PRODEX, 
and ProSTEP. 

PDES inc. HP has concentrated on three major areas of 
PDES Inc's STEP activities: mechanical design of 3D product 
data, associative drafting for CAD data, and electronic data 
definition and exchange. 

The mechanical design initiative of the U.S. aerospace and 
aircraft industries, the automotive industry, and the computer 
industry resulted in STEP application protocol 203. HP, a 
PDES Inc. member in the U.S.A. and an ESPRIT CADEX 
member in Europe, contributed to the 3D geometric design 
definition of AP203 in a joint effort of PDES Inc. and 
CADEX. The AP203 3D geometries cover solid models, sur- 
face models, and wireframe models and are shared by other 
application protocols, thereby promoting interoperability 
between different application areas. 

HP has also been actively supporting the U.S. initiative to 
define a good-quality standard for associative drafting docu- 
mentation in STEP. Associative drafting, covered by AP202, is 
considered an integral portion of the product data for con- 
tractual, archival, and manufacturing reasons. For example, 
government contracts and ISO 9000 require that product 
data be thoroughly documented. This includes engineering 
drawing data of a produa in addition to the 3D produa data 
and the configuration data. Electronic design and printed 
circuit board design data are also covered in STEP. 

PRODEX. In 1992 participants in the ESPRIT CADEX ptojea 
demonstrated publicly the first B-Rep solid model transfer 
via STEP for mechanical parts in Europe. To develop this 
new technology the PRODEX project was founded in 1992 
with the goal of developing STEP data exchange for CAD 
design, finite element analysis, and robot simulation systems. 
Twelve European companies joined the project. So far, the 
project's achievements include die definition of a STEP 
implementation architecture, the development of a STEP 
toolkit, and the development of STEP preprocessors and 
postprocessors. 

Product data exchange between the different vendors is on- 
going and shows very promising results for CAD-to-CAD 
data exchange, CAD-to-flnite-element-system exchange, and 
CAD-to-robot-slmulation-system exchange. The STEP stan- 
dard has been further fostered by a joint effort with the 
ProSTEP projea to develop AP214, in cooperation with the 
U.S., European, and Japanese automotive industries. 

ProSTEP. ProSTEP is an automotive industry initiative for a 
highway-like STEP product model data exchange. In 1992 
the German companies Bosch, BMW, Mercedes-Benz, Opel 
(GM), Volkswagen, and Siemens launched an initiative to 
bring the major CAD vendors together with the goal of im- 
plementing the first harmonized set of STEP product data 
exchange processors (product dara translators) for indus- 
trial use in the automotive industry. TTie approach taken was 
to compile the user requirements, to build on the results and 
experiences of the ESPRIT CADEX projea, and to launch at 
the ISO level a STEP application protocol, AP214, which 
covers the core data for automotive mechanical design. 
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The following CAD/CIM systems are involved in the projea 
and have STEP data exchange processors either available or 
under development: Alias, AutoCAD, CADDS/CV-Core, 
CATIA, EUCUD3, HP PE/SolidDesigner, EMS-Power Pack, 
I-DEAS Master Series, SIGRAPH STEFIntegrator, SYRKO, 
Tebis, ROBCAD, and others. 

The initial focus in ProSTEP for STEP products is on design 
data exchange for 3D geometry: B-Rep solid models, surface 
models, and wireframe models. For migration from legacy 
systems, wireframe data needs to be supported, at least for 
data import. Communication with applications like numerical 
control (NC) programming systems today typically requires 
surface model data, although in the future more solid model 
data will be used. Initially, the HP emphasis is on bidirec- 
tional produa model exchange (input and output) of 3D 
B-Rep and surface models. 
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STEP Tools A 

In STEP implementation projects, standardization has been 
extended beyond the product data to the STEP implementa- 
tion tools. The CADEX, PDES Inc, PRODEX, and ProSTEP 
projects have all taken this approach. 

A standardized STEP tool architecture provides the following 
benefits. These include shareability of tools between differ- 
ent implementors, shortened development time for STEP 
processor implementations (software development produc- 
tivity gain), increased likelihood of compadbility between 
STEP implementations (differences in STEP definition inter- 
pretations are minimized), parallel development of tools 
(concurrent engineering), extendability of tools to track new 
standardization trends, increased flexibility (new STEP 
models requite fewer code changes), and centralized 
maintenance of tools. 

Fig. 13 shows the PRODEX STEP tools architecture. The 
functional blocks of a STEP toolkit or STEP development set 
are: 

• STEP Standard Data Access Interface (SDAT), 

• STEP Express compiler 



The main interface to the STEP application is the STEP Stan- 
dard Data Access Interface, which provides a computer pro- 
gramming language for dynamic access to the STEP data. 
Application-specific mapping and conversions are imple- 
mented on top of this interface. 

The Express compiler conveys the produa data descriptions 
contained in an Express schema (the metadata of the data 
model) to the toolkit. It contains an Express file reader and 
compiles the file contents to the internal representation of 
the data model. The SDAI is the recipient of the product 
data metamodel and uses the metamodel as a reference for 
the product instance data, wfilch is imported through the 
STEP file scanner/parser. 

The STEP file scanner/parser reads (scans and parses) the 
STEP instance data contained in a STEP data file and uses 
the currendy valid metamodel for checking the syntax of the 
imported instance data. 

The STEP file formatter formats the data to a part-21-confor- 
mant STEP file which is read from the SDAI by using the 
current valid metadata (e.g., a specific application protocol 
such as AP203). 

The STEP data checker is a validation tool that checks the 
instance data currendy in the SDAI based on the correspond- 
ing metadata model, which is also contained in the SDAI. 
The checking covers consistency checks like references 
between entities (e.g., existence dependency), and rule 
checking, which is covered in the metamodel. The checking 
is optionally applicable to the data in the SDAI. It is very 
helpful during die development of processors, for checking 
new metadata models, or for checking the first data 
imported from a new system. 
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The STEP conversion tool is a pool of conversion functions 
(a library) that includes all kinds of geometrical, topological, 
and other model conversions. The focus is on geometrical 
conversions which are heavily used for data exchange be- 
tween systems with different geometric modeling concepts. 
For example, one CAD system might use rational polynomial 
representations for its inherent geometric represenmtion of 
curves and surfaces (e.g., NURBS, nonuniform rational B- 
splines), while the other might use nonrational representa- 
tions (e.g., NUBS). In this case an approximation to the 
nonrational representation has to be applied, at the price of 
increasing the amount of data. For another example, a sur- 
face modeling system might export trimmed surface data 
with curve representations in 2D parameter space, whereas 
the receiving system might handle only 3D .space curves. In 
this case the 2D parameter curves have to be evaluated and 
converted to 3D trimming curves in 3D space. 

By using a STEP toolkit the requirements for the implemen- 
tation of a STEP processor might be reduced to just the na- 
tive data interface to the STEP tools, which consists of the 
data output to the SDAI (for the STEP preprocessor) and the 
data imported from the SDAI (for the STEP postprocessor). 

The main task in linking a CAD system to the toolkit consists 
of defining and implementing the mapping between the 
system internal representation and the standardized entity 
representation in the schema of the standard (e.g., an appli- 
cation protocol). 

HP PE/SolidDesigner STEP Implementatiofi 

The target application protocols for HP PE/SoIidDesigner are 
initially AP203 and AP214, in which both solid and surfece 
models are supported. In addition to the HP PE/SolidDe- 
signer internal data models, the solid and surface models of 
other CAD systems are of major interest. With the introduc- 
tion of STEP, B-Rep .solid model data exchange comes into 
industrial use, representing a new technology shift. HP PE/ 
SolidDesigner has its focus on solid models and is best 
suited for STEP-based bidirectional solid model exchange. 
However, surface models are also supported. 

In addition to the geometric specifications, product informa- 
tion and configuration are covered in the implei 



this article, the geometric and topological mappings are dis- 
cussed. The assembly, product structure, and administration 
mappings are not covered. 

STEP Preprocessor (STEP Output) 

The preprocessor exports the HP PE/SoUdDesigner model 
data in a STEP file. The preprocessor takes care of the map- 
ping of the HP PE/SolidDesigner model to the STEP model. 

The internal geometrical and topological model of HP 
PE/SolidDesigner is in many respects similar to the STEP 
resources of part 42 of the STEP standard. Hence the map- 
ping is often straightforward. On the other hand, there are 
data structure elements that are not mapped to the STEP 

HP PE/SolidDesigner uses the following geometric 3D 
elements: 

» Analytics: 3D surfaces such as planes, cones, cylinders, 
spheres, and toruses, and 3D curves such as lines, arcs, 
circles, and B-splines 

» Nonanalytics: typically 3D elements such as B-spline curves 
and surfaces, and linear and rotational swept surfaces. 

The topology used for the exchange of solid models is based 
on the manifold topology of STEP part 42. The elements 
used are manifold solid boundary representations, closed 
shells, faces, loops, edges, and vertices. The link between 
the topology and the geometry is given by references fiiom 
faces to surfaces and fmm edges to curves. The geometrical 
points are referenced by vertices. 

The HP PE/SolidDesigner STEP surface models ate also 
based on topological representations. Special elements are 
used for surface models, such as shell-based surface models 
and closed and open shells. The other underlying topolocial 
elements are the same as in the solid models. The geometric 
representations of the surfaces are typically the same as in 
the solid model representations. 

STEP Postprocessor (STEP Input) 

The HP PE/SolidDesigner postprocessor supports the import 
of B-Rep solid models and surface models along with the 
necessary produa structure data. The postprocessor is 




Fig. 14. Data exchange cycles 
between different CAD systems, 
including robot simulation sys- 
tems, in the ProSTEP project. 
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Fig. 15. Golf club solid B-Rep model imported into HP PE/Solid- 
Designer from CATIA (CAP-Debis). 



capable of covering at least the functionality of the prepro- 
cessor so that it is possible to store and retrieve HP PE/Solid- 
Designer data in a STEP file representation (this is called the 
^lort cycle test). 

The STEP postprocessor imports STEP files bom other sys- 
tems based on specifically supported application protocols. 
Postproce.ssing is one of the most difficult tasks in data ex- 
change, especially when the data imported comes from a 
system that is very different bom the receiving system. Po- 
tential problems arise in postprocessing if the sending and 
receiving systems have different accuracies, use different 
modeling techniques to generate the data, have different or 
missing surface connectivity, use different algorithms or cri- 
teria to determine surface intersections or connectivity, or 
use different model representations for similar model charac- 
teristics. 

When surface models are imported, it cannot be guaranteed 
that they can be migrated to solid models even with user 
interaction. However, in special cases imported surface 




Fig. 16. Clamp solid B-Rep model imported from Unigraphics II 
(EDS). 
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models can be migrated to solid models without problems. 
In many cases imported surfaces provide boundary condi- 
tions for the solid model. In most cases the data can be used 
as reference geometry to check interference or provide di- 
mensions for the solid models. For example, an imported 
surface set might represent the surrounding boundary geom- 
etry within which the final mechanical part has to fit without 
interference. 

Importing surface models into HP PE/SolidDesigner is con- 
sidered important and critical since many other CAD sys- 
tems, especially legacy systems, often support only surfaces 
or wirefame models, not solid models. Therefore, the post- 
proces.sing of STEP surface models needs to cover a broader 
scope than the preprocessing. Sometimes, different surface 
representations are used in different application protocols, 
such as AP203 and AP214. Hence, different external repre- 
sentations may need to be mapped to one internal represen- 
tation in HP PE/SolidDesigner. 

In the initial implementation of the HP PE/SolidDesigner 
postprocessor, topology bounded surface models are sup- 
ported. These provide the most sophisticated description of 
the connectivity of the individual surfaces used in a solid 
model. Geometrically bounded .surface models are sup- 
ported as a second priority. 

The Accuracy Problem 

When importing CAD data from other systems the accuracy 
of the data plays a key role and determines whether a co- 
herent and consistent CAD model can be regenerated to 
represent the same kind of model in the receiving system. 

Let's define the term accuracy. There are different accuracy 
or resolution values that must be considered in geometric 
modeling and CAD systems. For 3D space, a minimum linear 
distance value (a length resolution value) can be defined, 
which is the absolute distance between two geometric points 
that are considered to coincide in the CAD internal algo- 
rithms; this represents the zero distance. We'll call this value 

* Often, CAD surface models are not conslstert because the generating system lacks checking 
mechanisms or does not track connectiviiy. Very oftea consistency and accuracy are the 
responslbaitles of the user of the system rather than under system control. 




Fig. 17. Wheel solid model imported from SIGRAPH-3D (Siemens- 
Nixdorf). 



Fig. 18. B-Rep model imported from Unigraphics (EDS). 

the linear accuracy . A typical value could be 10 ~^ mm 
which is highly accurate for many mechanical design appli- 
cations. A similar value can be specified for ar^^lar accu- 
racy, parametric accuracy, and so on. The discussion here Ls 
limited to linear accuracy. 

If the sending system uses a higher linear accuracy (more 
precise data) than the receiving system, distina geometric 
points will be detected to coincide in the receiving system. 
Hiis might result in a change in the topology (which might 
cause further inconsistencies) or the geometry. If the sending 
system uses a lower linear accuracy (less precise data), the 
receiving system might complain that the topology is not 
correa or the geometry and the topology are inconsistent. 

To prevent or at least minimize these kinds of accuracy 
problems it should be possible to adjust the accuracy in the 
receiving system to the accuracy values of the data to be 
imported. For example, if the sending system uses a differ- 
ent accuracy for the model generation process, say a linear 
accuracy of 10"^ mm, then the receiving system should 
adjust its internal algorithms to the same accuracy. 

Experience with HP PE/SolidDesigner has shown that this 
kind of adjustable accuracy helps regenerate CAD models 
that were generated in different systems with different accura- 
cies. Also, for data models composed of components witli 
different accuracies, the components can be brought together 
on the assembly level to form a complete product model. 

In the STEP implementation of AP214 an adjustable linear 
accuracy value is conveyed in the STEP file to tell the 
receiving system the appropriate accuracy value for post- 
processing. 

User Features 

The user can select via the HP PE/SolidDesigner graphical 
user interface the objects (e.g., several B-Rep bodies) to put 
into a STEP file. For example, the user decides whether to 
send the data in a B-Rep solid model or a surface model 
representation. The user can choose some configuration 
parameters that help tailor the model data set for best com- 
munication to a specific target application. However, all data 
must comply with the STEP standard. 

When importing (postprocessing) a STEP file the user can 
define some parameters that ease the processing of data. For 
example, the user might set the accuracy value before 



importing a data set that was designed with a specified 
accuracy, or might choose to convert the imported data to a 
different representation. 

STEP Model Exchange Trials 

Various STEP file exchanges have been performed within the 
last 12 months, not always with satisfying results. This has 
resulted in more development work by the exchange part- 
ners. This process of harmonizing the STEP preprocessors 
and postprocessors of different CAD vendors Is considered 
to be of vital importance for the acceptance of the STEP 
standard and its application protocols. Within the ProSTEP 
project this process has worked particulariy well. Other work 
has been done with, for example, AP203 implementors to- 
gether with PDES Inc. 

At this time, solid model data exchange can be said to be 
working very well, especially compared with what was 
possible with existing standards. STEP-based surface model 
exchange has also reached a level that was not possible with 
existing standards like IGES or VDA-FS, especially with 
respect to topological coherence, which is easily conveyed 
with STEP between many CAD systems. Of course, the wide 
variety of surface models, with the resulting accuracy and 
connectivity problems, will need to be addressed by the 
different CAD system vendors to optimize data transfer via 
STEP In the meantime, STEP file exchange has matured to 
the point where STEP products are offered by various CAD 
vendors and system integrators. 

Within the ProSTEP project one of the broadest ranges of 
STEP-based data exchange trials have been performed 
between HP PE/SolidDesigner and odier CAD systems 
(see Fig. 14). Solid model industrial part data has been 
exchanged, for example, with CATIA (CAP Debis and 
Dassault/IBM), Unigraphics II (EDS), SIGRAPH Design and 
STEP Viewer (Siemens-Nixdorf), and others. Some of the 
successful results are shown in Figs. 15, 16, 17, and 18. 
Surface model industiial part data has been exchanged with 
CATIA, EUCUD, SYRKO (IMercedes-Benz corporate design 
system), and others. Some of the successful results are 
shown in Figs. 19 and 20. 

Next STEPS 

Future releases of the STEP standard covering product data 
categories such as materials, tolerances, form features, manu- 
facturing process data, and others are expected in the next 
few months. The expected release of AP202, associative 
drafting, will allow documentation of the product data in 
engineering drawings. Work is ongoing towards the parame- 
terization of product features, which needs further develop- 
ment in the STEP standard. 

The expected finalization of AP214 will make it possible to 
convey the product data categories in STEP files and will 
help to reduce design and manufacturing development 
cycles for simple as well as complex products. This process 
will be supported by further extensive use of data communi- 
cation networks in the various countries. The migration from 
existing standards is aided by several product offerings of 
IGES-to-STEP and VDA-FS-to-STEP data conved:ere. 

The STEP implementation technology based on the STEP 
Standard Data Access Interface will be broadened and used 



October 1995 Hewlett- Packard Journal 49 




Fig. 19. Surface model imported from SYRKO (Mercedes-Benz Hg. 20. Headlight reflector surface model imported from EUCLID 

corporate design system). (Malta Datavision). 

in database access implementations to allow concurrent ac- 
cess by product design and manufacturing development. 

However, for industrial use, die database technology and the 
STEP data access technology need to be extended and inte- 
grated. This process is expected to take several years. 
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